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ABSTRACT 


Chapter I. Historical Review of Studies of Geomagnetic Disturbances 


An introductory review of recent problems concerning geomagnetic disturbances is 
given. 


Chapter IJ. Progressive Change in the Average Current- 
System of Bays 


The end-point of the horizontal disturbing force of geomagnetic bays in the middle lati- 
tudes describes a loop during the whole course of bays. Statistical examination of the loop gives 
rather systematic characteristics concerning its sense of rotation, namely it rotates clockwise 
in the forenoon and counter-clockwise in the afternoon. ‘The above characteristic of the 
disturbing force of bays can be interpreted by a systematic progressive change in the 
average current-system of bays, namely the positive bay area becomes wider and the 
negative bay area contracts during the development and decay processes of bays in the 
middle latitudes. 


Chapter III. Dynamo Theory for Sp-Field and Bay Disturbance 


It: has been shown that an electric current-system similar to the idealized current- 
system of Sp-field can be produced simultaneously with Sg current-system by means of 
dynamo theory when the electrical conductivity of the auroral zone ionosphere becomes 
higher than that of the adjacent regions. However, the calculated current-system is far 
weak in its current intensity and the phase of the calculated current-system is almost 
reversed compared with the observed one, so far as we concern the dynamo-action in a 
single layer. In order to get rid of the above difficulties, double layer hypothesis is proposed, 
in which the presence of two concentric conducting layers in the upper atmosphere is as- 
sumed and the air motion in both layers is opposite with each other. 

A further extension of the dynamo theoretical treatment in the case of inhomogeneous 
distribution of the electrical conductivity is given. If the electrical conductivity along the 
auroral zone assumes large in the dark hemisphere than in the sunlit one, the auroral zone 
current in the resultant Sp current-system becomes more intense in the dark hemisphere 
than in the sunlit side, being in good agreement with observed results. 

It is also shown that a model experiment can give the distribution of stream lines over 
a spherical surface corresponding to various complicated distribution of electrical conduc- 
tivity along the auroral zone. Some experimental results are illustrated. 


Chapter IV. Development and Decay Processes of 
Individual Geomagnetic Bays 


Four examples of geomagnetic bays were examined in detail with the aid of world-wide 
magnetogram copies, especially in regard to the progressive change in the current-system 
during the course of bays. It is seen that the disturbing force of bays in the middle and 
low latitude stations takes its maximum value nearly at the same time when the westward 
auroral zone current in the dark hemisphere reaches its maximum, and the eastward auroral 
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zone current in the sunlit hemisphere seems to become most intense a few tens of minutes 
later than the westward one does. The above characteristic of the auroral zone current 
will result in the general tendency of the broadening of the positive bay area during the 
course of bays in the middle latitudes mentioned in Chapter Ve 


Chapter V. Constitution of Polar Magnetic Storms 


Successive instantaneous aspects of the world-wide distribution of polar magnetic storms 
were traced with the aid of a number of magnetogram copies. It is concluded that the 
polar magnetic storms are composed of a number of elementary disturbances which take 
place intermittently or successively with duration from a few tens of minutes to about two 
hours. Sometimes, two or more of such elementary disturbances occur simultaneously. The 
elementary disturbance can be detected in every stage of magnetic storms. In the analysis 
of the geomagnetic variation by means of hourly values, the elementary disturbances are 
scarcely detectable. ; 

Although the average of all these elementary disturbances is of the form of Sp-field, 
the equivalent current-system corresponding to each elementary disturbance is generally 
much simpler, being approximated by the current-system produced by an electric doublet 
situated along the auroral zone. ; 

This fact will suggest that each elementary disturbance is caused by the production of 
e.m.f. in a rather small area in the auroral zone, and it is proved that the elementary 
disturbance can be interpreted consistently by assuming the production of small highly 
ionized area in the upper atmosphere over the auroral zone. 


Chapter VI. Relation between Sp-Field and D,:;-Field in the 
Middle and Low Latitudes 


Some examples of typical magnetic storms are illustrated, in which a fairly distinct 
development of Sp-field immediately after the sudden commencement of magnetic storms is 
noticed. The corpuscular stream from the sun can be considered to arrive at the upper 
atmosphere over the auroral zone almost simultaneously with the sudden commencement of 
magnetic storms. 

Geomagnetic bays are considered to be caused by a similar mechanism for polar 
magnetic storms, since the presence of the hypothetical equatorial current-ring at the time 
of bays is suggested. 


Chapter VII. Geomagnetic Variations and the Ionosphere 


Some problems concerning the electrical conductivity and air motion in the upper 
atmosphere, the height of the electric current-system responsible for geomagnetic variations, 
ionospheric perturbation accompanying geomagnetic disturbance, and the ionization power of 
impinging extraterrestrial charged particles are surveyed. 


Chapter VIII. Conclusion 


Concluding remarks throughout the present study are summarized, and the physical 
mechanism of polar magnetic storms is discussed. 


Chapter I. Historical Review of Studies of 
Geomagnetic Disturbances 


Since regular magnetic registration by photography was first introduced in 
1847, at the Royal Observatory, Greenwich, England, continuous recording of 
geomagnetic elements has been carried out at various localities over the earth. 
By means of the world-wide data, the spherical harmonic analysis has been 
tried by many investigators for the earth’s main magnetic field as well as its 
transient magnetic variations. The transient magnetic variations include the 
solar daily variation (denoted by S), the lunar daily variation (Z) and the 
disturbance field (D), these three variations comprising the chief changes in 
the earth’s magnetic field together with the secular variation. 

The spherical harmonic analysis of the solar daily variation, S, was carried 
out by many investigators. S. Chapman, whose analysis was more detailed 
than the earlier analyses, summarized and critically reviewed the works of 
his predecessors (See 1 of References at end of the paper). He also analysed 
the lunar daily variation, L [1]. The result of the spherical harmonic analyses 
of S and L shows that the origin of these transient phenomena must exist 
mainly in the exterior part of the earth, and a part under the earth’s surface. 
The ratio of the amount of the external origin part of S or Z to that of the 
internal one is roughly 2 to 1. Since there is no magnetic matter in the earth’s 
atmosphere, the equivalent current systems flowing in the upper atmosphere 
were introduced as the external origin parts of S and L, without any knowledge 
on the ionosphere at that time. Recent observation with rocket flight gave an 
evidence for the presence of electric current flow in the ionospheric region [2]. 
On the other hand, the internal origin parts of S and L were proved to be 
caused by the electromagnetic induction within the earth [3]. Theories of the 
solar daily variations were also proposed, they are dynamo theory, diamagnetic 
theory and drift-current theory. 

Since the ionosphere was found directly by the experiments with radio waves 
in 1925, the studies on physical properties of the upper atmosphere have been 
rapidly developed, and distribution of number density of electrons, ions and 
neutral particles, and other elements with height has been experimentally and 
theoretically estimated by a number of researchers. Their results serve the 
study of external origin part of geomagnetic variations. 

After some criticisms of these theories, the inadequacy of latter two 
theories was shown [4]. The dynamo theory, which attributes the electric 
current in the upper atmosphere to the dynamo action there caused by the 
horizontal movement of the air transverse to the vertical component of the earth’s 
permanent magnetic field, seems to be adequate also in the interpretation of 
the lunar daily variation, because the electromotive force producing the current 
system of L can be considered to be generated by the movement of the air 
due to the tidal motion of the upper atmosphere. Various merits are found 
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in the dynamo theory, for example, longitudinal inequality of Sq-field can be 
interpreted by taking the obliquity of the earth’s magnetic axis with its rota- 
tion axis into account [5,6]. Moreover, the magnetic variations at the time 
of solar flare [7] or solar eclipse can be interpreted from the standpoint of the 
dynamo theory, under the probable assumption of the increase or decrease in 
the electrical conductivity of the upper atmosphere [8,9]. Thus, dynamo 
action in the upper atmosphere seems to be the most important phenomenon 
for the causation of the solar and lunar daily variations of geomagnetic field, 
though we have still little information on the movement of the upper atmos- 
phere even at present. 

The simple and satisfactory definition of the disturbance field, D, observed 
at a point and at an instant is 

D=F—F-S,—L, 
where F is the observed instantaneous value, F the average value of the 
magnetic field over fairly long interval of time (generally more than a month), 
and Sj, is the solar diurnal variation on magnetically quiet days. Since L is 
so small compared with the magnitude of Sq at most stations, it can be 
usually ignored in this equation, which therefore becomes 
D=F—F—S,. 

The disturbance field D is usually small on quiet days, while during periods 
of disturbance its intensity changes rapidly as a function both of position and 
time, or in other word, the magnetic field shows complicated variation at the 
time of magnetic disturbance. The disturbance field is generally larger and 
more complicated in high latitudes than in lower latitudes. 

The studies on the disturbance field have been carried out by averaging 
a number of disturbances, because of the extreme complexity of individual 
disturbances. The disturbance field D can be separated, for convenience, into 
two parts, namely one depending on local time and the other. The former is 
the diurnal part of the D field and called the solar disturbance-daily variation 
and is denoted by Sp. The latter one include the mean world-wide variation 
Dm, which is the zonal component from spherical harmonic analysis, and the 
irregular variation, D;. Namely, 

D=Sp+Dm+D;. 

The solar disturbance-daily variation can be expressed by the overhead cur- 
rent system such as given by Fig. 2 in Chapter II, which is obtained by S. 
Chapman. However, Sp is not simply represented by the above-mentioned 
current system throughout the period of disturbance, but the current intensity 


seems to vary time by time, so that Sp may be written from the statistical 
standpoint in the following form 


Sp=k(T)-Sh(t), 
where Sh) is an idealized disturbance field such as represented by Chapman’s 
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Sp current system, depending upon local time ¢ alone, and k(T) means the 
magnitude of Sp varying its value with universal time 7. According to M. 
Hasegawa’s detailed analysis, the current system of Sp can also be found even 
in the average magnetic field variation on international quiet days, though the 
current intensity is far less than that observed on disturbed days, showing that 
R(T) in the above expression does not vanish even on magnetically quiet days 
[10]. It must be deeply mentioned that the above simple and idealized 
expression will be applied only from the pure statistical standpoint. 

Dm and D; are non-diurnal components of the disturbance field according 
to the above definition. In the average of a number of disturbance fields, the 
irregular variation D; will be averaged out or becomes very small, and Dm 
alone will remain. It has been remarked that the typical magnetic storm 
begins with sudden commencement, noticed by the abrupt increase in the 
horizontal component of magnetic field, and then after some several hours the 
horizontal intensity begins to decrease, and gradual recovery follows after its 
minimum in the middle and low latitudes. The mean world-wide variation of 
the geomagnetic field at the time of magnetic storm is especially called Ds, 
the storm-time variation, its time being reckoned from the commencement of 
magnetic storm. In other world, the average disturbance field at the time of 
typical magnetic storms will be written as 


D=Spt+Dst. 


The above separation of the disturbance field into Sp, Dm and D;, is only for 
convenience as mentioned before, and it must be remembered that Sp, Dm and 
D; in the above definition do not correspond to different physical mechanisms 
respectively. 

The dependency of magnetic activity upon the solar activity have been 
examined and established by various statistical investigations. It will be no 
doubt that the magnetic disturbance is caused by the corpuscles emitted from 
the sun, after the ultraviolet theory of magnetic storm failed in the satisfactory 
interpretation of the disturbance field. Early ideas by K. Birkeland [11] and 
C. Stormer [12] have a difficulty that the emitted charged particles of one sign 
will be diluted and diverged away during the travel from the sun to the 
earth, and consequently it is concluded that the solar stream, which causes 
geomagnetic disturbances, must be ionized but electrically neutral. The mean 
velocity of the solar corpuscles is considered to be order of 108cm/sec. from 
various results of the statistical relation on the geomagnetic disturbance and 
the solar activity. According to C. Stormer’s result, the emitted solar particle 
was suggested to be Ca*. Many efforts were made to detect the presence of 
approaching Ca*, but no conclusive result was obtained. 

It has been well known that the intense auroral display can be seen at 
the time of magnetic disturbance, and L. Harang found the close parallelism 


298 N. FUKUSHIMA 


between auroral luminosity and geomagnetic activity [13]. L. Vegard found 
the Balmer lines of hydrogen atom in the auroral spectrum in 1939 (14). 
Recent results by A.B. Meinel [15] and C. W. Gartlein [16] show a strong 
evidence of the presence of protons incoming into the upper atmosphere over 
the auroral zone along lines of the earth’s magnetic force, at the time of intense 
auroral display. The velocity of impinging protons at the level of auroral 
height was estimated by M. Sugiura after a rough analysis of their data, and 
it was concluded that the distribution of velocity of protons has its maximum 
around 1500—2000km/sec., and it spreads even over 4000 km/sec. [17]. One 
may safely say that the materials of impinging charged particles are mainly 
protons and electrons at the present stage of investigations [18]. 

On the other hand, three representative theories of magnetic storms were 
presented, in which magnetic disturbances are considered to be caused by the 
neutral and ionized corpuscular stream from the sun; they are theories by 
S. Chapman and V.C.A. Ferraro [19], H. Alfvén [20] and by D.F. Martyn 
[21]. The mechanism of the emission of electrically neutral corpuscular stream 
from the sun is not deait with in the theories of magnetic storms, but has 
been studied by astrophysicists, for example by T. Hatanaka [22] and others. 
The fundamental problem in every theory of magnetic storms is the motion of 
the neutral ionized solar stream of the velocity of order of 108cm/sec. under 
the influence of the earth’s permanent magnetic field. In Alfvén’s theory, the 
influence of the solar magnetic field is also taken into account. 

In every theory, sudden commencement of magnetic storm is interpreted 
as the increase of magnetic energy near the earth’s surface at the approach 
of the neutral and ionized solar stream to the earth, which was first suggested 
in Chapman and Ferraro’s theory. V.C. A. Ferraro showed in his recent work 
a further detailed mathematical treatment on the problem [23]. 

It is proved that there arises a forbidden space around the earth for the 
solar stream, and the stream should flow with a carved hollow in the vicinity 
of the earth [20, 24]. Then the charged particles near the surface of the for- 
bidden space flows to surround the earth forming a circular electric current flow. 
Dsz-field in the main phase of magnetic storm or Dm-field is interpreted by 
the magnetic field of the equatorial current ring, which is situated at the distance 
of several earth radii in the geomagnetic equatorial plane and electric current 
flows westward, though its mechanism of formation is different in respective 
theory. 

The magnetic field of the equatorial current ring is a uniform magnetic 
field from north to south along the geomagnetic axis in the vicinity of the 
earth. Various statistical results concerning Dm-field, for example by B. Cynk 
[25], show that the mean magnitude of the depression of horizontal intensity 
of geomagnetic field at the time of magnetic disturbances over the world is 
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nearly proportional to the cosine of the geomagnetic latitude in the low and 
middle latitudes. Hence, Dm-field in the low and middle latitudes is quite 
naturally interpreted by the presence of the equatorial current ring, because 
the horizontal component of magnetic field of the equatorial current ring is 
also proportional to the cosine of the geomagnetic latitude on the earth’s 
surface. However, intense zonal current flow in the high latitudes especially 
near the so-called auroral zones, such as seen in the current system for Dse- 


field by S. Chapman remains unexplained from the abovementioned view 
point. 


As for the mechanism of the causation of Sp, Chapman and Ferraro made 
no attempt to account for it. Alfvén and Martyn discussed the leakage of 
charged particles from the equatorial current ring. Charged particles expelled 
from the equatorial current ring impinge into the earth’s upper atmosphere 
along lines of the earth’s magnetic force, because of the least resistance along the 
path. The line of magnetic force started from the distance of several earth 
radii in the equatorial plane intersects the earth’s ionospheric region at polar 
distance of about 20°—25°, which agrees with the so-called auroral zone. In 
Alfvén’s theory, positive ions impinge upon the auroral zone ionosphere of the 
sunlit side, while electrons impinge in the opposite side. Then, electric current 
flows along the auroral zone westward in the forenoon side and eastward in 
the afternoon side. Although he does not mention about the electric current 
flow in the polar caps and the equatorial region in the current system for Sp- 
field, which flows in opposite direction to the auroral zone current, he insists 
the adequacy of his theory after comparison of the actual magnetic data and 
the result of his model experiment measuring the magnetic field produced by 
his current system. In short, the magnetic variation in the polar caps and the 
equatorial region is attributed to the special form of the auroral zone in his 


theory. On the other hand, charged particles of positive sign reach the upper 
atmosphere at the outer boundary of the auroral zone in the forenoon side 


and at the inner boundary in the afternoon side, while negative charged 
particles to the inner boundary in the forenoon side and to the outer boundary 
in the afternoon side of the earth after Martyn’s opinion. Then, drift current 
flows under the presence of the vertical magnetic field of the earth, the direc- 
tion of the drift current being westward in the forenoon side and eastward in 
the afternoon side. The electric current flow of counter direction in the other 
regions is interpreted to be caused by the electrostatic field due to accumula- 
tion of positive charge in the dark hemisphere and negative charge in the 
sunlit hemisphere in the auroral zone. 

First attempt to interpret Sp current system by dynamo theory was 
presented by T. Rikitake [26]. He obtained an additional current system 
quite similar as Chapman’s one from the dynamo theoretical treatment, when 
the conductivity over the auroral zone becomes higher than that in the other 
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regions. This result seems to make us convince that the dynamo action is the 
most important process for the interpretation of various magnetic variations. 

Summarizing the above results, it is said that every theory of magnetic 
storms tried to interpret the physical mechanism of the causation of the dis- 
turbance fields of idealized form, Dz and Sp. The separation of the disturbance 
field into Sp, Dm and Dj is only for convenience as stated before, the physical 
view for the composition of the disturbance field is not taken into account. 
However, the detailed study of the individual magnetic disturbances, for 
example by T. Nagata [27] and E.H. Vestine [28], show that the disturbance 
field does not always appear in its idealized form, which is derived from the 
statistical study. For example, total amount of the auroral zone current is 
more intense in the westward flow than in the eastward one at most cases, 
and the mode of the current flow in Sp current system varies slightly time by 
time, after the analyses of the hourly mean aspect of the disturbance field 
with the aid of the world-wide hourly values. 

Polar part of Dm, which manifests itself as fairly intense zonal current 
flow near the auroral zones and in the polar caps, is not explained in every 
theory, while equatorial part of Dm is interpreted by the magnetic field of the 
equatorial current ring. If we combine the idealized Sp current system and 
the polar part of Dm, the combined current system reveals that the westward 
auroral zone current is much intense than the westward one, and the combined 
current system shows better agreement with the observed phenomenon. 

A reasonable separation of the disturbance field from the physical stand- 
point will be as follows R 

D=Dy4+Dr, 
where Dy means the world-wide variation of the disturbance field due to the 
equatorial current ring, while Dy includes Sp, polar part of Dm, and D;. Sp, 
polar part of Dm, and most part of D; are considered to be originated mainly 
in the upper atmosphere of the earth. In other word, Dy is of an extrater- 
restrial origin, while Dy a terrestrial origin. Strictly speaking, however, 
although the equatorial part of Dm is considered to be of an extraterrestrial 
origin, its solstitial inequality indicated by B. Cynk [25] should be attributed 
to a terrestrial origin. Terrestrial origin part of the disturbance field includes, 
of course, the phenomena originated in the earth’s upper atmosphere. Terrestri- 
al origin part Dy will be divided into equatorial origin part Dg and polar 
origin part Dp, according as the location of the origin of the disturbance 
is situated in the equatorial region or the polar region including the auroral 
zone, namely 

Dy=D;4+ Dp. 
The Sp current system is clearly of polar origin, because the current flow in 
the equatorial region as well as the polar parallel current are considered as 
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the leakage current flow from the intense auroral zone current. Polar part 
of Dm is, needless to say, of polar origin. Then, De includes only a part of 
Di, whose origin locates in lower latitude regions, denoted by Dzg;. Remaining 
part of D;, origin of which is situated in high latitudes is denoted by Dp; and 
belonged to Dp. Hence, 

Dr=D;;; 

Dp=Sp+(polar part of Dm)+Dp;. 
Since we know no systematic pattern of Dz;, Dz; will be eliminated when 
many storms are superposed and averaged. Dp; will be also averaged. out in 
the statistical study. Then, neglecting the irregular features of geomagnetic 
variations, the disturbance field D will be expressed as 

D=D) wtD Pi 

where the world-wide variation Dy is of an extraterrestrial origin, and Dp 
includes disturbances originated in the polar regions in the vicinity of the 
earth’s surface. 

Now, it must be again mentioned that the most part of preceding discus- 
sion deals mainly with the average characteristics of magnetic disturbances. 
Polar magnetic storms show very complicated variation of geomagnetic field 
in high latitudes in general. Because of the extreme complexity of individual 
disturbances, the chief works on the world-wide morphology of polar magnetic 
storms are carried out from the statistical standpoint, as already stated before. 
It is also of importance to consider disturbances individually, in order to 
determine the characteristics of the irregular movements, which are averaged 
out in the statistical examination. Irregular variations of magnetic field, 
which are not a world-wide phenomenon, are considered of course to be 
a terrestrial origin. Some details of individual magnetic storms were ex- 
amined by early researchers, for example by W.G. Adams [29], A. Schmidt 
[30], K. Birkeland and others, before the recent statistical works on the 
world-wide morphology of magnetic storms started. They discussed the 
similarity and difference of the records of the same magnetic storm registered 
at various stations after scrutinizing them. As an example of regional distur- 
bance field, A. Schmidt found a striking peculiarity of the storm curves, which 
is supposed to be interpreted by the movement of localized vortices of dis- 
turbing electric current [30]. K. Birkeland studied polar magnetic disturbances 
of moderate magnitude, and found some kind of polar elementary storms after 
the comprehensive analysis of his expedition results [31], but he did not refer 
to severe magnetic disturbances, such as great magnetic storms. Although 
Birkeland’s interpretation for the physical mechanism of causation of magnetic 
disturbances includes some doubtful points, his monumental work concerning 
the polar magnetic perturbations exhibits various important features of 
elementary disturbances. It is not surprising that little attempt has been 
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made since then to account for the irregular features, because trials for 
explaining the world-wide regular features of magnetic disturbances have been 
in progress. It must be remembered that any theory of magnetic storms is of 
no value if it fails in interpreting the mechanism of any disturbance of simple 
type, such as polar elementary storm, even if the idealized current systems 
for Sp and Dm seem to be well explained by the theory. 

In the present stage of investigation, the morphology of magnetic distur- 
bances must be re-examined with physical picture for the mechanism of 
disturbance fields. Systematic examination of dynamical behaviour of the 
disturbance field, namely the spatial distribution of geomagnetic variations 
and its progressive change over the world is one of the most important pro- 
blems for criticizing the theories of magnetic storms. Hence in the followings, 
progressive change in the equivalent current system for the average of a 
number of bays and for individual cases is analysed at first, because geo- 
magnetic bays are the simplest polar disturbance. Theoretical and experiment- 
al interpretations are also given. In the next place, the constitution and the 
development process of polar disturbance field are examined in detail with the 
aid of world-wide magnetogram copies. These results will give many useful 
suggestions for the physical mechanism of polar magnetic storms. Probable 
mechanism of polar magnetic disturbances is also discussed. 
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Chapter II. Progressive Change in the Average 
Current System of Bays 


§1. Geomagnetic bays and Sp-field 


The average characteristics of the bay disturbances in geomagnetic field 
have been investigated by many authorities, for example by L. Steiner [32], 
H. Hatakeyama [33], A.G. McNish [34], H.C. Silsbee and E.H. Vestine [35], 
J. Ma. Princep Curto [36] and others. Some of them introduced the overhead 
current system equivalent to the distribution of the disturbing force of bays 
with the aid of the world-wide data [33, 35]. The equivalent current system 
for bays obtained by H.C. Silsbee and E.H. Vestine is given in Fig. 1, while 


SCALES OF FORCE IN Ganmas 


vt 600 


24% GMT, 12 MONTHS, (932 -IF 


Fig. 1. Mean 3-hour disturbance-vectors of magnetic bays centred 
at 21h, 24h and 03h G.M.T. (tentatively corrected for induced cur- 
rents) and corresponding average current system of height 150km., 
viewed from above geomagnetic north pole. (after H.C. Silsbee and 
E. H. Vestine) 


the idealized current system for the disturbance daily variation, Sp, by S. 
Chapman [37] is shown in Fig. 2 for comparison. Detailed examinations of 
Sp-field by T. Nagata [27], M. Hasegawa [10], E. H. Vestine [28, 38], and L. 
Harang [39] show that the idealized current system by S. Chapman should be 
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. Fig. 2. Idealized current system in amp. for Sp-field, magnetic and geographic axes 
being assumed coincident ; (A) view from equator at noon meridian, (B) view from north 
pole. (after S. Chapman) 


advanced in its phase by a few hours, in order to be in better agreement with 
the observed geomagnetic variation in the high latitudes. Then, the current 
system for Sp-field shows quite resemblance to that for bays, and there seems 
to be no essential difference between these two current systems. Geomagnetic 
bays might be considered to be a temporary appearance of Sp-field within a 
few hours. In other word, bays are one of the most simple type of polar 
magnetic disturbances, and comprehensive studies on geomagnetic bays are 
desirable for investigating the nature of polar magnetic disturbances. 

The current system in Fig. 1 shows the average distribution of the overhead 
current flow at the maximum stage of bays, occurrence times of which are 
centred at 21h, 24h and 3h G.M.T. It will be not unreasonable in the sight of 
magnetic records, to consider that the development and decay of disturbing 
force of bays will correspond to the intensification and diminution of the 
current intensity of the current system such as in Fig. 1, as a first approxima- 
tion. It has been noticed, however, that the end-point of magnetic field vector 
during a bay disturbance generally describes a loop, and horizontal projections 
of the loops have been obtained at some observatories, as will be mentioned 
later. The fact that the end of magnetic field vector describes a loop during 
the course of a bay disturbance, indicates that the development and decay of 
disturbing force of bays are not only caused by the mere intensification and 
diminution of the current intensity of the current system, but are accompanied 
by some kind of movement or deformation of the current system itself. In 
this chapter, the progressive change in the current system of the bay dis- 
turbance mainly in the middle latitudes during its whole course was examined 


. 
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through a statistical study of distribution of disturbing force of bays at some 
observatories over the world. The distribution of magnetic observatories, 
whose data were dealt with in this study, is given in Fig. 3. 
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Fig. 3. Distribution of magnetic observatories, whos2 data were dealt with in the 
study of characteristics of geomagnetic bays. (San Juan must be dropped in the figure) 
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§2. Progressive change in the average disturbing force of bays ata station 
during the whole course of bays 


2.1. Disturbing forces observed at Toyohara 


The vector-diagram; illustrating the change in the horizontal disturbing 
forces during the whole course of bays observed at Toyohara (46758/N, 
142°45/E), Saghalien, are given in Fig. 4, which were obtained by H. Hatakeyama 
[33] with the data from August 1932 through August 1935. The method used 
for the analysis was as follows. The whole duration of a bay disturbance was 
divided into six parts; the time of its beginning being denoted by 7% and that 


of its end by fs, while the middle of the whole duration by #3. The disturbing 
ee eae 


ce 
AIF RELL 


0~2h 2~4h 4~6h 6~8h 8~10h 10~12h 12~14h 11~16h 16~18h 18~20h 20~%2h 22~24h 


Fig. 4. Mean vector-diagrams of the horizontal disturbing forces of bays in every 
two-hour time interval in local time, observed at Toyohara. (after H. Hatakeyama) 
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force of bays was defined approximately to be the deviation from the straight 
line connecting the values at f and ¢; on magnetograms for each element. 
Whole selected bays were divided into twelve groups according as #;s belonged 
to 0—2h, 2—4h, ----: , 22—24h in local time. Numbers of bays used for deriving 
the mean diagrams in Fig. 4 were 52, 30, 16, 20, 14, 23, 17, 20, 37, 52, 68 and 59 
for every two-hour interval mentioned above, respectively. The rectangular 
axes in the vector-diagrams are referred to the magnetic N-S and E-W direc- 
tions. Regular diurnal change in the direction of the major axis of loops can 
clearly be seen, and of course the current system such as shown in Fig. 1 well 
expresses the above regular diurnal change. The vector-diagram of the dis- 
turbing force forms a loop during the course of bay, and the sense of its rota- 
tion depends upon the local time, namely clockwise in the forenoon and counter- 
clockwise in the afternoon. If we assume a bay disturbance as the simple 
intensification and diminution of the current system without any deformation 
of its mode such as given in Fig. 1, no clockwise rotation of disturbing force 
‘can be expected, even if the rotation of the earth with respect to the current 
system is taken into account. The same tendency of loops of the vector- 
diagrams has been noticed since earlier research by L. Steiner, who examined 
the magnetic records registered at O/Gyalla (47°52’N, 18°11/E), Hungary, during 
1906—1917 [32]. 

Now we suppose a bay disturbance, whose vector-diagrams of horizontal 
disturbing forces are distributed as given by Fig. 4, sofar as it were observed 
at each point of different local time on the same latitude circle as that of 
Toyohara. We shall call the bay thus defined “the average bay” in the 
followings. In the writer’s previous paper, the progressive change in the 
current system of the average bay was discussed [40]. According to the con- 
clusion in the paper, the positive bay area becomes wider and consequently the 
negative bay area becomes narrower with the progress of time during the 
whole course of the average bay in the middle latitudes, where the 
positive or negative bay area means the region where positive or negative bay 
is observed. A kinematical interpretation of the relation between the mode of 
progressive change in the current system of the average bay and the above- 
mentioned rotation of disturbing force vectors, which is subjected to the local 
time, was also discussed in the paper. Outline of the result will be given in Rabe 


2.2. Disturbing forces observed at Cheltenham and Tucson 


For the purpose of examining whether the above-mentioned characteristics 
of the looping phenomenon of disturbing forces of bays can be seen everywhere 
in the middle latitude stations, the records of American stations, Cheltenham 
(38°44’N, 76°50/W) and Tucson (32°15/N, 110°50’W) were analysed in collabora- 
tion with H. Ono and others [41]. “Magnetograms” published by U.S. Coast 
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and Geodetic Survey were available in this study. The method of analysis is 
quite the same with H. Hatakeyama’s one, with only difference in classification 
by the local time in a day, namely a day is divided into eight groups, 0—3h, 
3—6h, --- , 21—24h in local time in the present case. Numbers of selected 
bays in each three-hour time intervals in the period of 1947 and 1948 are 41, 
42, 18, 12, 15, 26, 29 and 49 respectively. The disturbing force of bays is 
defined to be the deviation from the straight line connecting the values at 
beginning and end of bays on magnetograms for each element, as already 
adopted in the case of bays at Toyohara. The loop of average disturbing 
forces shows the weighted mean of modes of the change in the disturbance 
field vector of bays with time occurred in each time interval, where the magni- 
tude of individual bays itself is adopted as the weight for deriving the average 
mode. As will be seen in Fig. 5, magnetic north component of the horizontal 
disturbing force is denoted by 4H, while magnetic east component is expressed 
by H)4D, where Hy is the mean horizontal intensity of geomagnetic field, and 
4D is the change in declination (eastward positive) in radian. Let geographic 
north and east components of disturbance force denote 4X and AY respec- 
tively, and geomagnetic north and east components be 4Xm and 4Ym, the 
following relations hold; 
4X=4H cos D—(H)4D) sin D, 
AY=4H sin D+(A)4D) cos D, 
AXm=AH cos (D—W)—H)4D sin(D—-wW), 
4Ym=4H sin (D—W)+H)4D cos (D—p), 


where Y denotes the angle formed by the 
great circle joining the station and the 
geomagnetic pole with the geographical 
meridian of the station, as will be seen in 
Fig. 5. Declination at Cheltenham and 
Tucson during the period was —7°.1 and 
13°.6 respectively, while Y is 2°.4 and 10°.1 
for each station, so that D—y is —9°.5 for 
Cheltenham and 3°.5 for Tucson. Since D, 
wv, D—w are not so large in both stations, 
the geomagnetic or geographic north and 
east components do not much differ from 
the magnetic north and east components. 
The vector-diagrams illustrating the average mode of change in the 
disturbing force of bays during the course of bays are shown in Fig. 6 for the 
two observatories, where the rectangular axes in the figure are referred to the 
magnetic north-south and east-west. Diurnal change in the direction of the 
major axis of loops in the vector-diagrams shows quite the same tendency as 
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Fig. 6. Vector-diagrams illustrating the change in the mean horizontal disturbing force 
of geomagnetic bays in every three-hour interval observed at Cheltenham and Tucson. 


those reported hitherto by several researchers. The sense of the rotation of 
loops in Fig. 6 is clockwise in the forenoon and counter-clockwise in the after- 
noon at both observatories, with only an exceptional one for 0—3h at 
Cheltenham. It is seen that the sense of the rotation of loops just mentioned 
above, namely clockwise in the forenoon and counter-clockwise in the afternoon, 
is common to the results at middle latitude stations in Europe, Asia and also 
in America. Therefore, it might be said that ‘the above characteristic of 
the disturbing force of bays will be found everywhere in the middle latitudes. 

From disturbing force vectors at #2 and %, mean angular velocity of the 
rotation of horizontal disturbing forces near the middle of the whole duration 
of bay can be derived for each loop and will be denoted by w. The angular 
velocity is taken positive or negative according as the rotation of the loop 
is counter-clockwise or clockwise. ‘The relation between @ and local time » 
(measured in angle from midnight meridian) is shown in Fig. 7 for the two 
observatories. If w assumes to be expressed by the form 


@=aytc, sin(A+8) , 
the relations between » and X in these two cases are 


@=—3+17 sin (A+150°) °/hour for Cheltenham, 
@=—3-+24 sin (A+169° ) °/hour for Tucson, 


where the calculations were made under the assumption that the duration of 
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the average bay is just two hours. The 
Same expression for the data at Toyohara 
obtained by H. Hatakeyama [33, 40] 
shows 


w=—9+46 sin (A+197°) °/hour, 
as illustrated in Fig. 8. In every case, 
the minimum and maximum values of 


take place in the morning and in the 
A evening respectively. 


Rite Ul abana 
Toyohara 


oo am side pm side x pm. side 
Veo ec | n 
0° 90° 170° 270° 360° 5% ° 5 7 0% 
Fig. 7. Relation between w and Fig. 8. Relation between w and 
local time ?, at Cheltenham and Tucson. local time 2 observed at Toyohara, 


2.3. Disturbing forces observed at Sitka and the progressive change in the 
current system of the average bay disturbance 

The result of the similar analysis of the data at Sitka (57°04’N, 135°20/W), 

Alaska, as well as the simultaneous variation in the vertical component 

(4Z) there during 1946—1948, are illustrated in Fig. 9, where the number of 
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Fig. 9. Vector-diagrams of the mean horizontal disturbing force of bays and 
simultaneous variation in the vertical intensity observed at Sitka. (after N, Fukushima 


and H. Ono) 
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bays adopted in each three-hour interval of local time is 56, 50, 18, 12, 8, 10, 
39 and 80 for 0—3h, 3—6h, ----, 21—24h. We notice here at once that the major 
axis of the loop of horizontal disturbing force does not change its direction 
smoothly counter-clockwise with local time in a day as in the cases for middle 
latitude stations, but it directs almost southward throughout the forenoon, 
northward at 12—18h and. eastward at 18—24h in local time. The geomagnetic 
latitude of Sitka Observatory is 60°.0N. Disturbing forces of bays observed 
at such a station of high geomagnetic latitude as 60°, show distinctly different 
characteristics compared with those in the middle latitude stations. 

Now we suppose “the average bay disturbance” as before. In the average 
bay, the horizontal disturbing force of the bay at a locality on the earth’s 
surface is given by such vectors as given in Figs. 6 and 9, so far as the 
geomagnetic latitude and local time in those figures are taken as the coordi- 
nates of the earth’s surface to indicate the locality. Since geomagnetic 
latitudes of Sitka, Cheltenham and Tucson are 60”.0, 5)°.1 and 40°.4 respectively, 
the vector-diagrams in Figs. 6 and 9 assume to represent the variation of 
horizontal disturbing forces of bays at the positions of various local times on 
60°, 50° and 40° parallels. Equivalent overhead current arrow maps for three 
stages of the average bay are given in Fig. 10, where the current arrows are 
written with reference to geomagnetic coordinates. The three stages are the 
developing stage, the maximum stage and the decaying one, represented by 
the mean stage of ¢; and f2, the state at #;, and the mean of % and i, 
respectively. Stream lines in the maps are so drawn that 2/3 of observed 
north- and east-components are taken as external origin part of respective 
components, and the simultaneous variation in Z-component is also taken into 
account. The stream lines which divide the middle latitudes into positive 
and negative bay areas, are drawn with thick lines in the figure. In the 
middle latitudes, there can be noticed a tendency that the positive bay 
area becomes wider and the current system shifts eastward on the whole with 
the progress of time during the course of the average bay. This characteristic 
of the current system was interpreted kinematically by means of the relation 
between w and A [40], and the interpretation will be briefly given in § 3. 

One can see in Fig. 10 at once that the effect of the auroral zone current 
appears in the disturbing force at Sitka, though the centre of the auroral 
zone current seems to be situated higher than 60° in geomagnetic latitude. The 
overhead current system in the polar region could not be obtained in this 
study. However, a current system such as given in Fig. 11, which illustrates 
the average current distribution for polar magnetic disturbances (See Chapter 
V, §6), seems to harmonize well with the current system of the average bay 
in Fig. 10. Diurnal change in the direction of disturbing force of bays at 
Sitka, which is illustrated in Fig. 9, is strongly subjected to the longitudinal 
asymmetry of the auroral zone current. General aspect of the combined 
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300 

Fig. 10. Progressive change in the over- 
head current arrows and the current system 
of the average bay disturbance. 


current system of Figs. 10 and 11 
shows fairly good agreement with 
the current system by H.C. Silsbee 
and E.H. Vestine in Fig. 1. 

It has been noticed that bays 
occur most frequently at about local 
midnight, and much less in frequency 
in daytime. Hence, there arises a 
possibility that the disturbing force 
at daytime will be rather overesti- 


mated and that at night rather 
underestimated in “the average bay.” 


The mean current intensity should 
be somewhat large in the dark 
hemisphere and weak in the sunlit 
one in practice, compared with the 
current distribution of the average 
bay in Fig. 10. Nevertheless, the 
above conclusions on the progressive 
change in the current system of bay 
disturbance will remain without 
serious alteration. 


2.4, Disturbing forces observed at 
Honolulu 
As for the disturbing force of 


Fig. 11. A mean equivalent current 
system for the disturbance field of polar 
magnetic storms. (after T, Nagata) 
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Fig. 12. Vector-diagrams of the mean horizontal disturbing force of geomagnetic bays 
observed at Honolulu. (after N. Fukushima and H. Ono) 


bays observed at Honolulu (21°19/N, 158°04’W), Hawaii, the result of the 
similar treatment with the data during January 1946 through June 1949 is 
given in Fig. 12. Numbers of bays used for deriving the average disturbing 
force are 53, 20, 5, 38, 20, 40, 33 and 69 for 0—3h, 3—6h, -:---- , 21—24h in local 
time. Geomagnetic latitude is 21°.1 N there. The magnitude of disturbing 
forces at Honolulu is much small compared with those observed in the middle 
latitudes. The disturbing force vector of bays directs almost northward or 
southward, showing that geomagnetic bays can hardly be recognized in the 
record of declination. The rotation of loops of vector-diagrams at Honolulu 
does not always show such a systematic character as detected in the middle 
latitude stations. This may be perhaps due to an ambiguity in defining the 
disturbing force of bays, since the magnitude of disturbing force of bays 
itself is so small there at such a low latitude station. 


§3. A kinematical interpretation of the progressive change in the current 
system of the average bay 

The angular velocity of the rotation of the overhead current arrows on 

a latitude circle, w, is approximately expressed by 
@=a)+c,sin(A+8), >) 

as already mentioned in the preceding paragraph, where \ means longitude or 
local time. Let the angle of direction of the overhead current arrow measured 
counter-clockwise from the north denote by a. Then, 


eda hd)", (2) 
and oo =o0(r)=a+¢ sin(A+8), (3) 


where ¢ denotes time measured from the beginning of bay. Distribution of the 
direction of current arrows with longitude or local time on a latitude circle at 
t=0 will be given by 


aA, 0)=A+aA0+9A), (4) 
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where gA)=yA+277) —a (5) 


and is a constant. The region in which 0<a< 7 corresponds to the negative 
bay area, because the overhead current in such a direction reduces horizontal 
component of geomagnetic field there, while the region in which r<a<2r 
corresponds to the positive bay area. Since a@ can be considered to be a 
monotonously increasing function of A, 


Be =14+9(r)>0. (6) 
In general, @ changes with time during the 
whole process of the average bay disturbance, 
or in other world, if the position of a same 
value of a@ is traced, longitude X of the posi- 
tion should shift with time. In this case, the 
relation between X% and ¢ becomes 


OX ay te, Sin (A+B) 


ors 149) @) 


as will be seen in Fig. 13 from the relations 
(3) and (6). 

The longitudes of two boundary meridians 
between the positive and negative bay areas, 
Ay and Ag (Ay<Azg), are determined by the 
conditions 


Fig. 13. 
sin {Ay+Q+g(A1)}=0, ! (8) y 
sin {Az+Qo+9(Az)}=0. 
If the following conditions 
| ao|<|a sinQ,+A)|, 
(9) 
{ao |<|a sin (A2+8) |, 
are satisfied, then 
Ory / Oz 0 
LOt-|. OF us yD) 


This means that the shifting of two boundaries is opposite in direction, resulting 


in the phenomenon that the positive or negative bay area becomes wider with 


time according as ee is positive or negative. If the condition (9) is not 


satisfied, the mode of the change in the current system is generally rather 
complicated. The expansion of the positive or negative bay area is subjected 


to the condition that a is positive or negative. Since the average 
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Ox 
value of OE 
that the current system has a tendency to rotate on the whole westward with 
the angular velocity a), though the life of the current system is generally 


shorter than a few hours. 


along a latitude circle is —a) when |g/()| «1, it can be said 


Table I. Various elements in the average bay at Cheltenham, Tucson and Toyohara. 


22 sin(21+8) sinQdat+8) cysinQat+8) cisin(22+8) 


Observatory % 6 1 

Cheltenham =3°. 17150?) 5b? S290? ~0.423 0.985 -— 7.2 16.7 
Tucson =3~ 24 169) © 100me.285 ~1.000 0.998 -24.0 24.0 
Toyohara =9>-46 1977 02) =295 -0.485 0.743 -22.3 34.2 


In the case of the present analysis of the data at Cheltenham and Tucson 
as well as those at Toyohara, the result is summarized in Table I. In every 


ONEos we ONpe. F 
case, ae is positive and er is negative, showing the broadening tendency of 


the positive bay area with the progress of time. Negative value of a) at every 
station reveals the eastward rotation of the current system on the whole, 


which was already shown in Fig. 10. 
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Chapter III. Dynamo Theory for Sp-Field and 
Bay Disturbance 


§1. Simple dynamo theory for Sp-field and bay disturbance 


Among theories of the solar daily variation of geomagnetic field on quiet 
days, Sa dynamo theory seems to succeed in interpreting also the. disturbance- 
daily variation, Sp. T. Rikitake has previously shown that a current.system 
of similar mode to the current-system for Sp-field or bay disturbance can be 
derived from the dynamo theory under the assumption that the electrical 


conductivity above the so-called auroral 
zones alone is much high compared with 
that in the other regions [26]. A further 
extended treatment of the dynamo theory was 
carried out by the present writer [42, 43]. In 
the followings, outline of the results is sum- 
marized. 

As will be shown in Fig. 1, the electrically 
conductive layer of radius Ayo, which is con- 
centric with the earth, is divided into five 
parts with respect to geomagnetic latitude, 
and the total conductivity in each region is 
assumed as follows: 

I. northern polar cap (colatitude 6=0°~20°) conductivity aK) (a~1), 

II. northern auroral zone (@=20°~25°) conductivity bKo (b>1), 

Ill. equatorial region (@=25°~155°) conductivity Ko, 

lV. southern auroral zone (@=155°~160°) conductivity 5K, 

V. southern polar cap (@=160°~180°) conductivity ako, 


where Ay is a constant value, factors @ and 0b being also constant. For the 


sake of mathematical simplicity, the following conditions are assumed, i.e. 


(i) The earth’s magnetic field is represented by a centred dipole of 
magnetic moment M, with its axis coincident with the rotation axis 
of the earth. The obliquity of the earth’s magnetic axis with the 
rotation axis is ignored. The vertical component of the earth’s 


magnetic field at the level of the conducting layer is given by 


H.=— ae cos 02=—2G cos @. e>) 
0 


(ii) The lateral motion of the air in the conducting layer is assumed to 
be derived from velocity potential , which can be expressed by the 


surface spherical harmonics on the layer, namely 


pans Wes 3 RPM (cos @) sin (MA+a§) , (2) 
NM NM 


316 N. FUKUSHIMA 


where PM (cos @) denotes Schmidt’s normalized surface spherical 
harmonic function, and 2» is longitude reckoned from the midnight 


meridian. 
As in the ordinary dynamo theory, we introduce current function, J. The 


differential equation for J in each region is 


oA arin Ble (tinge Obadpesss OY po lDsilgg } 
inf @z 60 (sino )= 2KG{ cot 0 a ae (sin @cos 8 0) (3) 


or 


sin 6 
2N+1 


ah BLN OF gd ia M 
sinbane * 00 (sine )=-2 3 2KGRN 


[{CN-+12—-1} {(N+1)2—M?}? Pyif (cos 8) 
+(N2—1) (N2- M2)? PX, (cos 6)] sin (MA+a), (4) 


where K represents the conductivity in respective regions, namely @Ko in the 
polar caps, bKy in the auroral zones and Kp in the equatorial region. Then, 
the solution of J can be obtained as follows: 


in region I, J=—2aKoGkNl > SM@(6) sin (MA+a¥) 
+5) CTUm@) sin (mr.+a%)], 
in region II, J=—2KGkNS pay S¥ (0) sin(MrA+a%) 
+2 {CF Um(O) sin (m+. a3) +C3 Vin (8) sin (nr+ a¥)}I, 


in region III, J=~2K6[S pay s¥(0)sin (Mata 


+5{CT Um(8) sin (rA+aT)+CF Vn(O) sin (mA+a% 3], fe 
in region IV, J=—2KClS > St (@) sin (Mr+a¥) 
+E{CEUm(O) sin (d+ a8)-+ CF VinC 8) sin (mr>+aF)}1, 
in region V, J=—2K0GCS = SHO) sin(MA+a¥) 
+3 CVn (@) sin (m4 at) ), 
where 
Ni(N+1y2—iey? 
Mc 9)\—pM M 
} 
(V+) WV—M?)" 
+ NQN+1) Py“ (cos >| ; (6) 


6 
Um(@)=tan™ om (for m=1), log tan = (m=0), 


Vin(@) = cot” sd (for m=1), 1(m=0), 
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and C?, a? (i=1,2, ..., 8; m=0,1,2, ...) are constants. The first term in 
the right hand side of the equations in (5) is particular solution of J, and 
remaining terms correspond to general solutions in each region. Groups of 
arbitrary constants C7? and a” must be determined by the condition that the 
normal component of electric current and tangential component of electric 
field given respectively by 


a tees 


Rosin 60r ” kK R00 


are continuous along the four boundaries of the five regions. From the boundary 


condition, it arises that a’ (7=1, 2, ..., 8) must be identified with a¥ for 
M=m. 

It will be convenient to consider that the expression for the current func- 
tion J consists of two parts, that is, 

J=J q +J D> (7) 
where the current-system derived from Jg corresponds to S,-field, while that 
from Jp causes Sp-field. In other word, Jqz is the solution in the case of 
uniform conductivity over the whole conducting layer. The form of J is 


5 Eee. papa SN) sin(MarA+aX), (8) 


everywhere through the five regions as in the ordinary dynamo theory, while 
Jp is the remaining parts in eq. (5). 
In the most simple case where assumes to be 
=F} (cos 8) sin (A+a}), 
the calculated current-system for Sp-field under the assumption a=1, b=19 


7 
and ai=—— is given in Fig. 2(b), which shows a similar mode of current 


2 
flow as Chapman’s Sp current-system with its phase lag by 6 hours. If we 


adopt ai=0 instead of rae the difference in the phase of both current- 


systems disappears. The S, current-system derived from Jz in the above 
case is given in Fig. 2(a) for comparison. The phase of the S, current- 
system in Fig. 2 (a) should be advanced by 24° in order to identify with P} 
component of the result of the actual analysis of geomagnetic data in 1905 by 
S. Chapman [1]. In both figures (a) and (b), electric current of 2.34x103 
amp. flows between successive stream lines in the direction indicated by 
arrows in the figures, so that the current intensity for the S, current-system 
agrees with Chapman’s analysis. By comparison of the current intensity of 
the calculated Sp current-system with that of Chapman’s one, it will be 
noticed that the calculated current-system is far less in its current intensity. 
In the above calculated current-system, total amount of the eastward or 
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Fig. 2. Calculated current-system from a simple dynamo theory for Sp-field (b), and 
the simultaneous Sg current-system (a). 


westward auroral zore current is 2.93x104 amp., while it amounts to 27.5x 10! 
amp. in Chapman’s idealized current-system. The total amount of the auroral 
zone current increases with increase in the value of 0, the ratio of the electrical 

conductivity of the auroral 
x o%amp zone to that in the other 
regions. The relation between 
the current intensity and the 
electrical conductivity in the 
auroral zone is calculated in 


Fig. 3, in which it is seen 
that the auroral zone current 
does not become infinitely 
large in its amount, but ap- 


Total amount of the auroral zone current 
> a 


2 4 6810 20 40 6080100 200 S00 /000 b proaches to an asymptotic 


. p ; ; value. This result is attribut- 
Fig. 3. The increase in the maximum total 


amount of the auroral zone current with the ratio ed to the following circum- 
(b) of the electrical conductivity of the auroral stance. In general, when the 
zone to that of adjacent regions. conductivity of the auroral 


zone becomes larger compared 
with that in the other regions, additional current intensifying S, current there 


flows along the auroral zone, and in consequence electric current of opposite 
longitudinal direction flows by the electrostatic field due to heterogeneous 
distribution of electric charge along the auroral zone, which is caused by the 
additional current flow along the auroral zone. Thus, an additional current 
system, namely Sp current:system is formed. Since the wind velocity derived 
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from velocity potential yl and the vertical component of the earth’s permanent 
magnetic field are larger near the pole, the induced electromotive force by the 
dynamo-action is always larger near the inner boundary of the auroral zone 
than near its outer boundary. Current density of the auroral zone current is, 
therefore, larger near the inner boundary of the auroral zone than near its 
outer boundary. When the electrical conductivity of the auroral zone becomes 
larger and larger, the electric current flow along the auroral zone becomes 
more intense. At that time, the electrostatic field, which produces the electric 
current flow of opposite direction in the auroral zone, also increases simulta- 
neously. Then, the electric current along the auroral zone cannot increase 
proportional to the conductivity of the auroral zone. When the conductivity 
of the auroral zone becomes more than a finite value, electric current along 
the outer boundary of the auroral zone is reversed, because the electromotive 
force due to the electrostatic field exceeds the induced electromotive force 
near the outer boundary of the auroral zone, making a circulatory current 
flow along the inner and outer~boundaries of the auroral zone itself. If we 
define the total amount of the auroral zone current as the total electric 
current along the auroral zone in the same direction with the induced 
electromotive force by the dynamo-action there, the total amount of the 
auroral zone current does not become infinitely Jarge even in the extremum 
case b>co. The asymptotic value of the total amount of the auroral zone 
current is 5.99x104 amp., which is only 22% of the total auroral zone current 
in Chapman’s idealized Sp current-system, so far the inner and outer bounda- 
ries of the auroral zone assume 20° and 25° in geomagnetic colatitude respec- 
_ tively. 

In the next place, the change in the total auroral zone current with the 
breadth of the auroral zone was calculated in the case of a=1 and b=10. 
The result is summarized in Table I, where the amounts of the polar closing 


Table I. Change in the total amount of the auroral zone current 
Cunit in 104 amp.) with the breadth of the auroral zone. 


Sis ite of the auroral | 992_959 99°28 = 20°—30°_-17°.5—-25°_| Chapman’s 

Total amount of the 2.93 3.42 3.58 3.46 27.50 
auroral zone current 

Amount of polar clos- 1.14 1.20 117 1.22 22.50 
ing current 

Amount of equatorial | 4.79 222 2.41 29g VR BbO 
closing current 


current (the current circulating through the auroral zone and the polar cap) 
and the equatorial closing current (the current circulating through the auroral 
zone and equatorial region) are also given, and each elements in Chapman’s 
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idealized current-system are written for comparison. The amount of polar 
current does not so much change with the position of the outer boundary of 
the auroral zone, so far as the inner boundary assumes to be 20°. According 
to the recent analysis of the auroral appearance frequency by E.H. Vestine 
[45], the zonal area in geomagnetic colatitude 17°.5~25° nearly coincides with 
the oval region, in which the auroral frequency is more than 80% according to 
his definition. Breadth of the auroral zone, 17°.5~25°, seems to be reasonable 
also from the analysis of the intense auroral zone current by T. Nagata [44] 
and others. In the case of the above breadth of the auroral zone, the 
magnitudes of both polar and equatorial closing currents do not show any 
remarkable difference with those in the preceding cases, as will be noticed in 
Table I, showing little dependency of the mode of current flow on the position 
of the inner and outer boundaries of the auroral zone. 

In the above calculation, the air motion in the upper conducting shell was 
assumed to be derived from the most simple form of the velocity potential 
wi, and higher order terms were neglected. Here the influence of higher 
order terms was examined. In the analysis of S,-field by S. Chapman with 
the data in 1905 [1], the magnetic potential of its external origin part, We, 
can be approximately expressed as 


W.(0, ’)=R{ 7.7" P} (cos 8) cos (A +24°) +4.47P2 (cos 8) cos (2\’+202°) 
+2.27P? (cos 0) cos (3/+35°) +0.77 P! (cos 8) cos (4A/ +232°)}, (9) 
where 2/ denotes longitude measured eastward from the noon meridian, and 


R the radius of the earth. Then, the velocity potential of the air motion in 
the conducting shell of radius Ro can be derived as 


R 5 /R 
Sai KG {7.7 -( are (cos 8) sin (A+114°) 


V3 
1 R 
Lae Biv _( i) P3 (cos 8) sin (2X-+292°) 


; 
+2.2 sls eal P3 (cos 6) sin (3X+125° ) 


55 / Ro 
+0.7 Faas. Pi (cos @) sin can+-322°)}, (10) 
by the simple dynamo theory for the case of uniform conductivity over the 
whole conducting shell, where 2» denotes longitude measured eastward from 
the midnight meridian. The result of the calculation by means of the form of 
the velocity potential in (10) is illustrated in Fig. 4, where electric current of 
2.34103 amp. flows between successive stream lines. Since the maximum 
current intensity of Jp along the auroral zone derived from Ww and 3 is 40% 
and 3% of that from >} respectively, and the current intensity derived from 


A aye Votre ie c ‘ . 
Wi is much less, 4 is ignored in the numerical computation for the current- 
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system in Fig. 4. This current-system shows better agreement in the distribu- 


tion of stream lines with the current-system of Fig. 11 in Chapter II, if the 
phase of the current system be advanced about 120°. 


Fig. 4. Calculated Sp current-system when, the higher 
order terms in the velocity potential of the air motion is 
taken into account. 


§2. Double layer hypothesis in dynamo theory 


It was shown in the preceding paragraph that the calculated current- 
system for Sp-field is far weak in its current intensity and has considerable 
phase difference compared with the observed ones. Fig. 5 shows the direction 
of the polar parallel current obtained by many researchers. It directs to the 
meridian of 12h in local time in Chapman’s idealized current system, to that 


1. Calculated current-system from Jp (ot =90°). 


0m Calculated current-system from Jp (at =114°). 
6 ig 1 3. After Chapman’s idealized current-system for 
6 Sp-field. 
2 4. After Hasegawa. 
5 12 5. After Hatakeyama, Nagata, Silsbee and Vestine, 
and Harang. 
4 3 6. Average current-system of four bays shown in 
Fig. 11. 


Fig. 5. Direction of the polar parallel current. 
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of about 10h by M. Hasegawa [10], to about 7.5h meridian by T. Nagata [27], 
E. H. Vestine [28] and L. Harang [39]. In the current-system of geomagnetic 
bays, the polar parallel current directs to the meridian of about 7.5h by H. 
Hatakeyama [33] or H.C. Silsbee and E.H. Vestine [35], while to about 5h 
meridian in the case of the average current-system of four examples of bays 
by the writer shown in Fig. 11 in §3. On the other hand, the polar parallel 
‘current, which should be theoretically derived by the air motion responsible 
also for S,-field, directs to 16.4h meridian, when the velocity potential of the 


air motion assumes to be in the form of W}, and the phase angle at} is 114°. 


qr. e 
In the calculation shown before, al is taken to be > for simplicity, resulting 


that the direction of the polar parallel current directs to 18h meridian. It will 
be said that the direction of the polar parallel current of the phase of the 
calculated current-system is rather inverse with the observed ones. 

In order to remove the above difficulties concerning. the current intensity 
and the phase difference of the calculated current-system, the writer proposed 
a double layer hypothesis of dynamo theory [42]. In the hypothesis, the 
presence of two conducting layers concentric with the earth is assumed, total 
conductivity. of these two layers being /AK) and (/+1) > respectively, as 
illustrated by Fig. 6. The air motion assumes to be inverse to each other, 


The first layer. of total 
conductivity 1Ko 


The second layer of total 
conductivity (J+1) Ko 


Fig. 6. Assumed two conducting layers concentric with the earth. 


and the air motion in the second layer of conductivity (J+1) Ky be quite the 
Same as in the ordinary dynamo theory for Sg-field. Then, the observed 
magnetic variation on the earth’s surface corresponds to the superposition of 
the two current-systems generated in the two conducting layers. It is assumed 
further that, on magnetically disturbed days, the first layer of conductivity 
1K alone is disturbed, and its conductivity becomes very large in the auroral 
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zones. Then, the resultant S, current-system is quite the same one as that 
produced in one layer of conductivity Ky, while Sp current-system appears 
l-times intensified in its current intensity and reversed in the direction of 
electric current flow compared with that calculated before. From Table I of 
the preceding results, appropriate value of / will be estimated to be about 5. 

The integrated electrical conductivity of the upper atmosphere has been 
estimated to be about 5x10-8 e.m.u. in the recent studies by T. Nagata [8] 
and K. Maeda (47 and others. According to the above double layer hypothesis 
of dynamo theory, this value corresponds to (2/+1) Ky. If J assumes 5, Ky 
should b2 about 5x10-9 e.m.u. The magnitude of wind velocity in the upper 
conducting layer, which is derived from the ordinary dynamo theory for S,- 
field with the value of Ao as 3x10-§ e.m.u,, is several tens cm/sec. at most. 
When 5x10-9 e.m.u. is adopted as Ky, wind velocity amounts to the order of 
a hundred m/sec. to produce the same magnitude of g20magnetic variation. 
The above wind velocity of the air motion in the upper atmosphere seems to 
coincide well with recent various results of study on the movement in the 
ionospheric region. 

Double layer hypothesis in dynamo theory seems to interpret the causation 
of Sp current-system as well as Sg current-system, as mentioned before. 
However, reliable proof is not yet obtained in the present stage for the 
circulatory motion in the upp2r atmosphere, which is favourable for the present 
assumption in the air motion of the upper atmosphere in the double layer 
hypothesis. 


§3. Calculation with more plausible assumptions on the electrical conduc- 
tivity in the upper atmosphere 


The current-system for the disturbance field in high latitudes at the time 
of polar magnetic storms or geomagnetic bays from the actual analysis of 
observed data shows that the current intensity in the auroral zone is larger 
in the dark hemisphere compared with that in the sunlit hemisphere in 
general, as already shown in Fig. 1 in Chapter II. Such a characteristic 
might be due to an inhomogeneous distribution of the electrical conductivity 
along the auroral zone. Dynamo theoretical calculation with more plausible 
assumptions on the distribution of electrical conductivity was treated by the 
writer by means of a perturbation method [43]. In the followings, the calcu- 
lated result will be briefly given. : 

The assumption on the distribution of total conductivity in the calcula- 
tion is 

K=aky in the polar caps, 
K=bKy) (1+ycosX), in the auroral zones, Cd) 
K= K,(1—y'sin@cos’) in the equatorial region, 
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where a~l1, b>1, y and »7/<1. The conductivity in the auroral zones at 


midnight is = times that at noon. This assumption seems to be adequate 


from the current distribution along the auroral zone in the disturbed state, or 
from M. Hasegawa’s suggestion [48]. In the equatorial region, the conductivity 
is chosen to be proportional to 1++’cos X, where X denotes zenith distance of 
the sun. If the sun’s declination assumes zero, the form is written as above. 
Other assumptions are quite the same as in the preceding calculation. 
The current function J can be considered to be composed of 
J=JatJo; (2) 

as in the same way in the preceding case. Further, if Jgo and Jpo denote the 
expressions for Jz and Jp when y=1'=0, Jq and Jp will be considered to have 
the composition as follows: 

Ja=AJatyS' a 

Jp=J vot] pv + YI py 
where y'J’qy, and y/J’py are perturbation terms corresponding to /KqcosX, 
and yJ'py is also a perturbation term, which corresponds to rybKocosX. 
When wW assumes to have a form 


W=kIPi(cos 6) sin (A+a}), 
the expressions for respective components are as follows: 


P\(cos 6 
Jon=—2IGGk aheeel sin (A+a}) 2 


(3) 


(4) 
J gv=—2K GR aH { (cos 30—cos@) sin (24-++a}) + (cos 38—3cos 8) sin aj} , 
over the whole conducting layer, and 
a—1 Px cos 6) 


0 
—2aKyGki { a 2/3 +C, tan a sin (A+a}) 


in the northern polar cap, 


b—1 Pkcos 0) é 6). 
b 2/3 +C, tan 9 + C3 cot y sin (A+al) 


in the northern auroral zone, 


—2bK iG kh { 


) 
—2K)G ki (C tan > + Cs cot = sin (A+a}) 


Jpo= 2 (5) 


in the equatorial region, 
b—1 P(cos 8) 
b 2/3 


--2bKyG k} { +C, tan g +C7cot a sin (A+a}) 


in the southern auroral zone, 
1f@—1 Pk(cos 0) A oe 
—2aK)Gk; eS 2/3 — +Cgcot | sin (A+a}) 
in the southern polar cap, 


3 
Jp y=‘ 


/ 
Jpy= 
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6 
—2aKy Gk} {p, tan? 2 sin (2X+a})+B, sin af 


in the northern polar cap, 


—2bKo GH| {cot 0+ eats =) cosec 6+ < sin 20+ D, tan? g 


+Dz3 cot? 5} sin (2A+al)+ te sin 20+ B, log tan e 


+B;} sin at | in the northern auroral zone, 


—2K Gk} | {Patan £ +D;cot? ‘a sin (2X+al)+ {2 log tan 
(6) 
+Bs) sin at | in the equatorial region, 


—2bK Gri [ {cot 0+ as —@) cosec 0+ < sin 20+ Dg tan? 6 


+ D, cot? sa sin(2\+at)+ ee sin 260+ B; log tan + 


+Be sin at | in the southern auroral zone, 


—2aK Gk {Ds cot? ¢ sin (2+ a})4+Bs sin al} 


in the southern polar cap, 


6 1 
—2aKyGki | {> tan? ae oye 


(cos 36—cos @) \ sin(2A+a}) 


4 Be a3 (cos 30—3 cos ®} sin at | 


in the northern polar cap, 


—2bKy) Gk} Ki D,! tan? 4 + D,! cot? 2 = a5 (cos 38—cos a} sin (2X 


2 246 


+alh)+ {By log tan g. + B3!— aie (cos 30—3 cos 8) \sin al | 
in the northern auroral zone, 


—2K Gk} | 4 D4’ tan? o +D;! cot? “| sin(2Xr bad) +{ By log tan “ 
ai te 2 2 2 (7) 


+ B,/ } sin at| in the equatorial region, 


6 1 4 
—2bKyG ki I Dé tan? £ +D/ cot? Do Tas (cos 30—cos 0) \sincen 


ar aS Sr {By tog tan. +B/— ae (cos 39—3 cos 4) } sin al | 
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in the southern auroral zone, 


—2aK,G k} K Ds! cot? & a a (cos 380—cos 0) sin(2r-+a}) 


+{ By— ata (cos 38—3 cos@) sin at | 


in the southern polar cap, 


where D;, D,’, B;, B;! (i=1, 2,......, 8) are constants, which should be deter- 
mined from the boundary condition. In the practical numerical computation, 
the values of a and b were chosen to be 1 and 10 respectively, and aj was 


assumed to be ae according to the result in the preceding paragraph. 


When 1>¥7>0, the current intensity in the auroral zone is more intense 
in the dark hemisphere and less in the sunlit hemisphere, while the equatorial 
current intensity shows only a slight decrease in the dark hemisphere and 
increase in the sunlit hemisphere with the value of y. The ratios of the total 
amounts of the current flowing in the auroral zone and the equatorial region 
of the dark hemisphere to those of the sunlit hemisphere, denoted by f(A) 
and f(£) respectively, are given in Figs. 7 and 8. At that time, the current 


OP 502 e038 HOA OS 506 


y 
Fig. 7. Change in the ratio of the Fig. 8. Change in the ratio of the 
intensity of the auroral zone current in current intensity of the equatorial region 
the dark hemisphere to that in the sunlit in the dark hemisphere to that in the 
hemisphere with y. sunlit side with y. 


intensity in the polar caps does not show remarkable change, its change being 
less than 10%. In Fig. 9, the calculated current-system when y= and y/=0 
is shown. The current-system for geomagnetic bays obtained by H.C. Silsbee 
and F.H. Vestine [35] is also shown here for comparison. The mode of current 
flow in Fig. 9 shows good similarity with that in Fig. 10, although there is 
a discrepancy in phase of about 30°. On the other hand, the mode of current 
flow in the equatorial region is considerably deformed with increase in the 
value of y’, while the current intensity in both the auroral zones and the polar 


Polar Magnetic Storms and Geomagnetic Bays 327 


j2h 
Fig. 9. The calculated current-system Fig. 10. The current-system for the 
for Sp-field and the bay disturbance when bay disturbance obtained by H.C. Silsbee 
y= and y/=0. and E. H. Vestine. 


caps is scarcely affected. Therefore, the asymmetry of the conductivity in the 
equatorial region will be not so significant in the present problem, so far as 
the current distribution in the higher latitudes is concerned. 

In Fig. 11, the average current-system for the maximum stage of four 
bays, which will be dealt with later in Chapter IV, is illustrated. This current- 
system was compared with many calculated ones derived from Jp, which were 
obtained by varying the magnitudes of a@ and 8 and y. It was concluded then 
that the distribution of the electrical conductivity in the upper atmosphere 
corresponding to the current-system in Fig. 11 can be best expressed by 


@=1.2,,6=9-and y=0.6 


in the present model. In other words, the conductivity of the auroral zone is 
about 14 times that of the equatorial region near midnight and about 4 times 
near noon, while the conductivity in the polar cap is comparable or somewhat 
large compared with the equatorial one, at the maximum stage of geomagnetic 
bays. In Fig. 12, the distribution of stream lines of the calculated current- 
system under the assumption a=1.2, b=9 and y=0.6 is given for comparison. 
Since the current intensity in the calculation is determined by the product of 
ki and Ko, the value of ki Ky may be appropriately chosen, so that the numer- 
ical difference of current function between the pole and the equator in the 
both current-systems in Figs. 11 and 12 agrees with each other. The distribu- 
tion of the values of the current function with respect to geomagnetic latitude 
in the noon and midnight meridians for both observed and calculated current- 
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WORLD, NORTHERN HEMISPHERE IN GEOMAGNETIC COORDINATE WORLD, NORTHERN HEMISPHERE IN GEOMAGNETIC COORDINATE 
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Fig. 11. Average current-system of Fig. 12. Current-system calculated from 
four bay disturbances at their maximum dynamo theory under the assumption a=1.2, 


stage. Electric current of 1.0105 amp. b=9 and y=0.6. 
flows between successive stream lines in 
the direction indicated by arrows. 


systems is compared in Fig. 13, in which we see good agreement of the dis- 
tribution of current function with each other. The width of the auroral zone, 
which is taken to be 5° in latitude in the calculation, may be somewhat wider 
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Fig. 13. Distribution of the value of current function with latitude 
in the noon and midnight meridians. The top and the bottom of the 
figure correspond to the average current-systems given in Fig. 11 and 
to the calculated one in Fig. 12, respectively. 
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in the observed one. However, no serious alteration will be expected in the 
above result for the distribution of the electrical conductivity in the upper 
atmosphere, when the auroral zone is assumed to be wider than 5° in 
latitude. Thus, the asymmetry of the intensity of the auroral zone current in 
the geomagenetic disturbance field will be attributed to an inhomogeneous 
distribution of electrical conductivity along the auroral zone from the stand- 
point of dynamo theory. 


$4. A model experiment for the distribution of stream lines in Sp current- 

system 

In the numerical computation of Sp current-system, only the cases of 
rather idealized distribution of electrical conductivity in the upper atmosphere 
are able to be obtained because of mathematical difficulty, and the calculation 
under an arbitrary distribution of electrical conductivity or electromotive force 
in the conducting shell can be hardly dealt with in practice. In order to find 
the distribution of electric current flow in such a case, a model experiment 
was carried out with technical assistance by T. Oguti, and the results were 
compared with calculated ones. 

The electromotive force (FE) generated in the electrically conductive iono- 
sphere produced by the dynamo-action is written as 


rot E= rot(v-H] (1) 


were v is lateral velocity of the air mass and Hf is the earth’s. permanent 
magnetic field at the level of the conducting layer. E consists of induced 
electromotive force (F;) produced by the dynamo-action and electrostatic field 
(E;) due to heterogeneous distribution of electric charge at the time of current 
flow, as 
E= Ej;+E;.- (2) 
The resulting electric current (i) in its stationary state can be obtained from 
the equation 
ix=oE, (3) 
where o« denotes electrical conductivity of the ionosphere, and is a function of 
position over the conducting shell. At the time of polar magnetic disturbances, 
ao can be considered to be larger in the auroral zone compared with the 
remaining parts, polar cap and equatorial region. As already mentioned, the 
current function (J) may be separated as 
J =JatJIo; (4) 
where J, means the current function obtained in the case of uniform con- 
ductivity over the whole conducting layer, while Jp is the additional current 
function when the conductivity of the auroral zone becomes higher than that 
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of the adjacent regions. The induced electromotive force produced in the polar 
caps and the equatorial region contributes only for the formation of Sg, current- 
system, because no additional induced electromotive force appears there when 
the conductivity in the auroral zone alone becomes high. On the other hand, 
the increase in the electrical conductivity in the auroral zone results in the in- 
crease in the current intensity along the zone, and consequently additional 
electrostatic field arises there. Then, leakage current flows in the polar cap 
and the equatorial region by this additional electrostatic field, forming a closed 
current-system with auroral zone current, which is Sp current-system itself. 
Therefore, we can obtain a similar current distribution as Sp current-system 
in the polar cap and the equatorial region by a model experiment, in which 
an appropriate distribution of electrostatic field is given along the auroral 
zone, so far as we do not concern the current distribution along the auroral 


zone. ; 
A rubber ball of radius 7.7cm was applied for the model of the conductiv 


ionosphere, its surface being painted with carbon powder mixed in glue. The 
model of the ionosphere and a part of apparatus are seen in the photograph. 
Thin thread was spread over the spherical surface to indicate locality. Electro- 
static potential was supplied by dry cells by means of dividing resistances at 
16 points along both the northern and southern auroral zones, as illustrated 
in Fig. 14 and the photograph. Arbit- 
rary choice of electric potential along 
the auroral zone can be allowed by 
changing the contact points of the divid- 
ing resistances. Nichrome line of resis- 
tance 20 was adopted as the dividing 
resistance. Equipotential lines over the 
whole ‘spherical surface were obtained 
with the aid of a potentiometric method. 
Then, stream lines of electric current 
flow could be derived graphically so 
that equipotential lines and stream lines 
cross perpendicularly with each other. 
At first, the case in which the 
electrical conductivity in the polar caps 
and the equatorial region is a constant 
value throughout these regions, was 
examined. Mixture of carbon powder 
and glue was painted about 3mm. thick 


The model of the electrically conductive pon the rubber ball as uniformly as 
ionophere and a part of appiratus for the s 


experiment. 


could. Electric resistance between two 
points at lcm apart on the sphere 
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Fig. 14. Supply of electric potential along the auroral zone in the 


model experiment, and a potentiometric method of measuring the 
equipotential line over the spherical surface. 


showed about 1522. Uniformity of the electrical conductivity was tested by 
comparing distributions of stream lines by the experiment with those already 
known by theoretical calculation. In practice, the distributions of stream lines 
when two electric doublets of the same moment were put symmetrically on the 
northern auroral zone, were compared with each other. Polar distance of the 
auroral zone was 22.5° in both experiment and calculation. Theoretical. cal- 
culation of stream lines in such a case will be mentioned in Chapter V, §7. 
As will be seen in Fig. 15, only a slight difference in the mode of current 


by experiment calculated 


Fig. 15. Comparison of experimental and theoretical distributions of stream 
lines for the test of uniformity of the electric conductivity over the spherical surface. 
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flow appears between them near the equator and the auroral zone. The difference 
in the distribution of stream lines near the auroral zone will be attributed to 
some increase in the electrical conductivity along the auroral zone, which 
originates in swell of carbon powder at the electric rodes along the auroral 
zone. The mode of current flow near the equator is not seriously dealt with 
in the present problem, so that the electrical conductivity may be considered 
to be satisfactorily uniform for the present purpose. 

It can be derived from the result in §2 that the current function Jp for 
the simplest case is expressed by 


Jp=—A tan = sin (A+a}) in the polar cap, (5) 
and — 
Jp=B cot es sin (A+aj) in the equatorial region, (6) 


when the northern auroral zone alone becomes highly conductive, and 
J p=2B cot @sin (A+ aj) = B(cot & — tan i sin (A+a}) (7) 


when both the northern and southern auroral zones become highly conductive 
symmetrically with respect to the equator. In the above equations, A and B 
are constants, which are determined relating to the position of the auroral 
zone. Let the electrostatic potential, from which the electrostatic field produc- 
ing. Jp current system is derived, denote by S, the following relations hold 
between Jp and S, namely 


Ofp os 


60 me orth O0Or ’ : 
We. yp ask ee 
sin @0X pe 
Then, S can be solved as follows: 
i 
S=A Ko tan S cos (A+a}) in the polar cap, (9) 


1 6 
S=B B&B cot cos (A+a}) in the equatorial region, (10) 


when the northern auroral zone alone becomes highly conductive, corresponding 
to (5) and (6) respectively, and 


ape ie - 
S=2B ie cosec 0 cos (X Fan =Bae- (cot air tan a cos(A+al) QD 


corresponding to the current distribution in the equatorial region given by (7), 
when both auroral zones become highly conductive. 
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On the other hand, the electric current flow over the conducting shell 
when the electric charge p of the form of 


p=po cos (A+) (12) 


is supplied along the auroral zone can be obtained as follows. First, the case 
when the electric charge p in its linear density distributes along a circle of 
radius 7, as shown in Fig. 16, is considered as a two dimensional problem. 
The charge distribution along the circle assumes to have the form of eq. (12) 
in (7, \) coordinate. Then, the electric potential S’ over the plane is given by 


‘ 2° : 
S!=2rpo “2 cos (A+a)) for r>%, 
(13) 


S!=2mrpor cos (A+a}) 5 for r<7. 


In the next place, S’ on the plane is projected on a spherical surface of radius 
Ro, which contacts with its top pole at the origin of (7, ») coordinate to the 
plane, as will be seen in Fig. 17. The straight line connecting the bottom pole 


Fig. 16. A cosine type distribution Fig. 17. Projection of a point on the plane 
of electric charge along the circle. to a point on the spherical surface. (¢ in the 
figure must be read as 2) 


of the sphere and a point P on (%, ) plane intersects the surface of the sphere 
at a point Q. A point P(%, ») is projected to a point @ (Ro, 6, %) in the 
spherical surface. The relation between Ro, 6 and r is given by 


vr=2R) tan ws (14) 


Do? 


Then, S’ on the surface of the sphere after the projection can be written as 


S/=4mRopy tan a cos (A+a}) for 0< 0, 
(15) 


S!=4rRopo tan? Bo: cot va cos (A+a}) for 0>0), 
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where 


%=2Ry tan a : (16) 


When the electric charge is given along two auroral zones symmetrically with 
respect to the equator, the electric potential in the equatorial region becomes 


S!=47rR. po tan? 2 cot o+ cot2 mb tan 5 cos (A+al) 
=A Ro po tan? -£6-(cot é + tan 2) cos (A+a‘l) 


=47 Ro py tan? cosec 6 cos (A+a}) for 7m—OQ)>0>0). (17) 


Comparing (15) and (17) with (9), (10) and (11), we see at once that the 
mode of current distribution on the conductive surface of the sphere given by 
(5), (6) and (7) can be obtained experimentally by supplying the electric 
charge of the form of (12) along the auroral zone. y?S=0 or p?S’'=0 is of 
course Satisfied in the above expression for S or S’. 

An experimental result under a cosine type distribution of the electric 
potential along the auroral zone of the northern hemisphere was compared 
with the calculated one in Fig. 18, while the distribution of electric potential 
supplied along the auroral zone in the model experiment is illustrated by Fig. 
19, in which the potential difference between two successive electric rodes is 


(a) by experiment (b) calculated 


Fig. 18. Comparison of experimental and calculated distribution of stream lines under a 
cosine type distribution of electric potential along the auroral zone. 
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chosen to be proportional to the cosine of the mean longitude of contacting 
rodes. It can be seen that the distribution of stream lines both from experi- 
ment and calculation shows nearly perfect agreement with each other. When 
the quite similar distribution of the electric potential is supplied also along 
the southern auroral zone, the comparison of the experimental and calculated 
results shows nearly perfect agreement with each other, too, the illustration of 
the result is omitted here. From these results, we find satisfactory agreement 
of the experimental current-systems with the calculated ones by the dynamo 
theory concerning the distribution of stream lines in the polar cap and the 
equatorial region. 


Fig. 19. Distribution of electric poten- Fig. 20. Distribution of electric poten- 
tial supplied along the auroral zone, tial supplied along the auroral zone, 
which results in the distribution of stream which results in the distribution of 
lines given in Fig. 18 (a). stream lines given in Fig. 21 (a). 


The experiment was extended for some cases under asymmetric distribution 
of the electromotive force along the circumference of the auroral zone. In 
those cases, the electric potential was supplied along both the northern and 
southern auroral zones symmetrically with respect to the equator. An example 
of the distribution of the supplied electric potential is illustrated in Fig. 20, in 
which the potential shows a consine-type distribution along each half circle, 
and the magnitude of the supplied potential along one side of half circumfer- 
ence is 2.7 times the other side. The distribution of stream lines obtained 
under the potential distribution of Fig. 20 is illustrated in Fig. 21 (a), while 
the calculated current-system when the conductivity along the auroral zone 
varies as 1+4cos(A+A)) with longitude, resulting the ratio of the maximum 
conductivity to the minimum one is 3 in this case, is given in Fig. 21 (b) for 
comparison, showing fairly good agreement with each other. In Fig. 22, the 
distribution of stream lines in the case when the magnitude of the supplied 
potential along the auroral zone of a hemisphere is 5 times as much as that 
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(a) by experiment (b) calculted 


Fig. 21. Comparison of experimental and calculated distributions of stream lines, where ~ 
the magnitude of electric potential along a half circle is 2.7 times as much as that in the 
other half circle in the experimental ‘case, while the maximum conductivity of the auroral 
zone is 3 times as Much as the minimum one in the calculated one. 


Fig. 22. Distribution of stream lines under Fig. 23. Distribution of electric po- 
the supply of electric potential along the auroral tential along the auroral zone, which 
zone as given in Fig. 23. results in the current-system of Fig. 22. 


in the other hemispere, while the supplied electric potential along the aurorai 
zone in the experiment is shown by Fig. 23. As will be noticed in «these 
figures, the distribution of stream lines in the current-system changes with the 
supplied electric potential along the auroral zone. The experimental current- 
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Fig. 24. Distribution of stream lines under Fig. 25. Distribution of electric po- 
the supply of electric potential along the auroral tential along the auroral zone, which 
zone as given in Fig. 23. results in the current-system of Fig. 22. 


system in Fig. 24 is obtained under the supply of the electric potential shown in 
Fig. 25, in which the potential difference between two successive points shows 
a cosine-type distribution along half a circumference, while a potential difference 
is given along only a part of the remaining part of the auroral zone. Many 
examples of current systems, which have similarity with above experimental 
ones, will be encountered in the following chapters. y 

Thus, the current-system under the complicated distribution of electrical 
conductivity along the auroral zone can be obtained by a simple model experi- 
ment by supplying the appropriate electric potential along the auroral zone. 
In other words, distribution of electrical conductivity in the upper atmosphere 
along the auroral zone at the time of magnetic disturbance of polar origin will 
be roughly estimated, after the comparison of the experimental result »under 
various conditions and the actual current-system for the disturbance field. 
The experiment was further extended also for various cases of non uniform 
conductivity over the surface of the sphere. However, no serious alteration of 
the mode of current flow was noticed so far as the maximum regional con- 
ductivity remains a few times as much as the minimum one over the conduct- 
ing shell in the model experiment. 
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Chapter IV. Development and Decay Processes of 
Individual Geomagnetic Bays 


§1. Average disturbance field of geomagnetic bays and individual bays 


The equivalent overhead current-system corresponding to the average 
disturbance field of a number of bays was obtained by H. Hatakeyama [33], 
H.C. Silsbee and E.H. Vestine [35] and others as already mentioned before, 
and geomagnetic bays have been considered to be a temporary appearance 
of Sp-field within a few hours. It is already shown in Chapter III, that the 
current-system for bays can be derived from the dynamo theory, when the 
production of highly conductive area takes place along the auroral zone in the 
upper atmosphere. However, the average characteristics of the disturbing force 
of bays observed in the middle latitudes, which were already dealt with in 
Chapter II, suggest that the current-system for bays reveals some deforma- 
tion during its development and decay processes. Since geomagnetic bays are 
considered to be a disturbance of polar origin, the above-mentioned deformation 
of the current-system is supposed to be attributed to some change in the 
physical condition, say electrical conductivity, of the upper atmosphere in the 
high latitudes with the progress of time. The disturbing force of geo- 
magnetic bays in the high latitudes shows sometimes fairly complicated 
variation, even when that in the lower latitudes reveals itself as a rather 
typical bay-type trace on the magnetogram. ‘Therefore, individual bays will 
be not always interpreted by the development and decay of such a current- 
system as an idealized Sp current-system. Detailed study on the disturbing 
force along and near the auroral zone is required for investigating the physical 
mechanism for causation of bays, with respect to the statistical distribution 
and also to individual examples. 

The writer examined four examples of geomagnetic bays, which took place 
during the Second International Polar Year, with the aid of a number of 
magnetogram copies, among which sufficient number of polar observations were 
included. The analysis was especially laid emphasis on the world-wide change 
in the distribution of the disturbing force of bays with the progress of time. 
Some suggestive conclusions on the mechanism of geomagnetic bays were 
obtained by the study [46]. The results will be briefly summarized in the 
followings. 


§2. Individual examples of development and decay processes of geomagnetic 
bays 
Geomagnetic bays dealt with in the present study were those occurred at 


April 3, 1933, 10—13h G.M.T., 
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Table I. Abbreviation and position of observatories. 


Observatory Abbreviation ® A g ? 

Thule TH 88.0 0.0 | 765 | 2011 
Godhayn GO 79.8 32.5 69.2 305.5 
Scoresby Sund Ss 75.8 81.8 70.5 338.0 
Angmagssalik AN - 743 52.7 | 65.6 37.6 
Sveagruvan SV 73.9 130.7 779 16.8 
Calm Bay CB 71.5 153.3 80.3 52.8 
(Bjorndja) BJ 71.0 | 1247 | 74.5 19.2 
Julianehaab JU 70.8 35.6 60.7 314.0 
Fort Rae FR 69.0 290.0 62.8 243.9 
Tromso TR 67.1 116.7 69.7 18.9 
(Petsamo) PE 64.9 1258 69.5 31.2 
(Matotchkin Shar) MS 64.8 146.5 73.3 56.4 
College co €4.5 259.4 64.9 212.2 
(Sodankyla) so 63.8 120.0 | 67.4 26.6 
Dickson DI 63.0 161.5 73.5 80.4 
(Kandalaktcha) KN 62.5 124.2 67.1 32.4 
(Lerwick) LE 62.5 88.6 | 60.1 358.8 
(Dombéas) DO 62.3 | 100.0 | 621 9.1 
Meanook ME 61.8 301.0 54.6 246.7 
Sitka SI 600 275.4 57.0 224.7 
(Eskdalemuir) ES 58.5 82.9 55.3 356.8 
Lovo LO 58.1 105.8 59.4 17.8 
Agincourt AG 55.0 347.0 43.8 280.7 
Abinger AB 54.0 83.3 51.2 359.6 
(Gross Raum) GR 53.3 105.3 | 54.8 20.5 
Yakutsk YA 51.0 193.8 62.0 129.7 
Cheltenham CH 50.1 350.5 38.7 283.2 
San Miguel SM 45.6 50.9 37.8 334.4 
Ebro EB 43.9 79.7 40.8 0.5 
(San Fernando) SF 41.0 713 | 36.5 353.8 
Tucson TU 40.4 312.2 32.2 249.2 
(Toyohara) TY 36.9 203.5 47.0 142.8 
San Juan SJ 29.9 Sie 18.4 293.9 
(Teoloyucan) TE 29.6 327.0 | 19.8 260.8 
(Helwan) HE 272 | 106.4 | 29.9 31.3 
Kakioka KA 26.0 206.0 36.2 140.2 
Aso AS 22.0 198.0 32.9 131.0 
Honolulu HO 21.1 266.5 21.3 201.9 
Antipolo AT 3.3 189.8 | 146 121.9 
Huancayo HU — 0.6 353.8 |—12.0 284.7 

MO — 2.7 1143 2.0 45.4 


Mogadiscio 
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(Continued) 


Observatory Abbreviation D A ? vn 
Flisabethville EL (12% | 940 ou | 275 
Apia AP —16.0 260.2 |—13.8 188.2 
(Tananarive) TN —23.7 112.4 |—18.9 47.5 
Cape Town CT —32.7 79.9 |—33.9 18.5 
Watheroo WA —41.8 185.6 |—30.3 115.9 
Toolangi TO —46.7 2208 |—37.5 145.5 
(Orcades de Sud) Os —49.9 182 |—60.7 315.2 


Notes: ® Geomagnetic latitude (The position of north geomagnetic pole 
assumes ¢=78°.5N, 2=291°.0E) 
A Geomagnetic longitude 
¢ Geographic latitude 
i Geographic longitude 


April 9, 1933, 9—12h G.M.T., 
April 10, 1933, 14—17h G.M. T., 
April 23, 1933, 5— 8h G.M. T. 


Copies of magnetogram at 33 observatories were collected, 26 of which being 
in the northern hemisphere. Besides those records, tables of hourly values 
could be utilized for reference for 15 observatories. The abbreviation and 
the position of the stations are summarized in Table I, in which the latter 15 
stations are written in parentheses. The disturbing force of bays was defined 
to be the deviation from the straight line connecting values at the beginning and 
ending of the above-mentioned intervals on the magnetogram for each element. 

The variation in the disturbing force in the high latitudes with time 
is fairly complicated, while it seems that the disturbing forces at various 
localities in the middle and low latitudes increase and decrease almost 
simultaneously. However, there seem to be some instances at which the 
records at considerable number of stations show optimum values of the 
disturbing forces. These selected times are 


on April 3, 11h06m, 11h30m, 11h52m, 12h40m, 
on April 9, 09h40m, 09h50m, 10h09m, 11h00m, 
on April 10, 14h30m, 15h00m, 15h19m, 15h30m, 
on April 23, 05h48m, 06h00m, 06h10m, 06h20m, 06h40m, 07h20m, 


where the above underlined times are those at which the disturbing force in 
the middle and low latitudes reaches its maximum magnitude. Average cur- 
rent-system of four bays at above respective times was already. shown in 
Fig. 11 in Chapter III, and discussed from the standpoint of dynamo theory. 
Progressive change in the current-system of individual bays was looked over 
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Fig. 1. Distribution of the overhead current 
arrows and the equivalent current-system for 
the bay disturbance on April 3, 1933. 
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Fig. 2. Distribution of the overhead current 
arrows and the equivalent current-system for 
the bay disturbance on April 9, 1933, 
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Fig. 3. Distribution of the overhead Fig. 4. Distribution of the overhead 


current arrows and the equivalent cur- 
rent-system for the bay disturbance on 
April 10, 1933. 


current arrows and the equivalent cur- 
rent-system for the bay disturbance on 
April 23, 1933. 
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by comparing current-systems at some stages during development and decay 
of bays. Current-systems were roughly estimated from the distribution of 
current arrows over the world for three representative times for respective 
bays. The estimated overhead current-systems are shown in Figs. 1-4. One 
can see at once that the current-system varies its aspect during its develop- 
ment and decay processes. General characteristics of the progressive change 
in the current-system are summarized in the following paragraph. 


§3. Some remarks on the change of the disturbance field of bays over the 
world with time 


In examining the change of the distribution of the disturbing force of bays 
shown in Figs. 1-4, following characteristics are noticed. 
(i) Development and decay proces: 
ses of bays cannot be considered to be 
the simple intensification and diminution 
of such a current-system of the average 


bay, because the disturbing force of bays Sha 

does not always increase or decrease G Ae peat 

simultaneously over the whole world, 

especially in the high latitudes. @ Apr. 9 
(ii) In general, the disturbing force A i LS 

of bays at the middle and low latitude Ne 

stations takes its maximum value nearly ert pa Ble Se Ls 

at the same time when the westward 9 = ae a [2h 


auroral zone current in the dark hemi- 
sphere reaches its maximum, and the 
eastward auroral zone current in the 
sunlit hemisphere seems to become most 
intense a few tens of minutes later than 
the westward auroral zone current does. 
Change in the total amount of both the 
westward and eastward auroral zone 
currents with time is schematically 
shown in Fig. 5, where W and E denote 
the total amounts of the westward and 
eastward auroral zone currents in arbi- 


Fig. 5. Change in the total amount 
of the westward electric current in 


trary unit respectively. 


(iii) In the middle and low latitudes eT Nest per Pak tee ae 9 9 PRE 
the area in which positive bay is sphere (W) and the eastward one 
observed, becomes wider with time in in the sunlit hemisphere (E). 


the above individual examples of bays. . 
A kinematical interpretation for this phenomenon was already given in Chapter 
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II, §3. It was already shown in Chapter III, §3, that the positive bay area is 
narrower than the negative bay area when the current intensity of westward 
auroral zone is larger than that of the eastward one, as will be seen in Fig. 
9 there. The broadening of positive bay area:during the whole course of bays 
will be interpreted by the general tendency for the mutual relation between 
the change in the intensity of westward and eastward auroral zone currents 
illustrated in Fig. 5 here. : 

For the purpose of establishing the above-mentioned charactiristics of the 
progressive change in the current-system of geomagnetic bays, further exami- 
nation by means of a number of bays will be necessary, especially with respect 
to the disturbance field along the auroral zone and its neighbourhood. From 
the standpoint of the dynamo theory, the progressive change in the disturbing 
force of bays along the auroral zone means that the increase or decrease in 
the electrical conductivity with time does not occur simultaneously there. The 
anomalous increase in the electrical conductivity of the upper atmosphere in 
the auroral zone is considered to be caused by the intense ionization there by 
impinging protons, which approach along lines of the earth’s magnetic force 
[14-17]. Then, it is suggested here that the impinging particles do not always 
show uniform concentration along the auroral zone. 
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Chapter V. Constitution of Polar Magnetic Storms 


§1. A problem on the constitution of polar magnetic disturbances 


It has already been established from the statistical study on the disturbance 
field that the polar magnetic storms, in their average, have generally such the 
mode as that named Sp-field by S. Chapman and E.H. Vestine [37, 49, 507. 
It was also pointed out by M. Hasegawa that the said Sp-field exists even on 
quiet days, though its intensity is far feeble compared with that on disturbed 
days [10]. These conclusions of the analysis of Sp-field concern its average 
state. Being based upon the conception of the average Sp-field, a few theories 
of its formation have been proposed: they are the theories by H. Alfven [20], 
by D.F. Martyn [21] and the dynamo theory by T. Rikitake [26] and the 
present writer (42, 43], which were already mentioned in Chapters I and III. 

On the other hand, it has also been well known that the intensity of the 
disturbance field of polar magnetic storms changes markedly time by time, 
especially in the polar regions. T. Nagata [27] and E. H. Vestine [28] analysed 
the world-wide hourly mean aspects of the disturbance field at the time of 
some typical magnetic storms during the Second International Polar Year with 
the aid of hourly mean or instantaneous values over the world, and pointed 
out that hourly changes in the polar part of magnetic storms are generally 
complicated, though magnetic storms seem to be composed approximately of 
Dst- and Sp-components, the mode of the latter being subjected to the local 
time or the position of the sun. 

One may notice at a glance of magnetograms recorded at high latitude 
stations that the geomagnetic variation of fairly large amplitude, whose dura- 
tion is from several minutes to a few hours, superposes upon the general 
diurnal variation of the magnetic field, as will be seen later in several ex- 
amples. It can be considered that the so-called Sp-field of magnetic storms is 
composed of a number of elementary disturbances which take place intermit- 
tently or successively with duration from several minutes to a few hours, and 
Sp-field reveals itself as the average of these elementary disturbances. An 
important problem here will be whether or not each one of those elementary 
disturbances has also nearly the same mode as that of average Sp-field. This 
problem ought to be directly related to the physical mechanism of polar mag- 
netic storms. A.G. McNish and H.F. Johnston have already proposed this 
question in a different expression, and they pointed out an example of a severe 
magnetic storm where Sp-field in the middle and low latitudes was not stable 
with respect to its mode [51]. The study in the followings is also an answer to 
the above-mentioned question. A preliminary report of this study was recently 
published under joint authorship of Prof. T. Nagata and the present writer [52]. 
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For the present study, magnetic storms occurred during the Second Polar 
Year were analysed, because a number of magnetogram copies and hourly 
values were available for the storms in the period over the world. Through- 
out the Second Polar Year, however, the geomagnetic activity was generally 
quiet, so that only seven magnetic storms could be picked out as typical ex- 
amples. The dates of those storms are 


Oct. 14-15, 1932, May 1-2, 1933, 
Feb. 19-20, 1933, July 23, 1933, 

Feb. 23, 1933, July 24, 1933. 

Apr. 30, 1933, 


Among them the magnetic storm on May 1-2, 1933, was the severest one, 
Strictly speaking, the hourly values for the geomagnetic elements are not 
sufficient for deriving the real pattern of the progressive change of polar 
magnetic storms. In this study, photostat copies of magnetograms at the time 
of the magnetic storms on April 30 through May 2, 1933, at 33 observatories 
over the northern hemisphere of the world, which were prepared by Prof. T. 
Nagata at Department of Terrestrial Magnetism, Carnegie Institution of Wash- 
ington, were chiefly examined, while the hourly mean or instantaneous values 
were utilized only for getting an outline of characteristics of successive changes 
in polar magnetic storms and their behaviour. In Fig. 1, the distribution of 
observatories is given, where full circles mean the stations at which photostat 


WORLD, NORTHERN HEMISPHERE IN GEOMAGNETIC COORDINATE 


Fig. 1. Distribution of the magnetic observatories, whose data 
were used in the present study. 


tion of north geomagnetic pole 
assumes 9=78°,5N, =291°.0 E), 


-G 


Geomagnetic longitude, 
Geographic latitude, 


2 Geographic longitude, 
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Table I. List of magnetic observatories. 

Observatory Abbreviation © A ? 2, wv D |D-» 
Thule TH 88.0 | 0.0| 765) 2911] 0.0|—81.2|—81.2 
Godhayn GO 79.8 | 32.5) 69.2) 306.5 |—17.5 |—57.2 |—39.7 
Scoresby Sund SS 75.8 | 81.8] 70.5| 338.0 |—36.2 |—34.5 1.7 
Angmagssalik AN 74.3 | 52.7 | 65.6) 37.6 |—22.6|—39.8 |—17.2 
Sveagruvan SV 73.9 |130.7 | 77.9| 16.8|—46.2|— 4.9] 41.3 
Calm Bay CB fAlasy \\alsyes' 80/3) 52.8 |= 32:2) 211 53:3 
Bjornoja BJ 71.0 /124.7| 74.5] 19.2|-379|— 19] 360 
Julianehaab JU 70.8 | 35.6 | 60.7} 314.0 |—13.8 |—43.5 |— 29.7 
Fort Rae FR 69.0 |290.9 | 62.8) 243.9| 24.1) 37.5] 13.4 
Tromso TR 67.1 |116.7 | 69.7} 18.9|—20.8/— 3.7} 27.1 
Petsamo PE 64.9 |125.8] 69.5} 312|-27.6| 58] 33.4 
Matotchkin Shar MS €4.8 |146.5 | 73.3| 56.4|—22.4| 21.7] 44.1 
College-Fairbanks CO 64.5 | 255.4 | 64.9) 212.2} 27.0] 30.4 3.4 
Sodankyla SO 63.8 |120.0 | 67.4) 26.6 |— 26.7 3.0)| 29.7 
Dickson DI 63.0 |161.5 | 73.5] 80.4|—12.8| 28.5) 41.3 
Lerwick LE €2.5 | 88.6 | €0.1| 358.8 |—23.6 |—13.6; 10.0 
Dombas DO 62.3 |100.0 | 62.1 OF 23/0) (O.olenlosd 
Meanook ME 61.8 |301.0 | 54.6] 246.7] 17.2} 264] 9.2 
Sitka SI 60.0 | 275.4 | 57.0] 224.7) 21.4) 30.1 8.7 
(Eskdalemuir) ES 58.5 | 82.9 | 55.3| 356.8|-20.4/—14.2| 6.2 
Lovo LO 58.1 1/1058 | 59.4) 17.8/—22.1|— 2.5] 19.6 
Sloutzk Sil, 56.0 |117.0 | 59.7}; 30.5|—20.6 4.4| 25.0 
(Rude Skov) RS 55.8 | £85) 55.8} 12.4)—20.6|— 5.5} 15.1 
Agincourt AG 55.0 | 347.0 | 43.8] 280.7} 3.6|— 7.6|—11.2 
Abinger AB 54.0 | 83.3 | 51.2} 359.6/—18.4|—11.8| 6.6 
Val Joyeux VJ 51eS 04.0 48.8) em 210i ail ei awl Oil em reO 
Yakutsk YA 51.0 |193.8 | 62.0} 129.7 5.8 |—16.2 |—22.0 
(Swider) SW 50.6 |1046 | 52.1) 21.2/-183|— 1.5] 16.8 
Cheltenham CH 50.1 | 350.5 | 38.7] 283.2) 24|— 7.1/— 9.5 
San Miguel SM 45.6 | 50.9 | 37.8) 334.4 |—11.3 |—18.2|— 6.9 
(Ebro) EB 43.9 | 79.7| 40.8) 0.5|—15.0/- 9.9} 51 
(San Fernando) SF 41.0 | 71.3 | 36.5} 353.8 |—13.6 |—12.2 1.4 
(Zouy ) ZO 410 |174.4 | 52.5] 104.0 |— 1.8]-- 0.1 1.7 
Tucson TU 40.4 |312.2 | 32.2) 249.2} 10.1|) 13.9 3.8" 
Toyohara Eleva 36.9 |203.5 | 47.0} 142.8 7 |— 9.0 |—15.7 
San Juan SJ 29.9 3:2) 18.4)" 293:9|— 0.7 |— 5.2 4:5 
(Teoloyucan) TE 29'6 9327.0) || 19:8) 260:8)) = 16.6)| 9.5) ae 2:9 
(Helwan) HE 27.2 \106.4 | 29.9) 31.3 |—12.7 (OO) SA 
Kakioka KA 26.0 |206.0 | 36.2] 140.2} 62/— 5.8|—12.0 
Honolulu HO Qt (266.5 eee lei) 2019 Le oil Org 222 
Huancayo HU —0.6 | 353.8 |—12.0| 284.7; 1.3) 7.3] 6.0 

“Notes: & ‘Geomagnetic latitude (The posi- v The angle formed by the great 


circle joining the station and 
the geomagnetic pole with the 


geographical 


station (eastward positive), 


meridian of the 


Declination (for April, 1933). 
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copies of magnetograms were available in the present study, and hollow circles 
mean those, at which only hourly mean or instantaneous values were used for 
reference. The list of observatories is given in Table I, where the abbreviation, 
position and other elements of the observatories are summarized, and the sta- 
tions corresponding to hollow circles are written in parentheses. 


§2. Geomagnetic activity from April 30 through May 2, 1933 


K>z-indices after J. Bartels during April 30 May 2, 1933, are as follows [53]. 
The magnetic storm on April 30, 1933, beginning with a sudden commencement 


Date 0-3h 3-6h 6-9h 9-12h 12-15h 15-18h 18-21h 21-24h | Daily sum 
April 30, 1933 TS Pon te 1g ies Ag 5s ao AE 19_ 
May 1, 1933 aie) ZW ks 30 70 8_ 8 fe 400 
May 2, 1933 aie Ee Bh) She Ls: Da 2 30 22 


at 16h 28m G.M.T. (preceded by preliminary reverse impulse in some locali- 
ties [54]) seems to diminish its activity once about Ih on May 1, after a 
considerable activity during 20-24h G.M.T. on April 39. It seems likely that 
another disturbance took place from about 12h on May 1 to about 3h on May 
2. The magnetic activity during 15-21h on May 1 was the most severest one 
throughout the Second International Polar Year. According to T. Nagata’s or 
E. H. Vestine’s analysis, Ds:-field revealed its maximum negative value in the 
horizontal intensity of magnetic field near the equator, amounting to about 
100 y around 18h G.M.T. on the day [27, 28]. After 3h G.M.T. on May 2, 
Kp-indices in the above table show quite small magnitudes, indicating that 


only minor disturbances took place in the recovery stage of the above magnetic 
storms. 


§3. Equivalent current system of each steps of polar magnetic storm com- 
ponent of the storm on April 30, 1933 


As will be shown in Figs. 2 and 3 for example, the magnetic storm on 
April 30, 1933 was composed of a number of elementary disturbances, each of 
which continued from about 10 minutes to an hour and took place successively 
or intermittently. 

The geomagnetic variation field at any instant can be considered to be 
composed of the field of the quiet day’s variation (Sq) and the disturbance 
field (D). And further, the disturbance field D is composed of the polar storm 
component (Dp) and the non-polar component (Dy), the latter usually attri- 
buted to the magnetic field caused by the idealized equatorial ring current, as 
described in Chapter I. 


In the auroral and polar cap regions, however, the polar storm component 
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Dp, is extremely larger than 
the components of Sz and Dy 
as will be seen in Figs. 2 and 
3, So that the latter two may 
be practically neglected. In the 
present study, therefore, the 
polar storm component is de- 
fined to be the deviation of 


22 
= est tye ws n 


geomagnetic force from the 16 April ae 20 0 May 1 2 
straight line connecting the 
beginning and ending of indi- Fig. 2. Magnetogram registered at Sodankyla 


vidual elementary disturbances on April 30, 1932. 
or of groups of elementary dis- 
turbanees, in such a way as 
illustrated by broken lines in 
Figs. 2 and 3. The above de- 
finition of the polar storm com- 
ponent, which seems to be quite 
safe for the polar region data, 
was extended for the data in 
the middle and low lJati- 
tudes, where the times of 
beginning and ending of each 
elementary disturbance were 
common throughout all simul- 
taneous data over the world, 
being determined from _ the 


polar region data only. April 30, 1933 May | 
As typical examples of in- Fig. 3. Magnetogram registered at Matotchkin 
dividual elementary polar dis- Shar on April 30, 1933. 


turbances, five cases indicated 

by arrows in Figs. 2 and 3 were selected. The selected times corresponding 
to the maxima of those elementary disturbances are 17h 10m, 20h 52m, 21h 
15m, 21h 45m and 23h 15m G.M.T. Current arrows and equivalent. overhead 
current-systems for the estimated disturbing force at that times over the world 
are shown in Figs. 4-8, where 2/3 of observed horizontal components and twice 
of observed vertical component’ were taken as the external origin part of 


respective components. 
The equivalent current-systems of Dp-field illustrated here can be clas- 


sified into two groups: i.e. (A) those almost similar to the equivalent current- 
system of the average Sp-field or of the average bay-field, and (B) those 
markedly different from the average Sp-field. In Case (A), both the westward 
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and eastward auroral zone currents in the dark and sunlit hemispheres res- 
pectively appear definitely, while in most cases of Case (B), the westward 
auroral zone current only appears, or the westward current is much predomi- 


nant, the eastward current being very weak. 

In other words, the intense auroral zone current of individual elementary 
disturbances of polar magnetic storms is frequently localized around a certain 
locality in the auroral zone, though there are also such the cases where the 
distribution of the auroral zone current is not much different from that of 


Fig. 4. Equivalent current-system of Fig. 5. Equivalent current-system of 
elementary polar magnetic disturbance elementary polar magnetic disturbance 
at 17h 10m G.M.T. on Apr. 30, 1933. at 20h 52m G.M.T. on Apr. 30, 1933. 


Fig. 6. Equivalent current-system of 
elementary polar magnetic disturbance 


Fig. 7. Equivalent current-system of 
elementary polar magnetic disturbance 
at 21h 15m G.M.T. on Apr. 30, 1933. at 21h 45m G.M.T. on Apr. 30, 1933. 
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Fig. 8. Equivalent current-system of Fig. 9. Electric current-system produced 
elementary polar magnetic disturbance by an electric doublet on a -spherical 
at 23h 15m G.M.T. on Apr. 30, 1933. surface. : 


average Sp-field. It will be noticed here that the simplest pattern of equivalent 
current-system of elementary polar disturbance such as shown in Figs. 5-7 is 
almost similar to the electric current-system on a conducting spherical surface 
caused by an electric doublet situated on the auroral zone and directed west- 
ward, such as illustrated in Fig. 9. The theoretical calculation of the current- 
system in Fig. 9 will be mentioned later in § 7. 


§4. Disturbance field at the severest stage of the magnetic storm on May 
1, 1933 


Some examples of magnetograms during the severe magnetic disturbance 
on May 1, 1933, are reproduced in Fig. 10. It may be allowed to consider, for 
convenience, that variations of less than a few hours in their period superpose 
upon the general daily variation at every station. There are some superposing 
variations of duration from a few tens of minutes to about an hour, which 
should be considered to be far beyond limits of fluctuations in the general 
conception. These kinds of superposing variations were picked out on the mag- 
netogram copies. Definition of the superposing variation is as follows. When 
a peak variation of fairly large amplitude is recorded almost simultaneously 
at some observatories, beginning and ending times of the variation are deter- 
mined from these records. The deviation from the straight line connecting 
the values at beginning and ending times, which is illustrated by broken lines 
in Fig. 10, is defined to be the disturbance force of the superposing variation. 
The same definition of the disturbance force is adopted at every other station, 
where the beginning and ending times are common to all stations. 
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/933 May 1 


Fort Rae 


16 18 - 20 22 2h 16 18 22h 


Fig. 10. Examples of magnetic records during the severe 
disturbance on May 1, 1933. 


Four examples of this kind of superposing variations are illustrated in 
Figs. 11-14. The selected ones are those centred at 15h 42m, 16h 50m, 20h 15m 
and 22h 15m G.M.T. on May 1. The overhead current arrows corresponding 
to the maxima of the disturbance field of the superposing variations at above 
respective times are illustrated in the figures, while current-system for the 
hourly mean disturbance field for 16h G.M.T. by E.H. Vestine [28] is given 
in Figs. 11 and 12 for comparison, with its phase with respect to the earth ad-. 
vanced by 18 minutes and delayed by 50 minutes respectively. In Figs. 13 
and 14, current-systems for the hourly mean disturbance field around 20h and 
22h G.M.T. are illustrated with similar corrections for their phase as in the 
above cases. The correction for induced current within the earth was per- 
formed by assuming that the ratio of the external origin part of the horizontal 
disturbing field to the observed one on the earth’s surface is 0.9 in the polar cap, 
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Fig. 11. Disturbance field at 15h 42m Fig. 12. Disturbance field at 16h 50m 
G.M.T. on May 1, 1933, represented G.M.T. on May 1, 1933, represented 
by the combination of current-system by the combination of current-system 
and current arrows. and current arrows. 


Fig. 13. Disturbance field at 20h 15m Fig. 14. Disturbance field at 22h 15m 
G.M. T. on May 1, 1933, represented G.M.T. on May 1, 1933, represented 
by the combination of current-system by the combination of current-system 


and current arrows. and current arrows. 


0.7 in the auroral zone, and 0.6 outside the auroral zone in Vestine’s calcula- 
tion, while the ratio adopted by us throughout the previous and present studies 
is 2/3 all over the world. Although the definition of the disturbance field and 
the method of correction for induced current by him have some difference with 
those adopted here by us, general aspects of the current-system for the dis- 
turbance field at that time will never be altered so much qualitatively, since 
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the disturbance field itself is so large that the amplitude of the normal diurnal 
variation on quiet days is practically negligibly small compared with the 
magnitude of the disturbance field. Stream lines in the figures are so drawn 
that electric current of 2x105 amp. flows between successive stream lines. In 
- Figs. 11-14, we can clearly see that the current arrows for superposing varia- 
tions are not always parallel to the stream lines in the figures. 

The disturbance field of polar magnetic storms at an instant, Dp, will 
be expressed as 


Dp=Dpm+Dpi, o>) 


where Dp,, means the smoothed value of the polar disturbance field over some 
interval of time, while Dp; denotes the superposing variation of rather short 
duration upon it. In the cases of Figs. 11-14, Dp,, can be considered as the 
disturbance field represented by the distribution of stream lines, namely nearly 
perfect form of Sp-field in these cases, while Dp; the superposing variations 
given by the current arrows in the figures, as a first approximation. Strictly 
speaking, however, the smoothed disturbance is affected by the superposing 
variations when the magnitude of the latter is comparable or larger than that 
of the former. In such a region where the above condition of the disturbance 
field is revealed, then, the total disturbance field may not be considered as the 
superpoSition of disturbances represented by current lines and current arrows 
in the figures. In the above four examples, the current-system for Dp;, which 
is formed by the distribution of current arrows in the figures, will be nearly 
the same as the “doublet-field” type pattern of a simple elementary distur- 
bance. 

The superposing variations in the above examples are accompanied by the 


additional auroral zone current, which 


8 
& 66 ern: is remarked by long westward cur- 
3 rent arrows in the dark hemisphere 
S65 along the auroral zone. The position 
S 


of the centre of the additional auro- 
ral zone current of Dp; can be rough- 
ly estimated by taking also the simul- 


12 4 16 IS 20 22 oth 
May 1, 1933 


Fig. 15. Roughly estimated geomagnetic 


latitude of the centre of 
auroral zone current. 

(Line in the figure is the centre of the 
auroral zone current in European side 
after T. Nagata.) 


additional 


taneous change in the vertical com- 
ponent of geomagnetic field into ac- 
count. Fig. 15 illustrates the roughly 
estimated geomagnetic latitude of the 
centering position of the additional 
auroral zone current of. superposing 
variations. Large circles in the figure 
are of the four examples shown in 


Figs. 11-14, and small ones of the remaining examples, which are not illustrated 
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by figures. The centre of the intense auroral zone current at that time in ~ 
the European region determined by T. Nagata [27] after the detailed analysis 
of hourly mean values is also shown in the figure for comparison. According 
to his result [44], the auroral zone current flows some 10 degrees of latitude 
in its width, so that the centre of the additional auroral zone current is situ- 
ated within the width of the general auroral zone current. 

The four examples of the superposing variation shown here are the most 
distinct ones among many examples examined. Remaining examples show 
more or less complicated distribution of current arrows over the world, indi- 
cating that Dp; at those times will be expressed by some complicated combina- 
tion of elementary polar disturbances, whose origins are distributed at more 


than two points along the auroral zone simultaneously, with their magnitude 
different with each other. 


§5. Examples of geomagnetic variations in the later stage of the magnetic 
storm 


In Figs. 16-21, hourly mean departures of the magnetic field during 4-10h 
on May 2, 1933, from the average of monthly mean values of April and May 
in 1933, are represented by the distribution of equivalent overhead current 
arrows. The westward current arrows in the low latitudes may indicate the 
presence of the equatorial ring current in the later stage of the magnetic 
storm. Although K>-indices during the above period shows shall numbers of 
scale value as mentioned in § 2, magnetic records registered at the stations in 


ETIC COORDINATE 
WORLD, NORTHERN HEMISPHERE 1N GEOMAGNE WORLD, NORTHERN HEMISPHERE IN GEOMAGNETIC COORDINATE 


Fig. 16. Distribution of overhead current Fig. 17. Distribution of overhead current 
arrows corresponding to the hourly arrows corresponding to the hourly 
mean disturbance field at 4-5h G. M. T. mean disturbance field at 5-6h G. M. T. 
on May 2, 1933. on May 2, 1933. 
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Fig. 18. Distribution of overhead current Fig. 19. Distribution of overhead current 
arrows corresponding to the hourly arrows corresponding to the hourly 
mean disturbance field at 6-7h G. M. T. mean disturbance field at 7-8h G. M. T. 


on May 2, 1933. on May 2, 1933. 


Fig. 20. Distribution of overhead current Fig. 21. Distribution of overhead current 


arrows corresponding to the hourly arrows corresponding to the hourly 
mean disturbance field at 8-9h G. M. T. mean disturbance field at 9-10h G.M. T. 
on May 2, 1933. on May 2, 1933. 


the north of American continent during 6-10h G.M.T. on May 2 showa 
considerable activity of the geomagnetic field, while only slight variations are 
observed even around the Arctic ocean. Some examples of magnetograms are 
given in Fig. 22. 


It is clearly perceived that some variations superpose upon the general 
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Fig. 22. Examples of superposing variations upon the recovery 
curve of magnetograms on May 2, 1933. 


recovery curve of the magnetic field after the preceding severe activity. Equiva- 
lent overhead current arrows corresponding to the above-mentioned super- 
posing variations centred at 06h 30m and 09h 40m G.M.T. on May 2, are 
derived by subtracting the, general recovery curve, which is shown in Fig. 22 
with broken lines. Equivalent overhead current-systems are roughly estimated 
from the distribution of current arrows corresponding to the superposing 
variations at above times, and are given in Fig. 23 and Fig. 24 respectively, 
where electric current of 5x104 amp. flows between successive stream lines. 
Tt is illustrated here that the selected superposing variations at the time of 
the recovery stage of the magnetic storm are of the form of “ doublet-field ” 
type pattern. With the aid of the expression for the polar disturbance field 
as given by (1) in the preceding paragraph, it will be said that Dp; in the 
later stage of the magnetic storm in these cases can be also represented by 
the disturbance field of an elementary disturbance. In the present cases, Dpy, 
is not revealed as a definite form such as the average Sp-field, but shows some 
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Fig. 23. Current arrows and current- Fig. 24. Current arrows and current- 
system for the disturbance field of system for the disturbance field of 
superposing variation at 06h 30m superposing variation at 09h 40m 
G.M.T. on May 2, 1933. G. M. T. on May 2, 1933. 


complicated aspect, while Dp, in the preceding paragraph shows almost per- 
fect form of Sp-field. 


§6. Elementary disturbances and average Sp-field 


The hourly distribution maps of the disturbance field over the world 
during principal magnetic storms in the Second International Polar Year were 
obtained by T. Nagata, the number of hourly maps amounting to 216 in total, 
where the disturbance field was defined to be the deviation of geomagnetic 
forces at disturbed times from thos? at the corresponding times of median 
values on quiet days in the same month, as a first approximation. Throughout 
216 hours data of disturbances, the nearly perfect doublet-field type appeared 
in rather few cases, so far as concerning the hourly mean values. Only two 
examples are illustrated in Figs. 25 and 26, where the electric current of 105. 
amp. flows between successive stream lines. In many cases among the hourly 
maps, the westward auroral zone current took place strongly in the dark 
hemisphere, while a weak eastward auroral zone current appeared in the sunlit 
hemisphere. Sometimes, the latter being negligibly weak in intensity compared 
with the former, as will be seen in Figs. 25 and 26. 

This result will be partly due to that the resolving power of “ hourly 
mean value” is not enough with respect to time to pick out a single pattern 
of an elementary polar disturbance, and partly due to the definition of distur- 
bance field adopted in this paragraph. In short, the hourly mean aspect is 
the mean of varying aspects during an hour, as it is. 
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Fig. 25. Example of hourly mean equiva- Fig. 26. Example of hourly mean equiva- 
lent current-system of polar magnetic lent current-system of polar magnetic 
storm, April 30, 21-22h G.M. T., 1933. storm, February 19, 23-24h G.M.T., 1933. 
(In Figs. 25 and 26, electric current of (after T. Nagata) 
10° amp. flows between successive 
stream lines.) (after T. Nagata) 


On the other hand, Fig. 27 shows the average of hourly mean current- 
systems of intense parts of five polar magnetic storms. In the above-mentioned 
averaging, the data of the magnetic storm on May 1, 1933 were excluded, 
since the intense polar storm in this case was of the form of nearly perfect- 
Sp-field throughout the whole storm time, as already mentioned before. The 
hourly aspects used for deriving the average current-system in Fig. 27 are of 


1932 Oct. 15, 09-10h, 13-14h, 16-17h, 23-24h, 
1933 Feb. 19, 14-15h, 19-20h, 23-24h, 

Feb. 23, 15-16h, 

Feb. 24, 01-02h, 06-07h, 

July 23, 02-03h, 13-14h, 

July 24, 01-02h, 07-08h, 12-13h G.M. T. 


As will be seen in Fig. 27, the average of hourly current-systems of polar 
magnetic disturbance is of almost perfect form of current-system for Sp-field. It 
will be concluded, therefore, that the average of a large number of individual 
elementary polar disturbances becomes the well-known pattern of average Sp- 
field, while each elementary polar disturbance is much “ partial ” the simplest 
form being of the doublet-field type. An example of the appearance of nearly 
doublet-field type pattern at the time of geomagnetic bay will be noticed in 
Fig. 4 (a) or (6) in Chapter IV. 
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Summarizing all results obtained 
in §§ 3-6, we can have an idealized 
picture for the composition of polar 
magnetic storms. That is to say, the 
most elemental unit of polar mag- 
netic storms will be the disturbance 
field of the doublet-field type. This 
Gz 7 : unit will correspond to the condition 
that a small area in the ionosphere 
over the auroral zone produces elec- 
tromotive force, caused by impinging 
of corpuscular stream. This argu- 
ment does not exclude such the con- 
dition that more than two elementary 
disturbances take place simultaneous- 


Fig. 27. The mean equivalent current- ly at different locations, or that the 
system for polar magnetic storms during 
the Second Polar Year. (Electric cur- 
rent between successive stream lines is 


area producing e.m.f. is much ex- 
tended along the auroral zone, fre- 


2.9104 amp.) (after T. Nagata) quently almost the whole auroral 
zone around the polar cap being ac- 
tivated. 


With the aid of mathematical expression, the constitution of polar magnetic 
storms will be written as follows. Let the disturbance field observed at a 
point (8, X%) on the earth’s surface due to an elementary polar disturbance of 
unit magnitude and of its origin at (@, Xo) in the auroral zone of polar dist- 
ance @) denote by Dh(O, A, 6, »), general polar disturbance field at the point 
at an instant, Dp(@, %), will be given by 


Dp(B, %)= ZB kA) - Db (0, , Oo, Xo) » (2) 


where k(A)) is a factor concerning the distribution of origin and magnitude 
of elementary disturbances along the auroral zone at the instant. D> (6, X, 9, 
Xo) will be derived from such the current-system as shown in Fig. 9 in § 3, 
and it contains 6) and X»y as parameters, which indicate the position of the 
centre of the elementary disturbance. In the above expression, however, mutual 
influence between more than two elementary disturbances, which take place 
simultaneously at different locations, is ignored as a first approximation. In 
practice, the distribution of origin and magnitude of elementary disturbances 


will vary time by time. When elementary disturbances are distributed uni- 
formly along the auroral zone, namely 


k (Ay) = const.=ky (3) 
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along the whole auroral zone, resultant disturbance field will be written as 


2a 
Dp(O, =ho| DCO, 2, Oy Xo) Arp « (4) 
0 


If we denote the perfect form of Sp-field by S},(@, \), the relation 
22 0 
Sb, %) = | DCO, X, 8, %0) dro (5) 
0 


seems to hold from the above-mentioned actual analysis of geomagnetic data. 
Further, if we consider 


R(X )=R(Ap)+SR(AQ) , (6) 


where k(A)) means the average value of (Aj) over some interval of time, 
and dk(Ay) is the deviation of the instantaneous value of k(Ao) from R(Ay) at 
the time, Dp, and Dp; defined by (1) will be replaced by 


Den, = = RA) + DBO, 2, 8 0) » (7) 
Dp(0, %)= > Sk(Ao) - DPA, 2, Oo, Ao) - (8) 


When (Ao) along the auroral zone shows nearly uniform distribution, then 
nearly perfect form of Sp-field is revealed. On the other hand, if dk(Xo) dis- 
tributes like a 6-function along the auroral zone, superposing variation of 
nearly doublet-field type pattern will appear at that time, as in the cases 
described in the preceding paragraphs. 

Mathematical foundation for the above conclusion from the standpoint of 
dynamo theory will be given in the following paragraphs. 


§7. Electric current-system on a conducting spherical surface by the dis- 
tribution of electric doublets 


In the preceding paragraphs, it was concluded that the simplest pattern of 
the polar magnetic disturbance is a “doublet-field” type pattern, and the 
electric current-system on a conducting spherical surface for the stationary 
state caused by an electric doublet situated westward on the auroral zone was 
shown in Fig. 9 in-§3 for comparison. Here, the method of deriving such a 
current system will be described. 

The electric potential ¢ and the current function J at a point P(y, r) in 
a two dimensional case when an electric doublet of its moment py is situated 
at the origin of the coordinate are given by 


cos »% 
o=— Fo F a 
fy Sin A 
:. 5 


ipas 
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The equation of the stream line is given by 


Se eons (10) 
r 


or =const. (11) 


sin X 

Here, we consider a sphere of radius R, which contacts at the origin of 
(x, y) or (7, X) plane with its top pole just similarly as shown by Fig. 17 in 
Chapter III. Let the intersection point on the surface of the sphere with the 
straight line connecting a point P(7, \) on the plane and the bottom pole of 
the sphere denote by Q. A point P is projected to a point Q on the spherical 
surface. The coordinate of Q is (R, 9, A) referring to the spherical coordinate 
with its origin ‘at the centre of the sphere. Then, the following relations hold, 


7y=2R tan + ; (12) 


A=A, (13) 
as will be seen from the comparison with Fig. 17 in Chapter III. Substituting 


the above relations into (9), the current function becomes 


J= Pik cot ~. sinA. (14) 


On the other hand, the electric potential is written as 


o=— ae cot cosA. (15) 
It is easily proved that 
7 ==()e 
rh (16) 
r?jJ=0, 
on the sphere, and moreover, the equation 
0 0 0 
J Op aj “o (17) 


R00 ROO Rsin9 dA RsinOdA 


is satisfied on the sphere, which shows the orthogonality of equipotential lines 
and stream lines. Thus, lines of J=const. are proved to express the stream 
lines on the spherical surface. 

In the next place, a new spherical coordinate system (CR, 8,2) is intro- 
duced, where the position of the electric doublet is CR, 8, 7) as is illustrated 
by Fig. 28. Mutual relations between two spherical coordinate systems are 

sin 0 cos A=—sin@ sind, 
sin 9 sin A=cos 6 sin @ cos X+sin 09 cos 0 2 (18) 


cos 9 =cos 0) Cos O)>—sin 6) sin @ cosa. 
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Fig. 28. (¢ and ® in the figure must be read as ) and 4 respectively.) 


From the above relations, 


O SipAcsesn OsinA ue Si Oo cos O+C08 8p sin 0 cos x - 9) 
1—cos 9 1—cos 6) cos 8+sin 6) sin 8 cosX 


cot 


Then, the current function J on the spherical surface of radius R under the 
presence of an electric doublet of moment py at the point (R, A, 7) is given 
by 


fz fod) sin @)cos @—cos G sin 8 cos xX 
2R 1—cos@ cos 6—sin @)sin @cosvr ’ 


(20) 


where oo denotes the electrical conductivity of the surface of the sphere. The 
distribution of stream lines is obtained by the numerical calculation by mean 
of (20), where 0) is taken appropriately. A case in which @)=22.5° is illus- 
trated as Fig. 29, which was already given by Fig. 9 in §3 in this Chapter. 
If the electric doublet is situated at CR, Go, Ao) along the latitude circle as in 
the above case, the current function becomes in general 


(eae sin 9) cos 8—cos Gp Sin 8 cos (A—Ay) (21) 
~ 2R  1—cos 6 cos 6—sin 6) Sin 9 cos (A—Ay) ~ 


If two electric doublets of the same moment fy are situated symmetrically with 
respect to 90°-270° meridian plane along the circle of colatitude 6, the cur- 
rent function in this case is 


jes sin 0p cos 9 +cos Oy Sin 9 cos A 
~ 2R \ 1i—cos @ cos +sin Oy sin 6 cos X 


sin 0) cos 8—cos 6 Sin 9 Cos A ! 
1—cos 6) cos@—sin 69 sin @ cos »® 
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_ ooo (cos O9>—cos @) sin 6 cos X : (22) 
~ R (cos @9—cos 0)2+sin? Ag sin? 6 sin?» 


The distribution of stream lines in the above case is shown in Fig. 30, which 


TRA I) 
SAY 


Fig. 29. The distribution of stream lines Fig. 30. The distribution of stream lines 
under the presence of an electric doub- under the presence of two electric doub- 
lets along the auroral zone. lets of the same moment symmetrically 

along the auroral zone. 


was utilized in the experimental study already mentioned in Chapter III, §4 
in order to test the uniformity of the distribution of electrical conductivity 
over the model conducting layer. The mode of current flow in the lower 
latitudes and near the pole is quite similar as that of the idealized current- 
system for Sp-field. 

Here we consider the spherical surface as a model of the conducting 
ionosphere over the earth. The current distribution in the conducting iono- 
sphere accompanies the magnetic field, which will be observed on the earth’s 
surface. The magnetic field produced by the electric current-system in the 
ionosphere observed at a point (@, X) on the earth’s surface of the middle and 
low latitudes can be roughly estimated as follows. If x- and y-axes are taken 
southward and eastward with its origin at the point on the earth’s surface, 
the x- and y-components of magnetic field Hz and Hy, produced by the over- 
head current-system are generally given by 


3 ay 3 as 
H-=— ie = 
ate R00’ Hy oo ee Oe Ga) 


BM: : 
where the factor “9 is attributed to the assumption that the induced current 


within the earth contributes half as much as the external current to the mag- 
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netic field observed on the earth’s surface. For example, in the case of the 
current-system of Fig. 29, the magnetic field produced by the current-system 
along 0° and 180° meridians, where Hy vanishes, will be given by 


~_. Pot 1 2 ee 
x 2R2  1—cos (646) along 0° meridian, 
(24) 
by Oo 1 


2R2 1—cos (@—0) along 180° meridian. 


In the case of the current-system of Fig. 30, the corresponding magnetic field 
is given by 


H.=— Hoo 1—cos 0) cos @ s ee 
x Rist Ceosbowcos 0% along 0° meridian, 
fy oy) 1—cOS Oy cos @ oi) 
Ax= along 180° meridian. 


R2  (cos@y—cos 0)? 


If we adopt R=6.47x108cm., namely the height of current-system assumes 
100 km. above the earth’s surface, g7==1x10-7 e.m.u. and p)=102 e.m.u., the 
calculated intensity of magnetic field due to the current-system on the spheri- 
cal surface can be obtained as summarized in Table II. 


Table II. The horizontal magnetic field (northward positive) produced 
by the current-systems of Figs. 29 and 30. 


In the case of Fig. 29 In the case of Fig. 30 
Colatitude 
0° meridian 180° meridian 0° meridian 180° meridian 

30° —28.97 13197 —13487 13487 
35 —24.4 476 — 500 500 
40 —20.9 244 — 264 264 
45 —18.3 148 — 166 166 
50 —16.1 99.8 — 116 116 
55 —14.4 72.0 — 86.4 86.4 
60 —13.0 54.6 — 67.6 67.6 
65 —118 43.0 —54.8 54.8 
70 —10.8 34.7 —45.5 45.5 
7s, —10.0 28.9 — 38.9 38.9 
80 —9.3 24.4 —33.7 Sonn 
85 —8,.7 20.9 —296 29.6 
90 —8.2 18.3 — 26.4 26.4 


§8. Current distribution over the conducting ionosphere caused by the dy- 
namo action with anomalous conductivity along the auroral zone 


As has been already mentioned, the idealized Sp current-system can be 
derived from the dynamo theoretical calculation when the electrical conductivity 
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over the auroral zone ionosphere is much large compared with that in the 
other regions. On the other hand, the examination of individual elementary 
disturbances reveals us that the idealized Sp current-system is the average 
of a large number of elementary polar disturbances, which show more “ partial” 
distribution of disturbing forces, and their simplest pattern is of a doublet- 
field type. Hence, it will be reasonable to consider that the elementary polar 
disturbance corresponds to the partial increase in the electrical conductivity of 
the upper atmosphere along the auroral zone due to the anomalous ionization 
there. Fig. 31 shows two extremum cases of the distribution of anomalously 
high ionization in the ionospheric region of the auroral zone, where the whole 
auroral zone is highly conductive in Case (A), while only a part of it has 
anomalously high conductivity in Case (B). Case (A) has already been theo- 
retically examined in detail in Chapter JJJ, as mentioned before. The result 
of applying dynamo action for Case (B) is summarized in the followings. 

Assuming first that the highly conductive area is very small compared 
with the whole surface of the earth, the problem will be treated as a two 
dimensional problem with the plane earth model, so far concerning the be- 
haviour around the conductive area. The small highly conductive region along 
the auroral zone assumes an elliptic area, major and minor radii of @ and b 
respectively. The electrical conductivity there being assumed to be kao, where 
oo denotes the conductivity of the spherical surface outside the elliptic area, 
and & is a constant factor larger than unity, as will be seen in Fig. 32. 

Since there is the ionospheric wind, its velocity being v(vz, vy), electromo- 
tive force CE) in the electrically conductive ionosphere produced by the dynamo 
action in the earth’s magnetic field (H ), and the resulting electric current (7), 
both in a steady state, are given by 


rot H=rot(v-A], (26) 

i=cE, (27) 
where the conductivity o is ko in the elliptic area and oy in the whole remain- 
ing part. If we assume that 

Vx=0, Vy=—0v, (28) 
the induced electromotive force is only in x-direction with its magnitude vHz, 
where Hz denotes the vertical intensity of the magnetic field directing into 
the paper in Fig. 32. Then, the problem will be reduced to obtain the current 


flow under the above distribution of electrical conductivity, when a uniform 


electric field E is affected in x-direction. We introduce the elliptic coordinates 
(é, 7), in the relation 


x=Icos n cosh E, y=lsin n sinh é. (29) 


The boundary of the elliptic area is given by é&=£) while a and b can be 
written by 
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(A) (B) 
Fig. 31. Highly conductive area in the auroral Fig. 32. An idealized model of ele- 
zone. mentary highly conductive area 
(A) Perfect ionization of the auroral zone, in the ionospheric region along 
(B) Partial ionization of the auroral zone. the auroral zone. 
a=I cosh &9, b=I/ sinh &o. (30) 


Complex electric potential X forthe uniform electric field outside the elliptic 
area is given by 
X=—Ez=—El(cos n cosh €+7sin 7 sinh €), (31) 
while the complex electric potential X; inside the elliptic area will be given by 
Xi=— Eiz =—Ejl(cos n cosh E+7 sin n sinh €), (32) 
where F; is the electric field inside the elliptic area, because the electric field 
there should be uniform also in x-direction. The electric field due to the 


polarization of electric charge at the boundary of the elliptic area is derived 
from the complex electric potential X’, the form of which is 


x — en = om e-(cos n—isin n), (33) 


where m denotes the electric moment directing also to x-direction. 
For the continuity of electric field and electric current along the boundary 
of the elliptic area, the relations : 


RCX)AERXIY=RCX;) , 
é TI MX)ABX') } Hho (Xi) 5 

must be satisfied at the boundary, or =&). From this, the following relations 

can be obtained : 


(34) 


14+ 

pit tanh £9 pO 
Ei-E yRtaoh ed |" py 4. eP) 

b 


cosh €) sinh £9 jp + (k—1) BG (36) 


toss to r¢ cae 
Ueto oach Eo +k sinh &, a+kb 
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In the case where the electric field E is in y-direction, the boundary conditions 


are 


SOX+IOXI=WX), ak 
o{ RCX)ERCX) § =Ra oh Xi) , 
instead of (34). E; and m in this case are derived as 
b 
Me eer 
1+ tanh Eo a 
E;=E = ; (38) 
; k+ tanh &o k+ oO 
a 
m= 5 (h—1) eo oh busin fo _ pF (g1) AGED B. .(39) 


sinh £9+k cosh &o ka+b 


In general, when k becomes larger than unity, an additional current flows 
circulating in the elliptic area and the surrounding area. The density of the 
additional current flow in the elliptic area I is given by 

I= kook j—ook = a) RE;—E) (40) 
in general. In the general case where the velocity components of the air mo- 
tion are vx and vy in x- and y-directions respectively, the components of the 
additional current density are written by means of (35), (38) and (40) as | 


+ (k=) 
oa —HzawWy meri tar WRT ’ 
ke 
(41) 
© (R=1) 
Iy=Hza vx b 
nae 


On the other hand, the components of electric moment due to polarization at 
the boundary of the elliptic area are 


cee 1 as Vy hing 1 
 ecacaate’ Hzab(a+b) (k—1) ser cay b(a+b) Ix, 
1! Ux 1 (4) 
MNy= a4 Hzab(a+b) (k—1) ka+b = ee a(at+b) Ty. 


Outside the elliptic area where e§e-£, the production of a highly conductive 
elliptic area results in the additional current flow, just same as that caused 
by an electric doublet at the centre of the elliptic area, its moment J being 
expressed by 


1 2 2 
Ih=(n?+-n2)t = — Hab ad-b)Ck—1 { v; ve \ 
o=( i= 4 HeaWa+b)(k-1) |S + ete} ee 


The mode of distribution of stream lines in a distant region from the elliptic 
area will be later mentioned about. 
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The distribution of electric current flow inside the elliptic area will be 


examined in the followings. Relative magnitude of J, and Ty varies with the 
values of a/b and k. If we put 


Vz=vVSINg, 
(44) 
Vy=—vCOS ~, 
I, and Iy become as follows from (41), 
§ (k-1) 5 (k-1) 
Ty = Hzaov = ei cos P =a Soa gees 
6 b a? 
* (k—1) al (k—1) 
Iy=Hzoyv b SB O= 90 spe Ey, 
k+ — k+ — 
a a 
because of 
Ex=Hzv cos @, 
x z Pp (46) 
Ey=Hzv sin Pp. 


av ites 
2 
values of Ix and Jy will be denoted by Iymax and Jy max- In Table III and Fig. 


33, relative magnitudes of [ymax and ZJymax under various values of a/b and k 


I, and Iy take their maximum values when y=0 and respectively, and the 


Table IIf. Relative magnitudes of Jymax and Iymax. 


k=5 k=10 k=30 k=100 ko 
a/b 
Temax Tymax Temax Tymax Temax Tymax Iemax lymax Jemax Tymax 
1 0.667 0.667 | 0818 0.818 | 0.936 0.936 | 0.980 0.980 | 1.000 1.000 
2 1.143 0.364 | 1500 0.429 | 1.812 0.476 | 1.940 0.492 | 2.000 0.500 
3 1.500 0.250 | 2.078 0.290 | 2.636 0.314 | 2.883 0.329 | 3.000 0.333 
5 2.000 0.154 | 3.000 0.176 | 4143 0.192 | 4.714 0.198 | 5.000 0.200 


are summarized, which can be obtained from the relation (45). In Fig. 34, 
relative magnitudes of J, and Jy are illustrated for the case of #=30. Rela- 
tions of J, and Iy to m are almost the same as Fig. 34 so far as & is larger 
than 10, as will be noticed in Table IJ and Fig. 33. This result indicates 
that the current density along the major axis of the ellipse is practically much 
predominant compared with that along the minor axis. If we put 


ae eqitan gy, 
4 (47) 


SSeS 
ie =tany, 
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Arbitrary unit 
Arbitrary unit 


Fig. 33. Change in ymax and Iymax with Fig. 34. Intensity of electric current 
the value of # for the cases of a/b=2, produced by the dynamo action 
Sando: ; within highly conductive area of an 


elliptic form, in the case of k=30. 


the relation between gp and ¥ can be derived by means of (45), namely 
b a 


a b b bk+a 


tan p= —— 
b a a-ak+b 


tan p= tang. (48) 


The calculated result of the relation between gy and wW for various cases is. 
given in Table IV. The result for the case of k=30 is illustrated in Fig. 35. 


Table IV. Relation between ? and wv. 


bas k=10 k=30 k=100 

Reon a ee 2. RR 5 goneek gs. roa o  regek 5 
0° | 0° 0? 0° 0° 60° = «0? 0° 09° 0° 0° 860° ~— 0? 
io"! 32°" 17 08 29 14 06 26°12 "os 26 ‘12'- 04 
20°| 72 35 “LT 59 29 12 5594 To 53 34° 09 
4h"! 104 55 °° 25 94 ‘46°. 1'8 86 39° 15 83 38 14 


40 150 80 37 3 SOL (anamenS 12:4, S169 92:2 12:0). ooo 2.0 
“50 2018) 1122) "5:2 leisy SoG ED TA = 80" 230 HOSS 20. (aeee co 
60 Con los 4.6 264 136 5.8 245° 115 46 23.7° 118 °742 
65 34.4 °197 94 S15 16778 ZO EE ONT. 286 13.8 51 
70 41.2. 24.6 11.9 SB.te 210 9 SoS MIG rerye3 34.9 17.4. 6.6.) 
75 499° 31.9 .16.0 ABOU 2512.4 44.4 23.7 9.7 43.4° 23.1 °89 
80 61.0 43.4 23.6 58.3 38.4 18.5 56.1 33.7 aA os Saizider 
85 74.6 62.3 41.3 73.0 58.0 33.9 (16 534 27.9 71.0 52.5 25:6 
87 80!6 72!5 55.7 796 69.4 483 78.7 66.0 41.5 78.3 65.3 38.6 
89 86.9 84.0 77.2 86.5 829 73.5 86.2 81.6 69.4 86.1 813 67.4 
90 90 90 90 90 90 90 90 90 90 90 90 90 
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These results also indicate the predominance of electric current flow along the 
major axis of the ellipse. In short, the electric current is susceptible to flow 
along the major axis of the elliptic area. In the practical geomagnetic problem, 
the above result suggests that the auroral zone current flows along the zone 
in general, so far as the highly conductive area is produced as an elliptic area 
with its major axis along the auroral zone. 

The total amount of the auroral zone current will be estimated to be 
5x105~8x105 amp. for the actual analysis of geomagnetic data, when the 
current intensity is fairly large, as will be seen in Figs. 4-8 in the preceding 
paragraph or other examples. On the other hand, assuming that v,x=0 and 


Vvy=—v in (41) for example, the total current intensity inside the conductive 
elliptic area, I, becomes 


Petite 2 (49) 

k+ 2 

b 

Putting Hz=0.5 e.m.u., op=10-7 e.m.u, v=10! cm/sec. and k=30, for example, 
we can get the following values of J. Then, it may be said that the total 
electric current of 5x105~8x10° amp. can be reasonably produced, provided 
that k& is fairly large, amounting to 10~100, and that ¢/b amounts to 2~5 as 


Table V. The total electric current in the conductive 
elliptic area for various cases. 


a/b b=100 km. b=200 km. b=300 km. b=400 km. 
2 1.8x 105 amp. 3.6 x 105 amp. 5.4 105 amp. 7.3 x 105 amp. 
3 2.6 x 105 5.3 x 105 7.9x 105 10.5 x 105 


5 4.1 x 105 8.3 x 105 12.4 x 105 


16.6 x 105 


will be presumed in practical cases as probable values. In the above Table 
V, 2b corresponds to the width of the auroral current zone. The moment of 
electric doublet corresponding to the above values of the auroral zone current 
is shown in Table VI. 


Table VI. The electric doublet moment produced by the presence of 
highly conductive area of various size. 


a/b b=100 km. b=200 km. 6b=300 km. b=400 km. 
2, 6.8 x 1017 e.m.u. 2.7 1018 e.m.u. 4.9x 1018 e.m.u. 8.7 x 1018 e.m.u. 
34 1.3 « 1018 5.3 x 1018 95x 1018 1.7 x 1019 


5 3.1 x 1018 1.2x 1019 2.2 x 1019 


4.0 x 1019 
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— 


Fig. 35. Relation between ? and Ww for Fig. 36. Resultant electric current-system 
the case of k=30. caused by superposition of a large num- 
ber of elementary current-systems of 

doublet-field type. 


Concerning the air motion in the upper atmosphere, it must be noticed 
that the ionospheric wind along the major axis of the ellipse or along the 
auroral zone is appreciably ineffective for producing the electric current com- 
pared with that along the minor axis, because in can be seen from (45) that 


1 § (aad OF 


i ghia re (50) 


when k»1. The above characteristic has been already mentioned about in 
the equations (42), (45), (47) and (48). Ionospheric wind derived from its 
velocity potential of the form of P}(cos@) cos’ is nearly uniform in the polar 
region. Summarizing abobe-mentioned results, we can get approximately the 
mean current function (J) of a large number of elementary polar disturbances 
uniformly distributed along the auroral zone as follows, 


= \a= laa Sin @) cos A—cos Gy sin 8 cos (Ayp—2) Rsin Od» 


A Byes SA A 1—cos 6) cos 8—sin 0 sin 8 cos (Ay—A) 


__ THoa0(COS Ay—cos 0) |a= {1 __ sin? 69 sin 0 pa ee (51) 


sin 6) sin 0 . (1—cos 89 cos 6): 


The distribution of stream lines given by the above mean current function J 
is illustrated in Fig. 36. The current-system in Fig. 36 is almost equal to the 
idealized current-system for Sp-field calculated from the dynamo theory under 
the assumption that the electrical conductivity in the auroral zone is uniformly 
large compared with that in the other regions. Thus, we can conclude that 
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a summation of a large number of elementary polar disturbances of the 
doublet-field type almost coincides with the so-called Sp-field. 

The mode of current flow in the lower latitudes produced by the dynamo- 
action in a small elliptic area along the auroral zone will be finally considered. 
In a two dimensional problem, the distribution of stream lines in a distant 
space from the elliptic area is approximated by the distribution of stream lines 
when an electric doublet is situated at the centre of the ellipse. The mode of 
current flow there should be symmetric with respect to the line AB in’ Fig. 37, 
whose equation is expressed by 

My 
I= q, * = tana x Cy 
referring to x- and y-axes in the figure. In the above expression, mz and my 
are x- and y-components of the electric moment due to polarization in the elliptic 
area, Directions of the air motion and the induced electromotive force are 
given in the figure by v and £ respectively. When the stream lines on the 
plane are projected on the spherical surface, in the way as mentioned before, 
the straight line AB in Fig. 37 will be projected as shown in Fig. 38, which is 


Fig. 37 


the view from above the pole. The circle of radius 7 in Fig. 37 is projected 
as the circle of polar distance 6) in Fig. 38. The line APB intersects the 
the equatorial plane at A and B, respectively. With the aid of trigonometric 
formulae, 4 and zg in Fig. 38 are expressed by : i | 


tan A.4= cos eae . 
(53) 
Mx wel 


=—Cos 
tan AB cosh ae 


Then, 44 and 4g are obtained from 
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Table VIL Values of A? for various cases of & and a/b. 
R=5 k=10 

PP ee 3 3 
0° 0° 0° 0° 0° 0° 0° 0° 
10 287! se9. 9dr SAGA 10 os 8.55) R145. 98.0 
20 = 31:8 9 182 a7 20 +14 Gia nine 
30 1d Sion 2B ne POA 30 AIS em IGA ye 
40 i f49> “219%, 270 40 197, 262k St 
50 137 '-1ag. [ood 50 -161 -220 —272 
60 te Is te “Siee 60 = ese ee 
70 ES pat 1 eee 70 ~~ OTe e120 ee 
80 Bree Cs TE 80 £7 eee tS cares 

90 0 0 0 90 0 0 0 
100 41 53 6.4 100 4.7 61 73 
110 7.9 10.5 12.6 110 91 12.0 144 
120 11.2 15.1 18.4 120 13.0 17.4 21.6 
130 13.7 18.8 23.4 130 16.1 22.0 272 
140 149 21.2 27.0 140 19.7 25.2 32.0 
150 145 213 28.4 150 175 26.1 34.8 
160 11.8 18.2 25.7 160 14.6 23.2 33.3 
170 6.7 10.9 16.4 170 8.5 14.5 23.0 

180 0 0 0 180 0 0 0 

h=30 k=100 

roe 2 3 5 Pas 2 3 5 
0° 0° 0° 0° 0° 0° 0° 0° 
10 “100 184, 1298 10 -106 -189 —323 
20 -169 —282 —395 20 ~178° $290" a0 
30 —20.0 305  —Se6 30 =205~ Lan ans 
40 200, “—288 . 38.7 40 —208 -293 —370 
50 179m» +247 208 50 -185 —-251 —308 
60 -144 -193 23.0 60 ~148 6-196... =35 
70 =100; - =192) . 4155 70 -103 =—194 ) 168 
80 ey Re Sy el 80 Be eltepioy) cegig 

9) 0 0 0 90 0 0 0 
100 51 6.7 7.8 100 53 63 8.0 
110 10.0 13.2 155 110 10.3 13.4 16.0 
120 14.4 19.3 23.0 120 148 19.6 23.6 
130 17.9 24.7 29.8 130 18.5 25.1 30.8 
140 20.0 28.8 357 140 20.8 29.3 37.0 
150 20.0 30.5 39.6 150 20.9 31.3 41.5 
160 16.9 28.2 39.5 160 178 29.0 42.0 
170 10.0 18.4 29.5 170 10.6 18.9 32.3 

180 0 0 0 180 0 0 0 
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AMa= 5 dat Pp, 


2 (54) 
ArNp=Azpgt+p— oe 


In the relation (53), it is obvious that 7X 4=7—)zp,, so that 4A4=—AAg=AAX. The 
values of 4X4 or 4Xg must vanish when the small conductive area is situated 
at p=0, a , 7 and = a, or the conductivity along the auroral zone is uni- 
formly distributed, as already shown in Fig. 36. 

In the practical geomagnetic problem, the line APB divides the lower 
latitudes into two part, in which the horizontal component of the dis- 
turbance field due to the overhead current flow is positive and negative res- 
pectively. In the case of geomagnetic bays, the above two parts are called 
as positive and negative bay areas. The values of 4A are calculated by means 
of (42), (52) and (53) for various cases of k and a/b. In Table VII, the result 
is summarized. 4X amounts to even about 40°. It is indicated here that the 
boundary of the areas of positive and negative horizontal disturbing force 
near the equator shifts about 40° at most according to the position of the 
centre of the highly conductive area along the auroral zone. In other words, 
the mode of current flow in the equatorial region has an uncertainty of its 
phase about 3 hours in local time, even if the mode of air motion in the polar 
region is strictly subjected to local time. 
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Chapter VI. Relation between Sp-Field and Dst-Field 
in the Middle and Low Latitudes 


$1. Development of Sp-field at the time of magnetic storms 


It has been generally accepted that, on magnetically disturbed days, there 
appears disturbance-daily variation, Sp. On the other hand, there exists also 
Dst-field of rather systematic form, after the commencement of magnetic storms. 
In general the superposition of Sp- and Ds-fields with some irregular varia- 
tions of the magnetic field is observed on the earth’s surface. Ds;-field is 
interpreted mainly by the equatorial current ring, the position of which being 
supposed to be several earth radii distant in the earth’s magnetic equatorial 
plane, while Sp-field can be considered to be originated from the current- 
system produced in the upper atmosphere, some 100km above the earth’s 
surface. 

The development of Ds;-field after the sudden commencement of magnetic 
storms in the middle and low latitudes has been fairly well clarified for 
various grades of magnetic storms from a statistical standpoint. On the other 
hand, the development of Ds;-field in the polar region, however, has not yet 
been dealt with, because of of its small magnitude compared with Sp-component 
or other irregular variations in high latitudes. If Ds:-component of the distur- 
bance field be attributed to the equatorial current-ring, as is ordinarily accepted, 
the examination of the development of Ds:-component of the disturbance field 
in the polar region is unnecessary. 

As for Sp-field, on the other hand, it has been detected from a statistical 
study by means of a number of magnetic storms that the intensity of Sp-field 
is most large in the first day of magnetic storms, and it decreases in the 
succeeding days. However, the development of Sp-field with storm-time has 
been never clarified. The writer is told that S. Chapman detected recently 
that Sp-field in the lower latitudes begins to develop within a few hours after 
the sudden commencement of magnetic storms, and attains its maximum 
magnitude about half a day in storm-time in a statistical examination of 
moderate magnetic storms [55]. 

Since the geomagnetic disturbance is considered to b2 caused by the solar 
corpuscle at any rate, there should be expected some mutual relation between 
the development of Sp- and Ds;-fields, though the mechanism for them is differ- 
ent with each other. T. Nagata showed an example of the comparison of the 
development of Sp-field and of Ds;-field with reference to the magnetic dis- 
turbance on May 1, 1933 without sudden commencement [27]. He remarked 
that the development of Sp-field took place in the early stage of the distur- 
bance. In his paper, southward shifting of the auroral zone current in the 
northern hemisphere with the development of Ds-field was also illustrated. 
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The above phenomenon was interpreted by himself by the effect of the magetic 
field of the hypothetical current-ring upon the motion of charged particies 
impinging upon the auroral zone ionosphere [27, 56, 57]. The above results 
must be appreciated for the study of the development of the disturbance field, 
because a mutual relation between Sp- and Dg-fields is presented. 

The examination of the mutual relation between the development of Sp- 
and Ds;-fie'ds is one of the most important problems for clarifying the general 
feature of magnetic storms. The study of this problem by means of a number 


Table I. Mean deviation of the geomagnetic field during 16h 28m—17h 00m 
G.M.T. from the value just before the sudden 
commencement, on April 30, 1933. 


Observatory Abbreviation Ad HAD AZ 
Thule , TH —Ef5Y —207 — 157 
Godhavn GO —35 20 —45 
Scoresby Sund ss 40 —20 —55 
Angmagssalik AN 50 10 20) 
Sveagruvan SV —10 * =—30 = 5 
Bjorndja BJ —40 —30 —50 
Fort Rae FR —25 15 10 
Tromso iis 50 —35 aa 
Petsamo PE 50 — 20 —20 
Matotchkin Shar MS 0) 25) —85 
College-Fairbanks CO —25 —10 5 
Sodankyla SO 50 0 -— 5 
Lerwick LE iS —10 — 3 
Dombas  _ DO 30 5 5 
Meanook ME —15 15 
Lovo LO 40 -15 5 
Sloutsk SL 25 —20 
Agincourt AG 5 20 0 
Abinger AB 50 —15 10 
Val Joyeux vJ 35 —10 5 
Yakutsk YA 70 0 15 
Cheltenham CH 10 13 = 8) 
San Miguel SM 20 ot) 3 
Tucson AUG) 20 8 —5 
Toyohara TY 30 10 2 
San Juan SJ 13 0 ee 
Kakioka KA 28 6 15 
Honolulu HO 15 5 7 
Antipolo AL 20 0 3 
Huancayo HU 50 8 0 
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of polar region data will be desired. In the followings, some remarks concern- 
ing the development of Sp-field in the initial phase of magnetic storms are 


summarized. 


§2. Some examples of Sp-field in the initial phase of magnetic storms 


2.1. Magnetic storm on April 30, 1933 

The sudden commencement of the magnetic storm on April 30, 1933, occur- 
red at 16h 28m G.M.T. In Table I, the mean deviations of geomagnetic field 
during 16h 28 m—17h 00m from the value just before the sudden commence- 
ment are summarized, where the magnetic north, east and downward com- 
ponents are denoted by 4H, H)4D and AZ respectively. Equivalent overhead 
current-system corresponding to the distribution of above-mentioned geomag- 
netic variations in the northern hemisphere is given in Fig. 1. The eastward 
current flow in the lower lati- 
tudes will be attributed to the 
initial increase of the horizontal 
intensity with sudden commence- 
ment. In the high latitudes, 
near the so-called auroral zone, a 
fairly distinct westward current 
seems to flow in the dark hemi- 
sphere, and eastward current in 
the sunlit hemisphere. It looks 
like as if Sp-field of weak intensi- 
ty exists in the polar region at 
that time. Moreover, the current- 
system for the disturbance field 
at 17m 10m, which was already 


given by Fig. 4 in Chapter V, Fig. 1. Equivalent current-system corres: 
ponding to the additional variation during 
16h 28 m—17h 00m G.M.T. an April 30, 1933. 


reveals that a distinct develop- 
ment of so-called Sp current- 
system has been performed by that time. According to the analysis of the 
Same magnetic storm by E.H. Vestine [28] with hourly mean values, the 
development of Sp current-system soon after the sudden commencement is 
also noticed. 
2.2. Magnetic storm on August 3, 1949 

In the next place, the magnetic storm on August 3, 1949, will be discussed. 
Although we have not sufficient polar data at hand for the present and follow- 
ing examples of magnetic storms as in the above case, a rough view of the 
disturbance field will be possible. The storm on August 3, 1949, was of a 
moderately severe one. Geomagnetic variation and ionospheric variation at 
that time were examined by the present writer [58]. 
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The sudden commencement of this magnetic storm occurred at 02h 05m 
G.M.T. on that day. In Fig. 2, records of the initial part of the magnetic 


Cheltenham 


eo ee SE Rt SEE ES ee 
Oh ed Ua Re GC ine Oh i OT Ssh 
August 3, 1949 August 3, 1949 


Fig. 2. Records of horizontal intensity and declination at 
Sodankyla, Dombas, College, Sitka and Cheltenham. 


Table II. List of magnetic observatories. 


Abbreviation Station La S ae Be eel aa - 
So Sodankyla 67.4°N 26.6°E 63.8°N 120.0°E 
Do Dombas 62.1 9.1 62.3 100.0 
SF San Fernando 36.5 353.8 41.0 71.3 
Co College €4.9 2122 64.5 255.4 
Si Sitka 57.0 224.7 60.0 275.4 
Ch Cheltenham 38.7 283.2 50.1 350.5 
Tu Tucson 32.2 249.2 40.4 312.2 
SJ San Juan 18.4 293.9 29.9 Bie 
Ho Honolulu 21.3 201.9 21.1 266.5 
Ko Kodaikanal 10.2 77.5 0.6 147.1 
Ka Kakioka 36.2 140.2 26.0 206.0 
As Aso 32.9 131.0 22.0 198.0 
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storm registered at Sodankyla, Dombas, College, Sitka and Cheltenham are 
reproduced. The list of observatories, whose data were available in the present 
study, is given in Table IJ. As will be seen in Fig. 2, the geomagnetic field 
showed a minor perturbation before the sudden commencement in the high 
latitudes for about an hour, and then, it was assumed that the geomagnetic 
field at 0h G.M.T. be perfectly in the quiet condition. Every 15 minute mean 
values of the disturbance field for each observatory are computed, where 
the disturbance field is defined as the deviation of the magnetic field from the 
value at 0h G.M.T. In the computation of the disturbance field, the magnetic 
variation due to Sg is approximately eliminated previously from the observed 
value. The equivalent current-system for the disturbance geomagnetic field 
at 02h G.M.T., namely 5 minutes before the sudden commencement, is given 
in Fig. 3, where the current-system corresponding to Sg is, needless to say, 


Fig. 3. Current-system at 02h G. M. T. on August 3, 1949. 
Electric current of 10,000 amp. flows between successive full 
stream lines in the direction indicated by arrows. 


eliminated. This current-system just looks like that for Sp or bay disturbance 
and it will be considered that a disturbance of geomagnetic field like a ioe 
disturbance took place just before the sudden commencement. 

In Fig. 4 (a)-(f), the overhead current arrow maps referred to geomagnetic 
coordinates at 02h 15m, 02h 30m, 02h 45 m, 03h 00m, 04h 00m and 05h 00m 
are illustrated. It will be seen in a series of figures that the intense current flow 
in the auroral zone, eastward in American side and westward in European 
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(é) v o4hogm os ke sicko oshogm 
Fig. 4. Progressive change in the overhead current arrows during 
2h-5h G.M.T. on Aug. 3, 1949. 
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side, just in the same mode as that of Sp, appeared soon after the sudden 
commencement, and the current became more intense with time. In the middle 
and low latitudes, on the other hand, the eastward current flow corre- 
sponding to the initial increase of the horizontal intensity remained until 03 h. 
Since then, the current arrows distributed in such a manner as can be identified 
with a typical Sp current-system, as is illustrated in Fig. 4 (e) and (f). The 
additional current-system during 02h 00m—02h 30m, i.e. the current-system 
corresponding to the difference in the geomagnetic field at 02h 30m and at 
02h 00m, is given in Fig. 5. One may say that this additional current-system 
is composed of Sp current-system and the eastward current flow corresponding 
to the initial increase in the horizontal intensity in the initial phase. The 
current-system for the disturbance field at 02h 30m is obtained by the super- 
position of those in Figs. 3 and 5, and shown by Fig. 6. It will be no doubt 
that Sp-field began to develop soon after the sudden commencement in the 
case of this magnetic storm. 


180° 
To Sun 
‘ 


270° 


1008 


Fig. 5. The additional current-system Fig. 6. Current-system at 02h 30m. 
for 02h00m-02h 30m. (Electic current of | (Electric current of 10,000 amp. flows be- 
10,000 amp. flows between successive stream tween Successive stream lines.) 
lines in the direction indicated.) 


2.3. Magnetic storms on January 24, 1949 and July 26, 1946 

There are some cases, in which a marked disturbance of fairly large 
amplitude in the earth’s magnetic field occurs immediately after the sudden 
commencement of magnetic storms. These marked variations in the initial 
phase were studied with respect to two severe magnetic storms, those occurred 
on January 24, 1949, and July 26, 1949 [597]. 

The sudden commencement of the magnetic storm on January 24, 1949, 
took place at 18h 27m G.M.T. The magnetogram recorded at Kakioka (36° 
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Fig. 7. Magnetogram at Kakioka on January 24-25, 1949. 
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Fig. 8. Copies of H- and D-components during the initial . 
phase of the magnetic storm on Jan. 24, 1949, observed at 
Sitka, Cheltenham, Tucson, San Juan, Honolulu, Dombas and 
Hermanus. 
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To the Sun — 


Fig. 9. The current arrow map in 
geomagnetic coordinates at 19h G.M.T. on 
Jan. 24, 1949. (Broken lines in the figure 
show the stream lines of the idealized cur- 
rent-system for Sp-field.) 


To the Sun 


Fig. 10. The current arrow map in 
geomagnetic coordinates at 19h G.MT. on 
July 26, 1946. (Broken lines in the figure 
show the stream lines of the idealized cur- 
rent-system for Sp-field.) 


Table III. Disturbing force of the geomagnetic field 
at 19h G.M.T. on January 24, 1949, 


Observatory Abbreviation pr ae Gos: AH AH, AD AZ 
Sitka SI 60.0° 275.4° — 8107 OY —390Y 
Cheltenham CH 50.1 350.5 —98 63 -8 
Tucson TU 40.4 3122 —162 €0 —27 
San Juan Si) 29.9 ee —34 Bs —2 
Honolulu HO Palak 266.5 8 —32 —15 
Dombas DO 62.3 100.0 67 185 -9 
Kakioka KA 26.0 206.0 120 —38 82 
Hermanus HR —32.7 79.9 57 —34. 8 


Table IV. Disturbing force of the geomagnetic field 
at 19h G.M.T. on July 26, 1946. 


Observatory Abbreviation Lat. ee one, AH H,AD AZ 
Sitka Sil 60.0° 275.4° —416Y —2137 oY 
Cheltenham CH 50.1 350.5 —116 239 23 
Tucson Ui 40.4 312.2 —191 120 —16 
San Juan Sy) 29.9 3.2 11 40. —11 
Honolulu HO 21: 266.5 -—17 0 -15 
Kakioka KA 26.0 206.0 101 8 64 
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14’'N, 140°11’/E) is given in Fig. 7. The magnitude of the sudden commence- 
ment was fairly large there, amounting to 52y in the horizontal intensity. The 
marked increase in the horizontal intensity of about 1.5 hours in its dur atian, 
which looks like a typical bay disturbance, followed immediately after the 
sudden commencement. The magnetic field did not show a large activity 
during following 3 hours. The copies of magnetic record during the initial 
phase of the magnetic storm at Sitka, Cheltenham, Tucson, San Juan, Honolulu, 
Dombas and Hermanus Observatories are given in Fig. 8. It will be noticed 
in Fig. 8 that the remarkable geomagnetic variation followed immediately 
after the sudden commencement everywhere over the world. The disturbing 
force and the current arrow map corresponding to the disturbance field at 19h 
G.M.T., namely 33 minutes after the sudden commencement, is shown in Table 
III and Fig. 9 respectively, where the magnetic field is assumed to be in a 
quiet condition just before the sudden commencement. Although there is no 
observational result in the north of the auroral zone, it will be said that the 
current-system for Sp-field or bay disturbance has developed by that time. 
The unnatural increase in the horizontal intensity observed at Kakioka can be 
interpreted as the development of the intense Sp-field or bay disturbance, 
which just succeeded the sudden commencement of the magnetic storm. 

Another example of the variation in the geomagnetic field during the 
initial phase of a large magnetic storm is found in the storm on July 26, 1946, 
with its sudden commencement at 18h 47m G.M.T. The magnetogram at 
Kakioka shows also an unusual increase in the horizontal intensity succeeding 
the sudden commencement of 114y in the magnitude of the horizontal intensity. 
The disturbing force and the current arrow map corresponding to the distur- 
bance field at 19h G.M.T., namely 13 minutes after the sudden commencement, 
are given in Table IV and Fig. 10 respectively, which are derived in the same 
way as in the former case. The map in Fig. 10 shows quite similar distribution 
of disturbing force as that of Fig. 9, though there are no data of European 
stations in the latter case. 

It will be concluded here that there are some examples of severe magnetic 
storms, in which the development of intense Sp-field can be clearly detected 
immediately after the sudden commencement. The above two magnetic storms, 
however, occurred nearly at the same time of a day. It will ‘be desirable to 
find another examples of storms distributing at various times of the day with 
respect to their sudden commencement. 


§3. Relation between Sp-field and Ds:-field in the initial phase of magnetic 
storms 


In the preceding paragraph, four individual examples of magnetic storms 
were shown, in which a fairly distinct development of Sp-field immediately after 
the sudden commencement of magnetic storms was noticed. On the other hand, 
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the causation of Sp-field is shown to be attributed to the anomalous increase 
in the electrical conductivity of the upper atmosphere over the auroral zone, 
so far. as the dynamo action there is concerned. The anomalous increase in 
the electrical conductivity will be produced by the precipitation of solar charged 
particles from outside the earth along lines of the earth’s magnetic force into 
the auroral zone ionosphere. Recent auroral investigation yields an evidence 
that the solar particles are mainly coraposed of protons, and their approaching 
velocity at the level of auroral display will be a few thousand km./sec. [17]. 
Therefore, it is pointed out that there are some typical examples of magnetic 
storms, in which solar particles are impinging upon the upper atmosphere over 
the auroral zone immediately after the sudden commencement. 

As for the mechanism of the sudden commencement of magnetic storms, 
S. Chapman and V.C.A. Ferraro proposed a theory in 1931, in which the initial 
increase of the horizontal intensity takes place over the world when the front 
of neutral ionized solar stream is in a distance of several earth radii in the 
earth’s magnetic equatorial plane [19]. No alternative theory from different 
standpoint of view has been proposed since then, so far as the sudden com- 
mencement of magnetic storms is considered to be caused by the corpuscles 
coming towards the earth from the sun. If the theory of the sudden com- 
mencement be accepted, the above results concerning the development of 
Sp-field in the initial phase suggest that there are cases where the arrival of 
the corpuscular stream in the auroral zone ionosphere takes place immediately 
after the sudden commencement of magnetic storms. 

It is not yet clear whether the development of Sp-field in high latitudes 
takes place soon after the sudden commencement of every magnetic storm. 
Magnetograms obtained at high latitude stations show some irregular variations 
succeeding the sudden commencement, in almost every magnetic storm. How- 
ever, the irregular variations are not systematically examined whether they 
are of the form of Sp-field or not, for a number of examples. Recent statis- 
tical examination by T. Nagata and H. Ono [60] shows a fairly distinct 
development of Sp-field within several hours in storm-time near the auroral 
zone in the case of large magnetic storms. In their analysis, about 40 severe 
magnetic storms with sudden commencements recorded at Sitka (geomagnetic 
latitude 60.0°N) during 1946-49 were dealt with, and a comparison of the pro- 
gressive development of Sp- and Ds-fields was also given, where the change 
in the horizontal intensity of geomagnetic field observed at Honolulu (geomag- 
netic latitude 21.1°N) was adopted as a mesure of Dy-field. The result shows 
that the magnitude of Sp-field attains its maximum simultaneously with the 
maximum depression of the horizontal intensity in lower latitudes, or the 
maximum of De-field. Their result also indicates the arrival of corpuscular 


stream in the auroral zone ionosphere in the early stage of severe magnetic 
storms, ‘ 
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In the analysis of SC* by T. Nagata [54], some observed facts are illus- 
trated, which suggest that some part of the corpuscular stream, which is 
responsible for the whole process of magnetic storms, are used to arrive at or, 
approach to the upper atmosphere in the polar region of the earth, even 
already at the time of the sudden commencement of magnetic storms. At 
any rate, the corpuscular stream can be considered to arrive at the upper 
atmosphere over the auroral zone almost simultaneously with the sudden 
commencement at the time of magnetic storms. 


§4. Geomagnetic bays and the equatorial current-ring 


Geomagnetic bays are considered as a temporary appearance of Sp-field 
within a few hours. It is often seen that bays occur in the later stage of 
magnetic storms, where the vestige of the equatorial current-ring still remains. 
Nevertheless, typical bays take place frequently on quiet days, at least appar- 
ently. The problem dealt with here concerns whether or not the geomagnetic 
condition before or after the occurrence of bays is really quiet, being free 
from Ds:-field or the presence of the equatorial current-ring. 

52 typical examples of positive geomagnetic bays observed at Kakioka 
(36°14/N, 140°11/E) at night times alone on quiet days in 1950 and 1951 were 
selected, where there was no preceding nor succeeding disturbance at least 
more than several hours before or after respective bays, because the geomag- 
netic field at night on quiet days is free from the diurnal change. The mean 
value of the horizontal intensity of geomagnetic field during three hours before 
the beginning of each bay assumes to represent the geomagnetic condition at 
the time of occurrence of the bay and be denoted by Hg. On the other hand, 
the standard value of the horizontal intensity on quiet condition (H ) in each 
month was obtained through the following procedure. In Fig. 11, the monthly 
averages of the daily mean of the horizontal intensity on five international 
quiet days are plotted for every month in two years. General increase of the 
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Fig. 11. Monthly averages of the daily mean of the 
horizontal intensity on international quiet days. 
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monthly average value with time will contain the secular variation of the 
earth’s magnetic field. Assuming that the secular change show a linear in- 
crease, the train of the above monthly averages can be expressed by means 


of the least square method as 

H=29994.2+2.2ny, 
where the equation gives the standard value of the horizontal intensity in the 
quiet condition in the m-th month from the beginning of 1950. The deviation 
of Hz from H, H 3-H, may represent the magnitude of Ds-component at the 
baytime. In Table V, the selected bays are summarized with their occurrence 
time, values of Hz, H and H 3-H, as well as their maximum magnitude. The 
magnitude of north component of the disturbing force of bays is denoted by 
4H, while the maximum of the horizontal disturbing force by 4F, which is 
obtained by the combination of the changes in horizontal intensity and declina- 
tion. In Fig. 12, the relation between Hg-H and 4F is shown, while the 
‘frequency distribution of Hg-H for the 52 cases is illustrated in Fig. 13. In 
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Fig. 12. Relation between Hg-H and AF. Fig. 13, Frequency distribution of Hp-H. 


these figures and Table V, it is clearly seen that the values of H;-H are 
generally negative, being mostly 0 to -40y. It is also noticed that the larger 
the disturbing force of bay is, the larger is the negative value of 1 3-H, in 
general. Bearing in mind that the night values of H at Kakioka is generally 
larger than its daytime value on quiet days, H should be rather an underesti- 
mation as the reference value for Hg, and the true values of Hp-H will have 
larger negative values than those given here. 
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Table V. Selected typical bay disturbances in 1950-51. 

Beginning Ending AH AF Hp H Hp3-H 

TGs Max, 4 13 a og 12 1g Uggs lage eee 
ie OY fe ge 23 23 29996 30901 — 5 

Apr. ae 2.0 (Os 47 48 29950 30003 —53 
1 A) Te ING 17 18 29974 30003 —29 

Tome 29 205 02 14 28 29996 30003 -—77 

30 =60.0 30° 20 32 33 29964 30003 —37 

May — 17 22.2 I Ass 13 14 29999 30005 — 6 
Sept. 10 23.6 is VAG 35 39 29999 30014 —15 
14 40 Ae 55 19 20 30014 30014 0 

17, -23.8 18 2.0 31 33 30007 30014 —7 

ifsh) PAG? 18 23.4 22 25 29989 30014 —25 

25) 23.5 PROS ws leg 28 32 30003 30014 —Il1 

25 205 28 214 10 11 30013 30014 -1 

Oct. 8920:2 8 21.8 13 16 29993 30017 —24 
18 0.0 Sele i5 23 29989 30017 —28 

Nov. Za 22.5 Tee 0:2 32 32 29989 30019 —30 
14 22.7 15 04 13 13 30012 30019 —7 

23 «(0.8 Hey) VHS) 34 35 30017 30019 — 2 

75) = PSY 20 47 26 ol 30026 30019 i 

Dec. ae. 2 Th this 17 Ze 29988 30021 —33 
6 19.3 6 21.0 28 52 29962 30021 —59 

10 23.6 O25) 16 LT, 30014 30021 — 7 

i) Tepes alsy  7syis) 49 49 _ 29994 30021 —27 

1951 Jan. 14 21.9 14 23.5 30 31 30015 30023 — 8 
LS 226 19 0.6 nis) 16 30031 30023 8 

29 20.0 29 213 ie 22 30015 30023 - 8 

Feb. 26 22.8 Zl Osl, 23 27 30017 30025 — § 
Mar. il OUY/ e219 PAD) 27 30001 30028 —27 
2) Ge) Gy (Oies} 16 I¢f 30014. 30028 —14 

4 OF 4 2.0 10 ia 30019 30028 -9 

TORE 16 10 24.0 47 47 30000 30028 —28 

24 20.9 24 22.9 25 26 30005 30028 —23 

Apr. 2-4 i aks} 16 ilv/ 30022 30030 — 8 
May 4 20.7 4 22.5 24 24 30007 30032 —25 
June a) PAG B Deys 9 17 30029 30034 — 5 
iy COLf ) Pall 16 19 30046 30034 12 

10 1.4 10 28 14 14 30039 30034 5 

July Sy) BAL S229) 45 A7 29978 30036 —58 
Sept. 13 23.4 Tah Glee 31 OZ 30023 30041 —18 
14 219 14 23.1 22 26 30016 30041 —25 

Oct. 11 20.4 pk, ANAS 24 34 30013 30043 —30 
15) 0.6 jb) PAS 24 24 30032 30043 —11 

20 19:5 20 20.5 20 29 30011 30043 —32 

Ope AMES) Pp Gash es) 21 25 30030 30043 —-13 
Zae2.o Peay WBE? Le, 19 30023 30943 —20 

Nov. 4 23.7 Lp V5) 39 44 30007 30045 —38 
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un 92.7" om os" 1 5 16 30044 30045" = 1 
24 «0.3 24 2.3 37 38 30038 20045 —7 
Dec. 12 0.4 Zee 40 41 30925 30048 —23 
19 22.1 19 24.0 29 31 30041 30048 == ff 
23 2.2 Dey eM 50 62 30012 30048 —36 
2922.0 30. (0.4 16 22 30010 30048 —38 


Then, we may safely conclude that the horizontal intensity before the 
occurrence of bays is appreciably smaller than the normal value on quiet days. 
This result may suggest that the negative value of H of Ds:-field exists usually 
at the time of occurrence of bays. In other words, it seems likely that the 
presence of the hypothetical current-ring is a necessary condition for the 
appearance of geomagnetic bays. Most part of the corpuscular stream imping- 
ing upon the upper atmosphere along the auroral zone is, in the present stage 
of investigation, considered to come from the equatorial current-ring [21, 61]. 
Thus, we may be allowed to presume that the corpuscular stream, which 
impinges on the auroral zone ionosphere and results in the bay phenomenon, 
comes also from the current-ring, just as in the case of polar magnetic storms. 
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Chapter VII. Ceomagnetic Variations and the Ionosphere 


§ 1. Height of the current-systems for geomagnetic variations in the 
upper atmosphere 


Direct measurement of the height of the current-system for Sq-field by a 
rocket flight gives a result that the electric current will flow at the height of 
90-105 km. above sea level near the magnetic equator [2]. Jonospheric fade 
, out at the time of a solar chromospheric eruption and the simultaneous geo- 
magnetic variation suggest that the current-system for S,-field exists about 
100km. or less in the height in the observation at middle or low latitudes 
(62, 7, 8]. Although there is no evidence for the height of the current-system 
responsible for S,-field in the high latitudes, the height will not so much differ 
from the above results for temperate and tropical zones. The above height of 
the current-system corresponds to the lower part of EF layer or D layer of the 
ionosphere. There are some Statistical relations between S,-field of geomagnetic 
variation and ionospheric condition, which also seem to support the above con- 
clusion on the height of the current-system for S,-field. The current-system 
for lunar diurnal variation in geomagnetic field is considered to be situated 
mainly also at the level of D layer [63, 64]. 

On the other hand, direct measurement of the height of the current-system 
for Sp-field or geomagnntic bays has not yet been succeeded. An estimation 
of the height of the intense auroral zone current during severe magnetic storms 
or bays, which was carried out by A.G. McNish [65] and T. Nagata [44], in- 
dicates that the height will be about 100km., not larger than 150km. There 
is no other direct proof for the height of the current-system for the polar dis- 
turbance field at present, and we know no conclusive result concerning the 
height of the current-system, especially in the middle and low latitudes. 
Summarizing the above-mentioned fairly reliable results, however, it will be 
not so unreasoeable to suppose that the electric current-system for S,-field flows 
mainly around the lower part of E layer or D layer, while the height of Sp 
current-system will be situated around the level of the maximum electron 
density of E layer, a little upper portion than the former. 


§ 2. Electrical conductivity and air motion in the upper atmospere 


Radio sounding of the ionosphere gives mainly the distribution of electron 
density with height, and it detects the presence of £ and F layers. On the 
other hand, the distribution of heavy ions with height cannot be directly ob- 
tained from the usual observation, and it should be estimated from theoretical 
considerations for the electrical equilibrium of the ionized state. It has been 
shown that the contribution of ions for the direct current conductivity of the 
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ionosphere is much predominant than that of electrons from the magneto-ionic 
theory, so far as the number density of ions and of electrons assumes to be 
not so much different with each other. Since the direct current conductivity 
of the ionospheric region is dealt with in the geomagnetic problem, a reliable 
distribution of ions with height is strongly desired to obtain the distribution 
of specific direct current conductivity with height. 

The estimation of the direct current conductivity based upon the present 
knowledges has been tried by many authorities, where it was compelled to 
introduce various assumptions on the equilibrium of ions and electrons, and 
other physical properties of the upper atmosphere. According to a result by 
K. Maeda [47], the distribution of specific direct current conductivity with 
height does not show two maxima corresponding to the levels of current-systems 
for Sg- and Sp-fields presumed in the preceding paragraph, but shows its 
maximum nearly at the levels of maximum electron density of & and F layers. 
The total or integrated electrical conductivity calculated by him is about 
2x10-7e.m.u. in the daytime. On the other hand, an estimated value of the 
total conductivity of effective thickness for Sg current-system is required to be 
about 5x10-8e.m.u. over the middle and low latitudes, after an analysis 
of Sgqg variation by T. Nagata [8]. This value required from geomagnetic 
problem seems to harmonize with that estimated from ionospheric standpoint. 

There have been presented some evidences that the density of heavy ions 
will be far larger than that of fre2 electrons in the lower part of the iono- 
sphere [66, 67]. Recent pressure measurement with rockets from the ground 
up to 160 km. enables to calculate the distribution of density and temperature 
with height up to the upper part of E layer [68]. Rough measurements of the 
effective electron density up to 120km. were carried out by H.E. Newell, Jr. 
[69] previous to the above observation, also by means of a rocket flight. With 
the aid of the above experimental results or further explorations, it will be 
expected to obtain more reliable distribution of electrical conductivity with 
height in the lower part of the ionospheric region in near future, though we 
have yet insufficient knowledges about the physical properties of the ionosphere 
at present owing to the poverty of experimental evidences. 

As for the air motion in the upper atmosphere, various kinds of results 
have been presented in the recent several years. They were found by many 
researchers over the world with different methods, for example, by examining 
the shift of ionic clouds or luminous parts of night sky light, the movement of 
irregularities in ionization of the ionospheric region, or by some special ex- . 
periments with radio waves, and so on [70]. Although the detail of the results 
is not mentioned here, it was noticed that the lateral motion in the so called 
ionospheric region has a velocity of the order of 100 m/sec. with some turbulent 
motion, so far as the shift of ionic clouds or parts of specially ionized state 
assumes to be carried by the general flow of the air mass. It is not yet pos- 
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sible to compose a consistent world-wide distribution of the air motion through- 
out the ionospheric region from the present’ accumulation of experimental or 
analytical results. 

On the other hand, tidal oscillation of the upper atmosphere has been 
studied by examining the periodical oscillation of the height of ionospheric 
layers. D.F. Martyn [71] reported the review of tidal phenomena in the iono- 
sphere, and proposed an opinion concerning the formation of the current- 
systems for geomagnetic variations. There have been proposed many other 
models of ionospheric circulation, and the relations between circulatory motion 
in the upper atmosphere and geomagnetic variations have been also discussed 
by S. Matsushita [72] and others. S. Chapman [73] also reviewed the present 
knowledge on the phenomenon of atmospheric tides and oscillations from the 
standpoint of atmospheric dynamics in detail, and emphasizes the necessity of 
further theoretical and observational investigations. Then, it is a little pre- 
mature to conclude about the world-wide circulatory motion in the upper 
atmosphere by means of the present insufficient and non-systematic data. 

The air motion in the electrically conductive ionospheric region, which is 
presumed in the ordinary dynamo theory for S,-field, has no established 
evidence from the theoretical or observational standpoint. However, owing to 
the utmost adequacy of the dynamo action for the interpretation of various 
kinds of geomagnetic variations, the hypothetical air motion is allowed to 
presume, without any serious objection. In Chapter III, §2, the writer pro- 
posed a double layer hypothesis for an interpretation of S,- and Sp-fields from 
the standpoint of the dynamo theory. In the hypothesis, the presence of two 
current sheets concentric with the earth is assumed. The hypothesis seems to 
be not so unreasonable, at least apparently. But there is no positive proof to 
approve of the air motion in the hypothesis, namely the movement in the two 
layers is nearly in opposite direction with each other. Further research on the 
conclusive distribution of the air motion with height in the upper atmosphere 
will give a criticism on the hypothesis. 


§ 3. Ionospheric perturbation in the middle and low latitudes accompany- 
ing geomagnetic disturbance 


Since routine observation of hourly ionospheric aspects has been started 
about twenty years ago, there have been best chances to experience typical 
ionospheric changes accompanying magnetic storms of severest class since the 
beginning of the continuous magnetic recording, say at the time of magnetic 
storms on April 16, 1938, March 24, 1940, and others. Jonospheric change at 
those times observed in middle and low latitudes was fairly well described in 
detail [74, 75]. Although various complicated features of the ionospheric dis- 
turbance at the time of magnetic disturbance have been known at present, 
only an outline of them will be summarized in the followings. 
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It has been generally recognized that F2 layer is more sensitive for geo- 
magnetic disturbance than E layer. Though the phenomenon of ionospheric 
perturbation accompanying geomagnetic disturbance has been noticed for a 
long time, systematic examination of its cHaracteristics has recently started, 
for example, by E.V. Appleton and W.R. Piggott [76], H. Uyeda and Y. Arima 
[77], D.F. Martyn [78] and others. #2 layer disturbance is dealt with in the 
first place in every work, because of its pronounced change with geomagnetic 
activity compared with the disturbance in the lower ionospheric region, and of 
its perpetual finding by the usual radio sounding throughout day and night. 

F2 layer is apt to be disturbed especially in the daytime of summer 
season. Most characteristic feature of F2 layer disturbance in the summer 
season is the diminution of maximum electron density and the upheaval of the 
layer. The change in the daily mean value of maximum electron density of 
F2 layer, however, shows rather systematic change with season in the 
middle latitudes, namely it decreases in summer, while it increases on the 
contrary in winter with geomagnetic activity [79]. The thickness of F2 layer 
also changes at the time of magnetic storms. T. Obayashi [80], K. Sinno [81] 
and D.F. Martyn [78] tried to separate the progress of F2 layer disturbance 
after the , sudden commencement of magnetic storms into two parts, namely 
the variation depending upon universal time and that depending on local time, 
just similarly as the separation of the geomagnetic disturbance field into Ds¢ 
and Sp components, by means of Japanese and world-wide observations. H. 
Uyeda and Y. Arima [77] classified F2 layer disturbance over the world into 
some types and clarified the appearance tendency of each types with season. 
H. Kamiyama [82], or K. Miya and N. Wakai [83] examined individual exam- 
ples of the development of F2 layer disturbance over the world. According 
to all of the above results, F 2 layer disturbance in the middle and low latitudes 
does not suddenly start from the time of commencement of magnetic storms in 
general, but it becomes remarkable some hours later. On account of the 
above tendency, the radio warning is ready in time if it is issued after the 
sudden commencement, so far as the wireless communication circuits through 
polar region are not dealt with [84]. Since the recent detailed analysis of 
the world-wide data is, regrettable to say, lacking in the polar region observa- 
tion, only the distribution of F2 layer storm in the middle and low latitudes 
can be looked over. Notwithstanding, K. Miya and N. Wakai [83] showed 
an example of the propagation of F2 layer disturbance from north to south 
in the northern hemisphere, which originated almost simultaneously with the 
sudden commencement of the magnetic storm in the local midnight region over 
the auroral zone. The propagation characteristics must be further examined 
by means of a number of examples, because statistical evidence is not yet 
presented. Various analyses by other researchers over the world concerning 
F2 layer disturbance in the middle and low latitude stations show a similar 
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tendency as mentioned above. 

As for the disturbance in the lower part of the ionospheric region, the 
occasional production of intense sporadic E layer and the intense absorption of 
radio waves in the lower part of the ionosphere have been remarked. Sporadic 
E Jayer, which is originated nearly at the level of normal E layer, often ob- 
scures F layer by reflecting the energy of electromagnetic waves transmitted 
into the layer. It sometimes happens that no echo from the ionospheric region 
can be obtained, when the absorption of radio waves in the lower part of E 
layer becomes extremely intense. Sporadic E reflection is considered now to 
be due to some cloudy or non-uniform distribution of electron density near the 
level of F layer, and there are some evidences that sporadic ionizations are 
localized over a relatively small area. In moderate latitudes, the sporadic 
ienization is found to be much stronger in summer than in winter with a 
maximum near midday. Systematic study for sporadic EF ionozation over the 
world is now going on. Although almost continuous observation of the iono- 
sphere has been succeeded and some characteristics of sporadic ionization have 
been clarified, the process of perturbation in the lower part of the ionospheric 
region accompanying magnetic disturbance is not yet fully understood. 


§ 4. Ionospheric perturbation in the high latitudes accompanying geo- 
magnetic disturbance 


Since there is a definitely established evidence that solar particles attack 
the upper atmosphere over the auroral zone, it is, needless to say, an important 
problem to clarify the ionospheric perturbation in the polar regions at the 
time of the geomagnetic disturbance. Here, an outline of the characteristic 
features of the ionospheric region in high latitudes will be reviewed in the 
following. 

The distribution of critical frequency or maximum electron density of F2 
layer over the world does not show remarkable anomaly in the polar region 
including the auroral zone [85], though it reveals geomagnetic control in the 
lower latitudes [85, 86, 87]. Recent study of the world-wide distribution of 
the total number of electrons below the level of maximum electron density in 
a column of unit cross-section in the F2 region instead of its critical frequency 
or maximum electron density alone, which was carried out by J.A. Ratcliffe 
[88], shows that the apparent geomagnetic control of the maximum electron 
density of F2 layer seems to be removed when the total electron number is 
‘dealt with, and the total electron number is closely related to the zenith angle 
of the sun’s ray, though not so simply as Chapman's law. As to the character- 
istics of the change in F2 layer accompanying geomagnetic disturbance, the 
decrease in the maximum electron density has been reported [89, 90]. Virtual 
height of the layer seems to increase with the decrease of maximum electron 
density. Recent study by T. Nagata and T. Oguti [91], in which the ionos- 
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pheric change in F2 layer after the sudden commencement of magnetic storms 
observed at College (geo.nagnetic latitude 64.°5N) during 1941-46 were ex- 
amined, clarifies that the maximum electron density decreases in summer and 
rather increases in winter on the whole, just similarly as in the case of 
middle latitudes, and the said change begins almost simultaneously with the 
sudden commencement of magnetic storms. Remembering that the solar 
‘ particles will impinge upon the auroral zone ionosphere almost from the begin- 
ning of magnetic storms, the detected F2 layer disturbance will be directly 
related to the precipitation of solar particles. Furthermore, the annual change 
of the diurnal variation of F2 layer disturbance there indicates that the 
maximum electron density decreases or increases according as the sun’s radia- 
tion beats upon the ionospheric region or not, resulting in the general tendency 
of decrease and increase of the maximum electron density in summer and 
winter seasons. In general, F2 layer disturbance is strongly subjected to the 
solar radiation. Individual examination of the ionospheric change by means of 
polar region data by J.H. Meek [92] shows also the same characteristis of 


disturbance. 
It is well known that the electron density of E and F layers of the iono- 


sphere varies with the sun’s zenith angle, and it is generally explained by a 
modified Chapman’s law. It was pointed out by J.C. W. Scott [93] that 
the sensitivity of E layer to solar angle is very low near the auroral zone, 
though it follows the Chapman’s law in the remaining part of the world, 
even in the polar cap, so far as the monthly mean values are concerned. 
Since the monthly mean values will be affected by the values on magnetically 
disturbed days, the said anomaly of the distribution of E layer over the 
auroral zone alone will be attributed to some geomagnetic effect of E layer 
there, in a statistical meaning. 

Sporadic £ ionization at the time of polar magnetic disturbances is one of 
the most pronounced phenomena. Early study by E.V. Appleton and others 
[94] showed a close relation of the occurrence of sporadic E ionization to the 
polar magnetic disturbance. L.V. Berkner and H. W. Wells [95] also showed 
some evidences for the relationship of the abnormal ionization to geomagnetic 
bays. Intense absorption of radio waves is also an important feature of the 
ionospheric perturbation accompanying geomagnetic activity. The absorbing 
layer seems to be produced in the lower part of E layer. Statistical and in- 
dividual examinations of the polar ionospheric perturbation have been accumu- 
lated, which were compared with not only geomagnetic activity, but also with 
wireless communication data. All of them show extremely close mutual rela- 
tionship among them. 

According to the recent opinion about the ionospheric disturbance in high 
latitudes, sporadic ionization or absorbing layer is produced in a limited region, 
though it sometimes spreads over. considerably wide area. For example, M. 
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Barré and K. Rawer [96] noticed that a high absorption occurring in the iono- 
sphere exists in a limited region from their experiment carried out around 
Adelie-Land (66°50/S, 141°25/E), or H.W. Wells [97] concluded that anomalous 
ionization shows patchy distribution over the auroral zone at the time of 
geomagnetic bays. There have been presented many other evidences for the 
local occurrence of anomalously ionized area. 

As for the relation between the auroral display and the geomagnetic dis- 
turbance, E.H. Vestine convinced the geographic incidence of aurora and 
magnetic disturbance both in the northern and southern hemispheres [98, 991]. 
L. Harang [13] showed remarkable parallelism of the auroral luminosity curve 
and the simultaneous geomagnetic activity. It has been noticed that an 
anomalous condition of wireless communication takes place at the time of 
auroral display. There have been also presented some experimental results 
concerning the reflection of radio waves from the polar aurora, where the 
frequency of transmitted waves is much higher than that used in the ordinary 
ionospheric exploration or usual wireless communication for a long distance 
(100, 101]. Summarizing all the results concerning the relation between aurora 
and geomagnetic disturbance, it may be allowed to consider that the luminous 
part of the aurora corresponds to the highly ionized space and the current- 
system for geomagnetic disturbance field will be originated around there. In 
other words, auroral display, sporadic ionization and production of current- 
system for geomagnetic disturbance field are three phenomena in the upper 
atmosphere caused by one physical origin, precipitation of charged particles 
into the upper atmosphere. 


§ 5. Anomalous ionization in the upper atmosphere over the auroral zone 
caused by the precipitation of charged particles 

It was shown in the preceding paragraph that the anomalous ionization in 
the lower part of the ionospheric region takes place over the auroral zone at 
the time of polar magnetic storms or intense auroral display. The ionization 
would be produced by the precipitation of solar stream there, and the stream 
is proved to contain protons [15, 16, 102]. Some attempts have been proposed 
recently concerning the physical mechanism for auroral display and geomagnetic 
disturbance caused by impinging protons, by L. Vegard [103] and others. M. 
Sugiura, M. Tazima and T, Nagata [104] gave a theoretical result on the ioni- 
zation in the upper atmosphere caused by extra-terrestrial protons. Their 
result will be briefly mentioned about in the following. 

For the distribution of pressure and density in the upper atmosphere, a 
model proposed by N.C. Gerson [105] was adopted, while the composition of 
the atmospheric materials assumed as the same as that near the ground. In 
practice, the following five cases of the upper atmospheric condition were dealt 


with, namely, 
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I Temperature at 300km. assumes 300° K 

I BS 500° Ac for night, 

Il rp 800° K 

V ye 1200° K for day in winter, 
Yy #5 2000° K for day in summer, 


respectively. The result of calculation about the depth of penetration of protons 
with energies up to 1 MeV is illustrated in Fig. 1 for the five cases of pressure 
distribution given above. It will be seen in Fig. 1 that the energy of protons 
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Fig. 1. Depth of penetration of protons. Curves I, II, Ill, IV and V correspond 
to the assumed atmospheric models I-V. 
(after M. Sugiura, M. Tazima and T. Nagata) 


penetrating to the height of 100km. above the earth’s surface is 680 keV, or 
1.1x 109 cm/sec. in the initial velocity. In the case of the practical problem, 
charged particles should impinge with a spiral motion along lines of the earth’s 
magnetic force. The reduction on the penetration range was proved to be the 
order of several km., and it may be practically neglected here. A rough esti- 
mation of the statistical fluctuation in the depth of penetration due to strag- 
gling showed that it amounts to 0.8~0.4 km. for protons stopping at 100~§0 km. 
height. It depends on the model of the atmosphere for protons of lower energy. 
Even in the case of the model V in Fig. 1, it amounts to about 4km. and 8 km. 
for protons penetrating into the level of 140 km. and 220km. respectively. At 
any rate, the fluctuation in the penetration range due to straggling is of the 
Same order of magnitude as the uncertainty in the model atmosphere. 
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The rate of ionization estimated by their calculation is shown in Fig. 2 for 
energies of 1 MeV, 100 and 10 keV for the cases of five models, respectively. 
The rate of ionization has a sharp maximum near the bottom edge of penetra- 
tion, and decreases rapidly with height. If the disappearance of free electrons 
assumes to follow a recombination law, the relaxation time of the ionized state 
is shown to be 10 seconds ~2 minutes. Accordingly, the anomalous concentration 
of ions and electrons will be diluted away in a few minutes after the removal 
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Fig. 2. Rate of ionization by 1 MeV, 100 keV and 10 keV protons. - Curves I, II, 
IIT, IV and V correspond to the assumed atmospheric models I-V. 
(after M. Sugiura, M. Tazima and T. Nagata) 


of ionizing source. When we discuss the geomagnetic variation of its period 
longer than several minutes in conjunction with simultaneous ionospheric vari- 
ation, it will be allowed to consider that there is incoming of ionizing agent 
throughout the polar magnetic disturbance in question is observed, as a first 
approximation. The density of precipitating protons required to maintain the 
ionized state was roughly estimated at 1proton/cm’ or less in their result, 
which can be considered as a reasonable value. 

When the proton stream of 1/cm3 in density and 109cm/sec. in initial velo- 
city or about 500 eV in energy is impinging, about 2x10’ ion pairs can be 
produced in a unit volume near the bottom of the penetration range, because 
one proton produces about 2x10-2 ion pair/cm. as will be seen in Fig. 2. 
The number of ions and electrons does not diminish so far as the proton stream 
continues to penetrate, since the recombination coefficient can be considered to 
have a value of the order of 10-8cm3/sec. The maintenance of the above 
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value of electron density, 2x10’/cm3, seems to be adequate in the interpreta- 
tion of the reflection of high frequency radio waves, which is found experi- 
mentally [100, 101]. The specific electrical conductivity transverse to the 
magnetic field in the upper atmosphere co, is given by 


Vi ry 
o. =e Ear + me iP EOP ) 

where N; and N,, mi and me, vi and ve, wi and we denote number density, 
mass, mean collisional frequency, gyro-frequency of ions and electrons respec- 
tively, and e their electric charge. If it assumes now Ni=Ne=2x107/cm; 
mi=3.3X 10-2 gr.; vj=103/sec., v,=105/sec. (near the level of E layer); wi=3x 
102/sec., we=107/sec. (corresponding to the magnetic field of 0.6/7"), a roughly 
estimated value of specific conductivity amounts to 1.7x10-e.m.u. If the 
thickness of highly ionized layer assumes 20km., the integrated conductivity 
amounts to 3.4x10-5e.m.u. The above value of the integrated electrical 
conductivity is 102~103 times as much as that required from the dynamo theory 
for Sg-field. Then, the anomalous increase in the electrical conductivity in the 
upper atmosphere, which was introduced for the interpretation of the mechanism 
of geomagnetic disturbance field, seems to be quite naturally explained by the 
above process. The penetration depth of the proton stream of 109cm/sec. in 
initial velocity is a little above 100km from the ground, as will be seen in 
Fig. 1. The lower boundary of the intense aurora is about 80~90km. in 
height. For the proton stream producing such an intense auroral display, the 
initial velocity is a little larger than 109%cm/sec., and the ionizing power is 
also a little high than that mentioned above. At any rate, the proton stream 
of about 1/cm or less in its density and about 109cm/sec. in its initial velocity 
causes a sufficient increase in the electron density and electrical conductivity 
in the lower part of the ionospheric region. 
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Chapter VIII. Conclusion 


§1. Constitution of the polar magnetic disturbance field and physical 
mechanism for polar magnetic storms 


Throughout the present study, the world-wide disturbance field of polar 
origin at the time of typical magnetic storms as well as geomagnetic bays 
was examined in detail, especially with respect to its spatial distribution and 
progressive change. The analysis was carried out by means of the world-wide 
magnetogram copies, so as to examine instantaneous aspects of the disturbance 
field comprehensively, which are hardly grasped by the analysis with hourly 
values only. The results of the analysis are as follows. 

The most simple pattern of the polar disturbance field is a “ doublet-field ” 
type one, which is nearly represented by the current distribution on the 
electrically conductive upper atmosphere when an electric doublet is situated 
along the auroral zone. The above elementary disturbance of the doublet-field 
type can be detected in every stage of magnetic storms. ‘There are some 
typical cases when almost perfect form of so-called Sp-field appears. Even 
then, superposing variations of shorter periods upon the smoothed curve of 
geomagnetic variation reveal a doublet-field type pattern in typical cases. 
Periods of these elementary disturbances detected in the present study are 
ranging from a few tens of minutes to about an hour or so. Therefore, these 
elementary disturbances will be generally smoothed out as irregular variations 
in deriving hourly values, and the analysis by means of hourly values cannot 
detect the elementary disturbances, unless typical ones of the duration longer 
than an hour take place over the world. 

The centering position of the elementary disturbance along the auroral 
zone originates at various locality, though it takes place most frequently 
around local midnight. The average current-system of a number of elementary 
disturbances, whose origins are distributed nearly uniformly along the auroral 
zone, reveals the mode just similar as the idealized Sp current-system. Since 
S. Chapman developed his eminent works on the statistical study of the 
distribution of the geomagnetic disturbance field, and derived the idealized Sp 
current-system, the disturbance field of polar magnetic storms in individual 
cases have been explained by the appearance of Sp current-system of appro- 
priate intensity with some deformation in respective cases. The present result 
clarifies that the idealized Sp current-system has rather a statistical meaning, 
namely the Sp current-system is the average of a number of elementary dis- 
turbances, and the most fundamental pattern for the polar disturbance field is 
no more than the doublet-field type pattern. 

. Historically speaking, a doublet-field type pattern of the polar magnetic 
disturbance was noticed in the early research by K. Birkeland. Since then, 
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considerations on such a localized magnetic disturbance have not been dealt 
with systematically, because the most important problem has been to clarify 
the world-wide regular feature of the disturbance field. It is shown in the 
present study that the regular feature of the polar disturbance field is some 
statistical result of the summation of a large number of elementary distur- 
bances. 

An elementary disturbance is explained to be caused by the production of 
a small highly conductive area in the upper atmosphere over the auroral zone 
from the standpoint of the dynamo theory. In the present stage of investiga- 
tion for geomagnetic variations, dynamo action in the upper atmosphere should 
be considered to be the most important cause, because of its applicability for 
all kinds of regular magnetic variations, though the theory must be checked 
by the further study on the ionospheric motion and the distribution of electrical 
conductivity with height in the upper atmosphere. The dynamo theory 
succeeds in interpreting the appearance of perfect and partial Sp-field by the 
perfect and partial increase in the electrical conductivity of the ionospheric 
region along the auroral zone respectively. In this view point, theories pro- 
posed by H. Alfvén [20] and D.F. Martyn [21] must be re-examined. Mathe- 
matical and experimental examinations of the dynamo theory for the polar 
magnetic disturbance field are give in the present study. Mathematical treat- 
ment in this study is, however, only applied for a quasi-stationary state. Since 
the electrical conductivity in the upper atmosphere along the auroral zone will 
vary fairly rapidly with time, the current distribution in the ionospheric region 
should be obtained through the electromagnetic dynamics for a_ transient 
phenomenon, in which the self-induction of the ionosphere must be taken into 
account. The picture for the constitution of polar magnetic storms in the 
present study is based upon the comparison of instantaneous aspects of the 
disturbance field with calculated current-systems for stationary states. Never- 
theless, it seems that no essential alteration will be required in the above 
conclusions for the constitution of polar magnetic storms. 

On the other hand, the anomalous increase in the electrical conductivity 
in the upper atmosphere over the auroral zone, about 100km. in height, is 
interpreted by the precipitation of protons of several hundred keV in energy, 
which come from outside the earth. Therefore, it happens that the extra- 
terrestrial protons impinge upon the auroral zone ionosphere now here, now 
there, sometimes more than two localized regions simultaneously, or sometimes 
almost uniformly along the auroral zone wi j i i 
the whole course of me magnetic storms en ee ee Seine 

r protons are 
proved to impinge upon the upper atmosphere along the line of the earth’s 
magnetic force, and they are considered to come mainly from the equatorial 
current-ring. The precipitation of charged particles can be seen even in the 
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early stage of magnetic storms, sometimes immediately after or almost simul- 
taneously with the occurrence of the sudden commencement of magnetic 
Storms. Since the life time of the anomalously ionized state is considered to 
be the order of a minute, it may be reasonable to presume that the charged 
particles are impinging whenever the polar disturbance field is observed on 
the earth’s surface, so far as we concern geomagnetic variations of longer 
than several minutes in their period. 


§2. Geomagnetic bays 


Geomagnetic bays are considered to be a temporay appearance of the polar 
magnetic storm within a short duration of time. The most simple distribution 
of the disturbing force of bays over the world reveals also the aspect of an 
elementary polar disturbance of the doublet-field type pattern. In general case, 
the development and decay processes of individual bays are not simply con- 
sidered as the mere intensification and diminution of the disturbance field such 
as the average current-system for bays obtained by many authorities. The 
intensity of the auroral zone current does not vary uniformly along the auroral 
zone with time, but the westward auroral zone current in the dark hemisphere 
reaches its maximum magnitude a little earlier than the eastward one in the 
sunlit hemisphere does, in a general tendency, and this tendency results in the 
so-called looping phenomenon of the disturbing force in the middle latitudes. 

Since the presence of the equatorial current-ring at the time of bays 
seems to be fairly expectable as examined in the present study, the mechanism 
for geomagnetic bays will be quite the same one as that for polar magnetic 
storms, namely charged particles impinge upon the upper atmosphere over the 
auroral zone from the equatorial current-ring along the line of the earth’s mag- 
netic force. Then, the said difference in the variation of the auroral zone current 
intensity in both hemispheres might be caused by some characteristic ejection 
of corpuscles from the equatorial current-ring to reach the upper atmosphere 
in the dark hemisphere earlier than the sunlit part. The fact examined by Y. 
Kato and J. Ossaka [106] that short period pulsations of very small amplitude 
are remarkable only before the disturbing force of bays reaches its maximum 
value on the magnetogram in the middle latitudes, will indicate that the cor- 
puscular stream is more rich in the fluctuation of its intensity before the 
occurrence of maximum intensity compared with the period of diminishing 
intensity. 

The occurrence of geomagnetic bays show a recurrence-tendency that they 
sometimes occur with interval of about 24 hours for some days. The recur- 
rence-tendency should be attributed to the recurrent ejection of charged 
particles from the equatorial current-ring from the above-mentioned standpoint 
for the mechanism of bays. Frequency distribution of occurrence of bays with 
respect to local time over the world shows an apparently pronounced diurnal 
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variation, namely bays occur most frequently around local midnight. This 
fact shows that charged particles are apt to impinge upon the earth’s upper 
atmosphere in the dark hemisphere. These two facts may be also characteristc 
features of ejection of charged particles from the equatorial current-ring. 

The recurrence-tendency of geomagnetic bays reminds us that charged 
particles are expelled from nearly the same position in the equatorial current- 
ring with nearly 24 hour intervals. Charged particles will be jetted out from 
a part of the equatorial current-ring, at which some instability for the circular 
motion of particles is produced. Radial stability of the equatorial current- 
ring has been discussed by V.C.A. Ferraro [107] and others. From the 
characteristics of occurrence of geomagnetic bays just mentioned above, the 
equatorial current-ring should have some characteristic instability, which 
results in the general tendency of geomagnetic bays. The obliquity of the 
earth’s magnetic axis with its rotation axis seems to be an important factor 
for the recurrence-tendency of geomagnetic bays, and a mathematical proof 
is desired. Anyhow, the instability of the equatorial current-ring must be 
examined in detail. 


§3. Further problems on the polar magnetic disturbance 


In the present study, the polar geomagnetic disturbance field is shown to 
be originated in the production of highly electrically conductive area along the 
auroral zone in the upper atmosphere, and the anomalous increase in the 
electrical conductivity is interpreted by the ionizing effect of precipitating 
protons from the equatorial current-ring. The energy of incoming protons is 
required to be several hundred keV in order to penetrate down into about 
100 km: from the earth’s surface. The initial velocity of protons of above- 
mentioned energy is about 109cm/sec., though the solar corpuscular stream has 
been estimated to travel from the sun towards the earth with its mean velocity 
of about 10®cm/sec. Then, charged particles must be accelerated by some 
mechanism before they impinge upon the earth’s upper atmosphere. S. Chapman 
[18] showed a possibility of getting such a high energy for extra-terrestrial 
protons. D.F. Martyn [21] proposed a mechanism of acceleration of charged 
particles by taking the effect of electric polarization field in the equatorial 
current-ring into account. Further examinations are required for clarifyinz 
the mechanism of acceleration of charged particles. 

In the present stage of investigation, charged particles impinging upon the 
auroral zone ionosphere along the line of the earth’s magnetic force are con- 
sidered to come mainly from the equatorial current-ring, except in the initial 
phase of magnetic storms. It is already shown that a part of the solar stream, 
which approaches the earth nearly along the line of magnetic force, is scarcely 
affected by the retarding force due to the induction within the stream. 
Even if the particles impinge upon the earth’s atmosphere over the so-called 
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auroral zone without retardation, their initial velocity should be about 109 
cm/sec. in order to cause Sp current-system in the initial phase when the 
current-ring is not yet formed. The above value of velocity contradicts with 
those hitherto presumed, namely about 108cm/sec. If the velocity of the solar 
stream during the travel from the sun towards the earth assumes about 108 
cm/sec., some mechanism for the acceleration must be considered in the vicinity 
of the earth, in order to interpret the development of Sp-field in the initial 
phase of magnetic storms. 


There are some unsolved problems concerning the nature of solar stream 
as the origin of geomagnetic disturbance. The above-mentioned point is one 
of them. The approaching velocity of the solar stream has been estimated 
from the correlation between solar phenomena and terrestrial magnetic dis- 
turbances. Solar chromospheric eruption was first considered to eject the 
solar stream, but it has been shown that there is not so distinct correlation 
between solar eruptions and magnetic storms in the statistical study. 27-day 
recurrence-tendency reminds us the presence of a long-lived source of the 
stream on the sun’s surface. It will be hardly expected that the chromospheric 
eruption is accompanied by the ejection of the solar stream of such a long 
duration. Since any other pronounced phenomenon on the sun’s chromosphere 
was not noticed to have direct connection with geomagentic disturbance, 
“M-region” has been assumed, from which the solar stream flows out. 
Identification of M-region with some solar phenomenon has been tried by C.W. . 
Allen [108],.M. Waldmeier [109], K.O. Kiepenheuer [110] and others, after 
the detailed examination of solar and terrestrial relationships. Their results, 
however, do not seem to give a final satisfactorily decisive conclusion. Nowa- 
days, coronagraph and solar noise observations are also available for the 
external and internal parts of the chromosphere, and some preliminary reports 
are presented. It is an urgent problem to determine what is the true source 
of the solar stream responsible for the geomagnetic disturbance, and what 
kind of characteristics does the stream show in the extra-terrestrial space. 
A peculiar 27-day recurrence-tendency of magnetic storms suggested by M. Ota 
[111], or occasional appearance of S.C.-type perturbations over the world report- 
ed by W. Jackson [112] and many other problems will be expected to be inter- 
preted by some characteristic feature of the solar stream. 

Since the ionosphere is approximated as an electrically conductive shell 
surrounding the earth, the change of the magnetic field originated outside the 
ionosphere, namely by the equatorial current-ring, should cause the induced 
current in the ionosphere, which shields the change of magnetic flux inside it. 
The shielding effect of the ionosphere when the periodical change of the 
external uniform magnetic field along the earth’s magnetic axis is applied, has 
been estimated. M. Sugiura and T. Nagata [113, 114] calculated the effect 
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for the case of uniform distribution of electrical conductivity in the ionoshere 
fairly in detail, while A.A. Ashour and A.T. Price [115] for a non-uniform. 
case. Both of them showed that the shielding effect is appreciably large, so, 
that the effect must be taken into account in practice, when we deal with the 
short period geomagnetic variation on the earth’s surface in connection with 
the external change in the current intensity of the equatorial current-ring.. 
When severe magnetic storms are going on, and both Sp and Des-fields are. 
present, the world-wide distribution of electrical conductivity of the ionospheric 
region is clearly different from those approximated in the above calculations, 
because the conductivity over the auroral zone will be far larger than that in 
the other regions. The induced current in the ionosphere due to the change 
in Ds-field at that time may show some complicated mode of flow, and may 
result in some special type of geomagnetic variation over the world. This 
effect must be theoretically estimated. 

Very short period pulsations are generally recorded at high latitude 
stations at the time of polar magnetic storms. There are some cases when 
giant pulsations are observed only near the auroral zone. These kinds of 
short period variations are considered to be originated in the rapid change in 
the current flow in the ionospheric region and not outside the region, when 
the shielding effect of the ionosphere is taken into account. From the stand- 
point of the dynamo theory, the above-mentioned short period geomagnetic. 
_variations will be attributed to the rapid change in the electrical conductivity 
of the upper atmosphere or some turbulent motion superposed on the general 
flow of air mass in the upper atmosphere. A mathematical examination is 
necessary, in which the above problem must be dealt with as a transient 
phenomenon. 


§ 4. Tonospheric problems 


Magnetic storms are caused by the neutral ionized solar corpuscular 
stream, as already mentioned above. The solar stream is considered to be 
composed mainly of protons and electrons, though the presence of other 
elements of material is not yet ascertained. Geomagnetic disturbances of polar 
origin are attributed to the intense ionization in the ionospheric region over 
the auroral zone by precipitating protons, as already mentioned about. Ioniza- 
tion of electrons, which are expelled also from the equatorial current-ring, 
however, is not dealt with at present. Though electrons do not seem easily 
to penetrate into the level of the auroral height and contribute to the forma- 
tion of the current-system for polar magnetic disturbance field, they will cause 
some ionospheric change in the upper part of the ionospheric region. Precipitat- 
ing protons also lose their energy during passing through the ionospheric region 
by inellastic and ellastic collisions with atmospheric particles. The lost energy 
of impinging protons and electrons is absorbed by atmospheric particles, resulting 


Polar Magnetic Storms and Geomagnetic Bays 407 


in the ionization and excitation of atoms and molecules, or other change of 
physical condition, say temperature change. The above change in the physical 
condition of the upper atmosphere will be revealed as the ionospheric disturbance 
accompanying polar geomagnetic disturbance. The concrete exploration of the 
change in the physical condition of the ionospheric region is a problem in 
future. Recently, T. Nagata and T. Oguti showed a theoretical opinion on the 
mechanism of F'2 layer disturbance in the auroral zone [91]. 

The ionospheric region in the middle and low latitudes is a forbid- 
den region for extra-terrestrial charged particles. Hence, the ionospheric 
perturbation accompanying polar geomagnetic disturbance there will not be 
caused by direct bombardment of charged particles. From the characteristics 
of F2 layer disturbance suggested by some investigators [81, 83, 91], it may be 
allowed to suppose that F2 layer disturbance propagates from the high 
latitudes with a finite velocity, or F2 layer in lower latitudes suffers some 
change in its physical condition by some mechanism of diffusive character. 
Although F2 layer disturbance seems to be an associated phenomenon with 
magnetic storms, #2 layer disturbance is one of the most important problems 
for the ionospheric research. In the present stage of investigation, although 
various interesting and important characteristics of F2 layer disturbance over 
the world have been presented, the study of mechanism for the disturbance 
belongs to the future problem. D.F. Martyn proposed an interpretation of 
ionospheric storm in the middle and low latitudes in connexion with geomag- 
netic disturbance, after comprehensive study on geo-morphology of F2 layer 
disturbance [78 ]. 

In the above discussions, only charged particles are considered. The 
neutral ionized solar stream will contain also neutral particles. Neutral par- 
ticles are free from the influence of the earth’s magnetic field, and are able to 
impinge upon the upper atmosphere in lower latitudes. M.Sumi [116] showed 
a calculated result concerning the penetration of neutral particles under several 
assumptions and approximations. According to his result, hydrogen-like atoms 
of about 108cm/sec. in the initial veloctity penetrate into about 100km. above 
the earth’s surface and they produce thin ionized layers or clouds, the electron 
density of which is capable to be larger than that of E layer. This result 
seems favourite for the occurrence of sporadic E layer at the time of magnetic 
storms in the lower latitudes over the earth. Recent continuous observation 
of h/-f curve by means of a movie camera by Y. Nakata illustrates various 
important features of sporadic ionization at the time of magnetic disturbance. 
Y. Nakata, M. Kan and H. Uyeda [117] suggested an extremely high correlation 
between the ionospheric change in & region and the simultaneous geomagnetic 


variation. in the middle latitudes. 
Generally speaking, the physical mechanism of ionospheric disturbance 


over the world must be discussed in future. Further systematic study is 


408 N. FUKUSHIMA 


required both in experimental and theoretical operations. The comparison of 
world-wide ionospheric perturbation with simultaneous geomagnetic data for 
individual cases of magnetic storms will give many suggestive results for the 
mechanism of ionospheric disturbance. Recent world-wide distribution of 
ionospheric stations makes us possible to look over a general view of ionospheric 
perturbation over the world to some extent, but it is difficult now to know the 
simultaneous geomagnetic condition over the world, because the accumulation 
of a sufficient number of geomagnetic data over the world is not possible, 
especially with respect to the polar region data. On the other hand, iono- 
spheric observation was not yet distributed over the world during the Second 
International Polar Year, 1932-1933, in which fairly sufficient number of geo- 
magnetic observations were co-operated. It is strongly desired to compare the 
world-wide ionospheric disturbance with simultaneous geomagnetic data over the 
whole world in near future. 
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